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Abstract 



This is the first graduate course on elementary quantum mechanics in Inter- 
net written for the benefit of undergraduate and graduate students. It is a 
translation (with corrections) of the Romanian version of the course, which 
I did at the suggestion of several students from different countries. The top- 
ics included refer to the postulates of quantum mechanics, one-dimensional 
barriers and wells, angular momentum and spin, WKB method, harmonic 
oscillator, hydrogen atom, quantum scattering, and partial waves. 
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0. FORWARD 



The energy quanta occured in 1900 in the work of Max Planck (Nobel prize, 
1918) on the black body electromagnetic radiation. Planck's "quanta of 
light" have been used by Einstein (Nobel prize, 1921) to explain the pho- 
toelectric effect, but the first "quantization" of a quantity having units of 
action (the angular momentum) belongs to Niels Bohr (Nobel Prize, 1922). 
This opened the road to the universalization of quanta, since the action is 
the basic functional to describe any type of motion. However, only in the 
1920's the formalism of quantum mechanics has been developed in a system- 
atic manner. The remarkable works of that decade contributed in a decisive 
way to the rising of quantum mechanics at the level of fundamental theory 
of the universe, with successful technological applications. Moreover, it is 
quite probable that many of the cosmological misteries may be disentan- 
gled by means of various quantization procedures of the gravitational field, 
advancing our understanding of the origins of the universe. On the other 
hand, in recent years, there is a strong surge of activity in the information 
aspect of quantum mechanics. This aspect, which was generally ignored in 
the past, aims at a very attractive "quantum computer" technology. 

At the philosophical level, the famous paradoxes of quantum mechanics, 
which are perfect examples of the difficulties of 'quantum' thinking, are 
actively pursued ever since they have been first posed. Perhaps the most 
famous of them is the EPR paradox (Einstein, Podolsky, Rosen, 1935) on the 
existence of elements of physical reality, or in EPR words: "If, without in any 
way disturbing a system, we can predict with certainty (i.e., with probability 
equal to unity) the value of a physical quantity, then there exists an element 
of physical reality corresponding to this physical quantity. " Another famous 
paradox is that of Schrodinger's cat which is related to the fundamental 
quantum property of entanglement and the way we understand and detect 
it. What one should emphasize is that all these delicate points are the 
sourse of many interesting and innovative experiments (such as the so-called 
"teleportation" of quantum states) pushing up the technology. 

Here, I present eight elementary topics in nonrelativistic quantum me- 
chanics from a course in Spanish ( "castellano" ) on quantum mechanics that 
I taught in the Instituto de Fisica, Universidad de Guanajuato (IFUG), 
Leon, Mexico, during the semesters of 1998. 

Haret C. Rosu 
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1. THE QUANTUM POSTULATES 

The following six postulates can be considered as the basis for theory and 
experiment in quantum mechanics in its most used form, which is known as 
the Copenhagen interpretation. 

PI.- To any physical quantity L, which is well defined at the classical level, 
one can associate a hermitic operator L. 



P2.- To any stationary physical state in which a quantum system can be 

|2 
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found one can associate a (normalized) wavcfunction. (|| \\'j'2= 1) 



P3.- In (appropriate) experiments, the physical quantity L can take only the 
eigenvalues of L. Therefore the eigenvalues should be real, a condition 
which is fulfilled only by hermitic operators. 

P4.- What one measures is always the mean value L of the physical quantity 
(i.e., operator) L in a state tpn, which, theoretically speaking, is the 
corresponding diagonal matrix element 

{lpn \ L\ tpn) = L. 

P5.- The matrix elements of the operators corresponding to the cartesian 
coordinate and momentum, Xi and pk, when calculated with the wave- 
functions / and g satisfy the Hamilton equations of motion of classical 
mechanics in the form: 

J^{f\p^\9)=-{f\§-\9) 

l(/l£-d5) = (/lgl5), 

where H is the hamiltonian operator, whereas the derivatives with 
respect to operators are defined as at point 3 of this chapter. 

P6.- The operators pi and Xk have the following commutators: 



[xi, xj,]=0 
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n, = h/2'K = 1.0546 X 10-2^ erg.sec. 

1. - The correspondence between classical and quantum quantities 

This Ccin be done by substituting x^, with Xi pk- The function L is 
supposed to be analytic (i.e., it can be developed in Taylor series). If 
the L function does not contain mixed products x^pfe, the operator L 
is directly hermitic. 
Exemple: 

T = (EbD/2m ^f= {Elf)/2m. 

If L contains mixed products XiPi and higher powers of them, L is not 
hermitic, and in this case L is substituted by A, the hermitic part of 
L (A is an autoadjunct operator). 
Exemple: 

w{xi,pi) = J2iPiXi — ^ w = l/2Y^l{piXi + XiPi). 

In addition, one can see that wc have no time operator. In quan- 
tum mechanics, time is only a parameter that can be introduced in 
many ways. This is so because time does not depend on the canonical 
variables, merely the latter depend on time. 

2. - Probability in the discrete part of the spectrum 

If n is an eigenfunction of the operator L, then: 

L =< n \ L \ n >=< n | A„ | n >= Xn < n \ n >= 5nn^n = ^n- 
Moreover, one can prove that L = (A„) . 

If the function (j) is not an eigenfunction of L, one can make use of the 
expansion in the complete system of eigenfunctions of L to get: 

Ltpn = Ktpn, ^ = En OnV'n 

and combining these two relationships one gets: 

L(l> = En ^nanlpn- 
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In this way, one is able to calculate the matrix elements of the operator 
L: 

(0 I L I 0) = Y.n,mKn(^nK{'m \ u) = J^m I P ^m, 

telling us that the result of the experiment is with a probability 

I I • 

If the spectrum is discrete, according to P4 this means that | am P, 
that is the coefficients of the expansion in a complete set of cigcnfunc- 
tions, determine the probabilititics to observe the eigenvalue A„. 
If the spectrum is continuous, using the following definition 

<^(T)=/a(A)V'(r,A)dA, 

one can calculate the matrix elements in the continuous part of the 
spectrum 

(</> I L I </>) 

= / dr / a*(A)V'* (r, A)dA / /xa(/^)^/)(T, n)dn 
= J J a*a{^ji)^ji J iP*{t, X)tp{tau, fj,)dXdfj,dT 
— II ci*(A)a(jLt)jLt(5(A — /x)dAd/x 
= /a*(A)a(A)AdA 
= / I a(A) p AdA. 

In the continuous case, | a(A) p should be understood as the probabil- 
ity density for observing the eigenvalue A belonging to the continuous 
spectrum. Moreover, the following holds 

L = {<t>\L\<P). 

One usually says that {fi \ <1>) is the representation of | in the 
representation /i, where | /i) is an eigenvector of M. 
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3.- Definition of the derivate with respect to an operator 



dF(L) _ F{L+d)-F(L) 

4.- The operators of cartesian momenta 

Which is the exphcit form of pi, and pjj, if the arguments of the 
wavefunctions are the cartesian coordinates Xi ? 
Let us consider the following commutator: 

IPii^i ] ~ Pi^i Pi 
= PiXiXi XiPiXi -\- XipiXj, XiXipi 
= {piXi X'iP^^Xi -\- Xii^P'iXi XiPi^ 

— \pi J Xi\Xi + Xi \pi , Xi] 

= —ihxi — ihxi = —2ihxi. 
In general, the following holds: 

Pi^i Pi ~ TliflXi 

Then, for all analytic functions we have: 

Piip{x) - ipix)pi = -ih§^^. 

Now, let Pi(j) = f{xi,X2, xs) be the manner in which acts on (j){xi,X2,X3) = 
1. Then: 

Pii/j = —i^-§^ + /iV' aiid similar relationships hold for X2 and X3. 
Prom the commutator \pi,Pk] = it is easy to get V x / = and 
therefore /, = V^F. 

The most general form of pi is pi = —i^-^ + §§-, where F is an ar- 
bitrary function. The function F can be eliminated by the unitary 
transformaton = ex.p{^F). 
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P^ = UH-^n£: + §§:)U 
= expt^(-i?l^ + ^)expi'-^ 

leading to pi = — — ^ P = —ihV. 

5.- Calculation of the normalization constant 

Any wavefunction '(/'(x) G of variable x can be written in the form: 

that can be considered as the expansion of ijj in eigenfunction of the 
operator position (cartesian coordinate) xd{x — = ^{x — ^). Thus, 

I tpi^x) p is the probability density of the coordinate in the state ip{x). 
From here one gets the interpretation of the norm 

II tp(x) ||^= / I il}{x) p dx = 1. 

Intuitively, this relationship tells us that the system described by il>{x) 

should be encountered at a certain point on the real axis, although we 

can know only approximately the location. 

The eigenfunctions of the momentum operator are: 

—ih^ = Pii^, and by integrating one gets ip{xi) = Aexpt*'*^'. x andp 

have continuous spectra and therefore the normalization is performed 

by means of the Dirac delta function. 

Which is the explicit way of getting the normalization constant ? 

This is a matter of the following Fourier transforms: 

f{k) = J g{x) exp-*^^ dx, g{x) = ^ / f{k) exp^*^^ dk. 

It can also be obtained with the following procedure. Consider the 

unnormalized wavefunction of the free particle 

ipx 

4>p{x) = Aexp~ and the formula 

S{x-x') = ^ exp^*^(^-^') dx . 
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One can see that 

= J_oo ^ exp n A exp n dx 

= J-oo I ^ r exp n dx 

I 

I ^ l9 TOO ^^(P-P ) 

= 1'Kh\A |2 - p') 

and therefore the normaUzation constant is: 

A = ^ 

Moreover, the eigenfunctions of the momentum form a complete sys- 
tem (in the sense of the continuous case) for ah functions of the jC.'^ 
class. 

^ 7^ "'^P^ expT dp 

"'^P^ ^ vkh ^ exp-^ dx. 

These formulae provide the connection between the x and p represen- 
tations. 

6.- The momentum (p) representation 

The explicit form of the operators pi and Xk can be obtained either 
from the commutation relationships or through the usage of the kernels 
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1 1/1 —zpx ipx 

^ip^P) = U^xU = 2^ / exp^^i" xexp" dx 

1 —ipx p. i(3x 

The integral is of the form: M(A,A') = J U\X,x)MU{X' ,x)dx, and 
using xf = J x{x,^)f{^)d^, the action of x on a{p) G is: 

xa{p) = J x{p, P)a{P)dP 
= /(2k/exp^(-^^^exp^)da;)a(/?)d/3 
= 2^ i i exp a ^ exp a a{l3)dxdf3 
= ^II exp"^ ^ exp^ a{P)d^dp 

= //exp^ ^a{P)d^dp 
= -inj'-^S{P-p)d(3 = -ih'-^, 
where (5(/? — p) = J exp a d|. 

The momentum operator in the p representation is defined by the ker- 
nel: 

p{p,P) = WpU 

—ipx p. if3x 

= 2^1 exp~ i-i^Q^) exp— dx 

1 —ipx ipx 

= 2^iexp R /3exp fi dx = (iX(p - P) 
11 



leading to 'pa{p) = pa{p). 

It is worth noting that x and p, although hermitic operators for all 
f(x) € are not hermitic for their own eigenfunctions. 
If pa{p) = PoO-ip) and x = p = p\ then 



<a\px\a> — <a\xp\a >= —ih < a \ a > 

Po[< a\x\a> — <a\x\a>] = —ih < a \ a > 

Po[< a\x\a> — <a\x\a>]=0 

The left hand side is zero, whereas the right hand side is indefinite, 
which is a contradiction. 

7.- Schrodinger and Heisenberg representations 

The equations of motion given by P5 have different interpretations 
because in the expression -^{f \ L \ f) one can consider the temporal 
dependence as belonging either to the wavefunctions or operators, or 
both to wavefunctions and operators. We shall consider herein only 
the first two cases. 

• For an operator depending on time O = 0{t) we have: 

ri. - -dH. r.-M. 

^» ~ dxi ' ~ dpi 

[p, f] = pf - fp = -iT^mi 

[x, f] = xf - fx = -ifi§i 
and the Heisenberg equations of motion are easily obtained: 

Pi = ^\p,H], Xi = ^[x,H]. 
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• If the wavefunctions are time dependent one can still use pi = 
^\pi,H], because being a consequence of the commutation rela- 
tions it does not depend on representation 

it< f\Pi\9>=^ < f\\p,H]\g>. 

If now Pi and H do not depend on time, taking into account the 
hermiticity, one gets: 

= ^{pf,Hg) + ^{Hf,Pig) 

The latter relationship holds for any pair of functions f{x) and 
g{x) at the initial moment if each of them satisfies the equation 

This is the Schrodinger equation. It describes the system by 
means of time-independent operators and makes up the so-called 
Schrodinger representation. 

In both representations the temporal evolution of the system is char- 
acterized by the operator which can be obtained from Hamilton's 
function of classical mechanics. 

Exemple: H for a particle in a potential U{x-i,X2-,x^) we have: 
H = ^ + U{xi,X2,X3), which in the x representation is: 

H = -^Vl + U{xi,X2,X3). 

8.- The connection between the S and H representations 

P5 is correct in both Schrodinger's representation and Heisenberg's. 
This is why, the mean value of any observable coincides in the two 
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representations. Thus, there is a unitary transformation that can be 
used for passing from one to the other. Such a transformation is of 

the form = exp~?i~. In order to pass to the Schrodinger repre- 
sentation one should use the Heiscnberg transform = s^ff with / 
and L, whereas to pass to Heisenberg's representation the Schrodinger 
transform A = s"l"L.s with tp and A is of usage. One can obtain the 
Schrodinger equation as follows: since in the transformation ij^ = f 
the function / does not depend on time, we shall derivate the trans- 
formation with respect to time to get: 

t = #/ = i(exp^)/ = f ^exp-# / = ^Hs^f = ^H^.. 
Therefore: 

iU^ = H^. 

Next we get the Heisenberg equations: putting the Schrodinger trans- 
form in the form sAst = L and performing the derivatives with respect 
to time one gets Heisenberg's equation 

t = fA^"^ + = t#exp# A.t - i|Aexp-F H 

= ^{HsAs^ - sAs^H) = ^{HL - LH) = ^[H, L]. 
Thus, we have: 

Moreover, Heisenberg's equation can be written in the form: 
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L is known as an integral of motion, which, if ^ < | L | ^ >= 0, is 
characterized by the foUowing commutators: 

[H,L] = Q, [^,A] = 0. 

9.- Stationary states 

The states of a quantum system described by the eigenfunctions of H 
are called stationary states and the corresponding set of eigenvalues 
is known as the energy spectrum of the system. In such cases, the 
Schroedinger equation is: 

The solutions are of the form: ipn{x, t) = exp rf^ 4>n{x)- 

• The probability is the following: 

5{x) =1 ll^n{x,t) P = | exp^f^ (l)n{x) ^ 

= exp R exp R I 1^=1 (pn{x) p. 

Thus, the probability is constant in time. 

• In the stationary states, the mean value of any commutator of 
the form [H, A] is zero, where A is an arbitrary operator: 

< n I hA - Ah I n >=< n \ hA \ n> - <n \ Ah \ n > 
=< n I EnA \ n> — <n \ AEn \ n > 
= En <n \ A\n> -En < n | 1 | >= 0. 

• The virial theorem in quantum mechanics - if is a hamiltonian 
operator of a particle in the field U{r), using 

A = 1/2 J2i=i{PiXi - XiPi) one gets: 



15 



<ilj\[A,H]\ip>=0=<ip\AH-HA\ip> 

= Ei=i < I PiXiH - HpiXi I ip > 

= Td=i < V' I Xi\pi + Xi[H,pi\ \'tp>. 

Using several times the commutators and pi = —ifiVi, H = 
T + U{r), one can get: 

< V I [i, ^] I V' >= 

= -in{2 <^|f|'0>-<V'|^- Vi/M I tp >). 

This is the virial theorem. If the potential is U{r) = Uor^, then 
a form of the virial theorem similar to that in classical mechanics 
can be obtained with the only difference that it refers to mean 
values 

T = ffJ. 

• For a Hamiltonian H = -^V^ + U{r) and [r, H] = =^p, calcu- 
lating the matrix elements one finds: 

(Ek-E^) <n\f\k>='^ <n\p\k>. 

The uour(;lalivitilic pixjljability currenl density 
The following integral: 

/ I i'nix) P dx = 1, 

is the normalization of an eigenfunction of the discrete spectrum in 
the coordinate representation. It appears as a condition on the micro- 
scopic motion in a finite region of space. 

For the wavefunctions of the continuous spectrum one cannot 

give a direct probabilistic interpretation. 

Let us consider a given wavefunction (j) £ C'^, that we write as a linear 
combination of eigenfunctions of the continuum: 
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(t> = J a{X)^x{x)dx. 



One says that (p corresponds to an infinite motion. 
In many cases, the function a(A) is not zero only in a small neighbor- 
hood of a point X = Xq- In such a case, cp is known as a wavepacket. 
We shall calculate now the rate of change of the probability of finding 
the system in the volume 

P = In \i'{x,t) \^ dx = i^i;*{x,t)i){x,t)dx. 
Derivating the integral with respect to time leads to 

Using now the Schrodinger equation in the integral of the right hand 
side, one gets: 

Using the identity fV'^g — gV"^ f = div[{f)grad{g) — {g)grad{f)\ and 
also the Schrodinger equation in the form: 

and subtituting in the integral, one gets: 

f = I /n[V'(-l^vv*) - ri^^m^ 

= -In div^ii/jVip* - ilj*Vip)dx. 

By means of the divergence theorem, the volume integral can be trans- 
formed in a surface one leading to: 
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The quantity J{tp) = f^('0V'0* — ip*V'4>) is known as the probabihty 
density current, for which one can easily get the following continuity 
equation 



f^+div{J) = Q. 

• If il){x) = AR{x), where R{x) is a real function, then: J{tp) = 0. 

• For momentum eigenfunctions tp{x) = (27rfi)^/2 ^xp~i oi^c gets: 

-((W72^^P ' ^{2ifeexp a )) 

ift / 2ip \ p 

~ 2m V h{2iTh)'-^ ' ~ m(2iTny-i ' 

which shows that the probability density current does not depend 
on the coordinate. 

Operator of spatial transport 

If H is invariant at translations of arbitrary vector a, 



H{r + a) = H{r) , 

then there is an operator T(a) which is unitary T\d)H{f)T[d) = 
H{r + d). 

Commutativity of translations 

f (a)f (6) = f (6)f (a) = f (a + 5), 

implies that T is of the form T = exp'*^", where ^ = f • 
In the infinitesimal case: 

f{6d)Hf{Sd) «i (7 + ik6d)H{i - ikSd), 
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H{r) + i[K, H]6a = H{r) + {VH)6a. 



Moreover, [p, H] = 0, where p is an integral of the motion. The sistem 
of wavefunctions of the form '^{p,r) = (2^^)3/2 ^^P^ the unitary 

ipa. 

transformation leads to exp~ ip{f) = '^{r+a). The operator of spatial 

-^j. —ipa J. —iHt 

transport TT = exp n is the analog of = exp s , which is the 
operator of time 'transport' (shift). 

- Exemple: The 'crystal' (lattice) Hamiltonian 

If H is invariant for a discrete translation (for exemple, in a crystal 
lattice) H{f + a) = H{r), where a = J2i(Hni, rii G N and Oj are 
baricentric vectors, then: 

H{f+ a)'4){f+ a) = Etp{f+ a) = H(f)il^{f+ a). 

Consequently, tlj{f) and tlj{f+a) are wavefunctions for the same eigen- 
value of H. The relationship between t/^(r) and ip{f+a) can be saught 
for in the form ip{f+ a) = c{a)ip{f), where c(a) is a gxg matrix (g is 
the order of degeneration of level E). Two column matrices, c(a) and 
c{b) commute and therefore they are diagonalizable simultaneously. 
Moreover, for the diagonal elements, Cii{a)cii{b) = Cii{a + h) holds 
for i=l,2,....,g, having solutions of the type Cii{a) = exp'*^*". Thus, 
V'fc(^ = ?7fe(r) exp**^", where is a real arbitrary vector and the func- 
tion Uk{r) is periodic of period a, Uk{r + a) = Uk{r)- 
The assertion that the eigenfunctions of a periodic H of the lattice 
type H{f + a) = H{f) can be written V'A;(^) = U^if) expika^ where 
Uk{f+ a) = Uk{r) is known as Bloch's theorem. In the continuous 
case, Uk should be constant, because the constant is the only function 
periodic for any a. The vector p = hk is called quasimomentum (by 
analogy with the continuous case). The vector k is not determined 
univoquely, because one can add any vector g for which ga = 27rn, 
where n & N. 

The vector g can be written g = J2i=i bifni-, where rrii are integers and 
hi are given by 
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bi = 2tt^ 



for i ^ j ^ k. bi are the baricentric vectors of the lattice. 
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IP. Problems 

Problema 1.1: Let us consider two operators, A and B, which commutes 
by hypothesis. In this case, one can derive the following relationship: 

Solution 
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Defining an operator F(t), as a function of real variable t, of the form: 

then: ^ = Ae^^e^^ + e^^Be^^ = {A + e^*Be-^*)F{t). 
Applying now the formula [A,F{B)\ = [A,B\F (B), wc have 
[e^*, B] = i[yl.S]e^*, and therefore: e^*S = 5e^* + t[A, 5]e^* . 
Multiplying both sides of the latter equation by exp~"^* and substituting in 
the first equation, we get: 

d£ = (A + B + t[A,B])F{t). 

The operators A , B and [A,B] commutes by hypothesis. Thus, we can 
integrate the differential equation as if A + B and [A, B] would be scalar 
numbers. 
We shall have: 

F{t) = F(0)e(^+^)*+^/2[A,B]t2 _ 

Putting t = 0, one can see that -F(O) = 1 and therefore : 

^ ^(A+B)t+l/2[A,B]t\ 

Putting now t = 1, we get the final result. 

Problem 1.2: Calculate the commutator [X, D^;]. 
Solution 

The calculation is performed by applying the commutator to an arbitrary 
function ■(/'(r): 

[X, Z5.]V(r) = {x£ - Ixmr) = x£^{f) - ^[x^(r)] 

= x^ip{r) - ijj{r) - x^i]j{r) = -tl^if). 

Since this relationship is satisfied for any '4){r), one can conclude that [X, D^] = 
-1. 

Problem 1.3: Check that the trace of a matrix is invariant of changes of 
discrete orthonormalized bases. 

Solution 

The sum of the diagonal elements of a matrix representation of an operator 
A in an arbitrary basis does not depend on the basis. 
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This important property can be obtained by passing from an orthonormal- 
ized discrete basis \ Ui > to another orthonormaUzed discrete basis \ tk >. 
Wc liave: 

'Ei<Ui\A\ui>=J2i<Ui\ (Efe I tk >< tk\)A\ui> 

(where we have used the completeness relationship for the states tk)- The 
right hand side is: 

<Ui\tk><tk\A\ui >= <tk\A\ui><Ui\tk >, 

(the change of the order in the product of two scalar numbers is allowed). 

Thus, we can replace J2i I >< Ui \ with unity (i.e., the completeness 
relationship for the states | Ui >), in order to get finally: 

^ < Ui\ A \ Ui >=^ <tk\ A\tk > . 

i k 

Thus, we have proved the invar iancc property for matriceal traces. 

Problem 1.4: If for the hermitic operator N there are the hermitic oper- 
ators L and M such that : [M, N] = 0, [L, N] = 0, [M, L] 0, then the 
eigenfunctions of N are degenerate. 

Solution 

Let iIj{x; /jl, u) be the common eigenfunctions of M and N (since they com- 
mute they are simultaneous obscrvablcs) . Let 'ip{x;X,i') be the common 
eigenfunctions of L and (again, since they commute they are simulta- 
neous observables) . The Greek parameters denote the eigenvalues of the 
corresponding operators. Let us consider for simplicity sake that N has a 
discrete spectrum. Then: 

fi^) = 51 ai^'<P{x; n,iy) =^ K'tp{x; A, u) . 

V V 

We calculate now the matrix element < f\ML\f >: 

< f\]^L\f >= / ^ Hvttuip* {x\ iJL,v)^Wiip{x\\,v )dx . 

If all the eigenfunctions of iV are nondegenerate then < f\ML\f >= iii^aiyXi/b^. 
But the same result can be obtained if one calculates < f\LM\f > and the 
commutator would be zero. Thus, at least some of the eigenfunctions of N 
should be degenerate. 
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2. ONE DIMENSIONAL RECTANGULAR 
BARRIERS AND WELLS 

Regions of constant potential 

In the case of a rectangular potential, V{x) is a constant function 
V{x) = y in a certain region of the one-dimensional space. In such a 
region, the Schrodinger eq. can be written: 



^l;{x) + ^{E-V)i^ix)=0 (1) 



One can distinguish several cases: 

(i) E>V 

Let us introduce the positive constant k, defined by 

Then, the solution of eq. (1) can be written: 

iIj{x) = Ae'^"" + A!e-'^'' (3) 

where A and A' are complex constants. 

(ii) E <V 

This condition corresponds to segments of the real axis which would be 
prohibited to any particle from the viewpoint of classical mechanics. In this 
case, one introduces the positive constant q defined by: 

(4) 

n 

and the solution of (1) can be written: 

V'(x) = 5e«^ + B'e-«^ , (5) 

where B and B' are complex constants. 

(iii) E = V 

In this special case, il){x) is a linear function of x. 
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The behaviour of il){x) at a discontinuity of the potential 

One might think that at the point x = xi, where the potential V{x) is 
discontinuous, the wavefunction ip{x) behaves in a more strange way, maybe 
discontinuously for example. This is not so: iIj{x) and ^ are continuous, 
and only the second derivative is discontinuous at x = xi. 

General look to the calculations 

The procedure to determine the stationary states in rectangular poten- 
tials is the following: in all regions in which V{x) is constant we write ^{x) 
in any of the two forms (3) or (5) depending on application; next, we join 
smoothly these functions according to the continuity conditions for ip{x) and 
^ at the points where V{x) is discontinuous. 

Examination of several simple cases 

Let us make explicite calculations for some simple stationary states 
according to the proposed method. 



The step potential 



V(x) 



V 



II 



Fig. 2.1 
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a. £^ > Vb case; partial reflexion 
Let us put eq. (2) in the form: 



ki = (6) 



The solution of eq. (1) has the form of eq. (3) in the regions I{x < 0) 
and II{x > 0): 

In region I eq. (1) takes the form 

2mE 

a, 

and in the region II: 



<(a;) + ^^ip{x) = ^"(x) + k^^{x) = 



V'"!^) - 5[^o - E]<p{x) = r{x) - q'^i^ix) = 
n 

If we hmit ourselves to the case of an incident particle 'coming' from x = 
— oo, wc have to choose A2 = and we can determine the ratios A'^/Ai and 
A2/A1. The joining conditions give then: 

• ipi = ipii, at a: = : 

Ai + A[ = A2 (8) 

• V/ = V'//) at a; = : 

Aiiki - A[iki = A2ik2 (9) 

Substituting Ai and A[ from (8) in (9): 

. Ml_M (,o) 

A, = (11) 



25 



From the two expressions of the constant A2 in (10) and (11) one gets 



A ^ h - fca 

and from (11) it follows: 

^ - (IS) 

ip{x) is a superposition of two waves. The first (the Ai part) corresponds 
to an incident wave of momentum p = Tik\, propagating from the left to 
the right. The second (the A'l part) corresponds to a reflected particle of 
momentum —hki propagating in opposite direction. Since we have already 
chosen A'2 = 0, it follows that ipii{x) contains a single wave, which is as- 
sociated to a transmitted particle. (We will show later how it is possible 
by employing the concept of probability current to define the transmission 
coefficient T as well as the reflection coefficient R for the step potential). 
These coefficients give the probability that a particle coming from x = —00 
can pass through or get back from the step at x = 0. Thus, we obtain: 

whereas for T: 

Taking into account (12) and (13) one is led to: 

R = I - -z (16) 

{h + k2)^ 

^ 4k,k2 

{h+k2f ^ ^ 

It is easy to check that R + T = 1. It is thus sure that the particle will 
be either transmitted or reflected. Contrary to the predictions of classical 
mechanics, the incident particle has a nonzero probability of not going back. 

It is also easy to check using (6), (7) and (17), that if £^ » Vq then 
T ~ 1 : when the energy of the particle is sufficently big in comparison with 
the height of the step, everything happens as if the step does not exist for 
the particle. 
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Consider the following natural form of the solution in region I: 



i = -^(0*V<^-^V^*) (18) 
with Aie*^i^ and its conjugate A\e'~'^^'^^ : 

j = - — [(.4*e-*^)(^izA;ie*^)-(ylie^'=i^)(-^tzA;ie-*^)] 
2m 

m 

Now with ^e~**^i^ and its conjugate yl*e**^i^ one is led to: 

m 

In the following we wish to check the proportion of reflected current 
with respect to the incident current (or more exactly, we want to check the 
relative probability that the particle is returned back): 



Similarly, the proportion of transmission with respect to incidence (that 
is the probability that the particle is transmitted) is, taking now into account 
the solution in the region II: 

b. ii^ < Vo case; total reflection 
In this case we have: 



k, = (21) 



= VMp^ (22) 
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In the region I(x < 0), the solution of eq. (1) [written as ^p{x)" +k1'ip{x) = 0] 
has the form given in eq. (3): 

V-/ = ^le^*^!^ + A[e-''''^ , (23) 

whereas in the region II{x > 0), the same eq. (1) [now written as il){x)" — 
q2'tp{x) = 0] has the form of eq. (5): 

■07/ = S2e«2^ + S^e-«2=^ . (24) 

In order that the solution be kept finite when x — > +oo, it is necessary that: 

B2 = . (25) 

The joining condition at x = give now: 

• ipi = i^ii, at x = : 

Ai+^'i = S^ (26) 

• ip'j = ip'jj, at a; = : 

Aiiki - A[iki = --62^2 • (27) 

Substituting Ai and A'^ from (26) in (27) we get: 

. fc-l PS) 

Equating the expressions for the constant B2 from (28) and (29) leads to: 

A[ iki +q2 ki — iq2 



Ai iki — q2 ki+ iq2 ' 

so that from (29) we have: 

B!^ _ 2iki _ 2ki 
Ai iki — q2 ki — iq2 



(30) 



(31) 
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Therefore, the reflection coefficient R is: 




ki - iq2 |2 ^ kf + g| 
fci + iq2 kl + g| 



= 1 . 



(32) 



As in classical mechanics, the microparticle is always reflected (total re- 
flexion). However, there is an important difference, namely, because of the 
existence of the so-called evanescent wave e~'^^^, the particle has a nonzero 
probability to find itself in a spatial region which is classicaly forbidden. This 
probability decays exponentially with x and turns to be negligible when x 
overcome 1 /q2 corresponding to the evanescent wave. Notice also that A'-^/Ai 
is a complex quantity. A phase difference occurs as a consequence of the 
reflexion, which physically is due to the fact that the particle is slowed down 
when entering the region x > 0. There is no analog phenomenon for this in 
classical mechanics (but there is of course such an analog in optical physics) . 

Rectangular bcirrier 



V(x) 



v„ 



1 



II 



111 







1 



X 



Fig. 2.2 



a. £^ > Vo case; resonances 

Here we put eq. (2) in the form: 



y/2mE 

n 



(33) 



k2 = 



^2m{E - Vb) 



(34) 
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The solution of eq. (1) is as in eq. (3) in the regions I{x < 0), //(O < 
X < a) and III{x > a) : 



If we hmit ourselves to the case of an incident particle coming from 

X = — oo, we have to choose ^3 = 

ipi = i^ii, at a; = : 



^1 + a; = A2 + (35) 

• = i^'iii at X = : 

Aiiki - A[iki = A2ik2 - ^2^A;2 (36) 

• tpii = tpiii, at X = a : 

Ase''^^" + A'2e-'''^^ = A^e'^^" (37) 

• i^'ii = V'///) at x = a : 

A2ik2e'''^'' - A!2ik2e-^^^'' = A^ikie'''^'' . (38) 

The joining conditions at x = a give A2 and A2 as functions of A3, whereas 
those at X = give Ai and A[ as functions of A2 and A'2 (thus, as functions 
of A3). This procedure is shown in detail in the following. 
Substituting A2 from eq. (37) in (38) leads to: 

Substituting A2 from eq. (37) in (38) leads to: 



Substituting Ai from eq. (35) in (36) leads to: 



(40) 



_ A2{k2-ki)-A'2{k2 + ki) 



30 



Substituting A[ from eq. (35) in (36) gives: 

A2{k2 + h) - A'2{k2 - h) 

= 2k[ • (^^^ 

Now, substituting the eqs. (39) and (40) in (41), we have: 

A[ = i ^^l~jJ\ smk2a)e''^'^As . (43) 
Finally, substituting the eqs. (39) and (40) in (42) we get: 

Ai = [cos k2a - i^lj^ sin feaaJe^'^iMs . (44) 
2kik2 

A'^/Ai and A^/Ai [these ratios can be obtained by equating (43) and (44), 
and by separating, respectively, in eq. (44)] allow the calculation of the 
reflexion coefficient R as well as of the transmission one T. For this type of 
barrier, they are given by the following formulas: 

n-\A' /A |2 - (fcf - fc|)^ sin^ k2a . . 

T = |A3Mi|2 = 4^2^2 + (^2^^'2)2sin2fc2a" ^^^^ 

It is easy to see that they check R + T = 1. 
b. £^ < Vb case; the tunnel effect 

Now, let us take the eqs. (2) and (4): 



= ^ (47) 

= «ZI). (48) 

The solution of eq. (1) has the form given in eq. (3) in the regions I{x < 
0) and III{x > a), while in the region 77(0 < x < a) has the form of eq. (5): 

i^iii = A^e'^^"" + A'ge"^*^!^ . 

The joining conditions at x = and x = a allow the calculation of the 
transmission coefficient of the barrier. As a matter of fact, it is not necessary 
to repeat the calculation: merely, it is sufficient to replace k2 by —iq2 in the 
equation obtained in the first case of this section. 
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Bound states in rectangular well 
a. Well of finite depth 

V(x) 



Fig . 2.3 Finite rectangular well 

We first study the case 0<£'<Vb(-E>Vbis similar to tlie calculation 
in the previous section). 

For the exterior regions I, {x < 0) and III, {x > a) we employ eq. (4): 



For the central region II (0 < x < a) we use eq. (2): 



k=^. (50) 

The solution of eq. (1) has the form of eq. (5) in the exterior regions and 
of eq. (3) in the central region: 
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In the region (0 < x < a) eq. (1) has the form: 

'ZttiE 

tp{x)" + -^tjj{x) = ^{x)" + k^^{x) = (51) 
while in the exterior regions: 

- 1? - ^]'^(^) = V'(^)" - = . (52) 

Because i/j should be finite in the region I, we impose: 

B[ = . (53) 

The joining conditions give: 
'tpi = ipii, at a; = : 

Bi = A2 + A'^ (54) 

ip'j = i/^'jj, at X = : 

Biq = A2ik - A2ik (55) 

"0// = i^iiij at x = a : 

A2e''"' + ^e-^'^'^ = ^369" + B'^e-'i'' (56) 

i/jjj = ipjjj, at X = a : 

A2ike''"' - A'^ike-''"' = B^qe^" - B'^qe'i" (57) 

Substituting the constants A2 and A'2 from eq. (54) in eq. (55) we get 

^ ^ Bi{q-ik) 
^ -2ik 

_ Bi{q + ik) 

- 2ik ' ^^^^ 



respectively. 
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Substituting the constant A2 and the constant A2 from eq. (56) in eq. 
(57) we get 

B'^e-'^''{ik + q)+ Bse'^'^iik - q) + A'^e-''"' {-2ik) = 
2ikA2e''"' + B'^e-'i''{-ik + q)+B3Ei^{-ik-q) = 0, (59) 

respectively. 

Equating B'^ from eqs. (59) and taking into account the eqs (58) leads to 
Since ^(x) should be finite in region III as well, we require B3 = 0. Thus 



q + ik e 



-ika 



Because q and k depend on E, eq. (1) can be satisfied for some particular 
values of E. The condition that ■0(x) should be finite in all spatial regions 
imposes the quantization of the energy. Two cases are possible: 

(i) if 

= _e«fca ^ ^g2) 



q + ik 

equating in both sides the real and the imaginary parts, respectively, we 



have 



Putting 



one gets 



tan(^) = I . (63) 



l2mVo 



ko = d—^ = ^k^ + q' (64) 



l + tan2(^) = ^ = (^)^ (65) 



cos2(^) ^ 2 ^ P ^ k 

Eq. (63) is therefore equivalent to the system of eqs. 



tan(y) > (66) 



34 



The energy levels are determined by the intersection of a straight line of 
slope ^ with the first set of dashed cosinusoides in fig. 2.4. Thus, we get 
a certain number of energy levels whose wavefunctions are even. This fact 
becomes clearer if we substitute (62) in (58) and (60). It is easy to check 
that i?3 = Bi and A2 = A'2 leading to ip{~x) = ^{x). 



(ii) if 



a similar calculation gives 



q — ik 
q + ik 



= e 



ika 



(67) 



|sin(^)| 



k_ 

ko 



,ka, 
tan(— ) < 



(68) 



The energy levels are in this case determined by the intersection of the 
same straight line with the second set of dashed cosinusoides in fig. 2.4. 
The obtained levels are interlaced with those found in the case (i). One can 
easily show that the corresponding wavefunctions are odd. 




2Jl/a 3Jl/a ^o^^/a 



Fig. 2.4 



b. Well of infinite depth 
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In this case it is convenient to put V{x) equal to zero for < x < a and 
equal to infinity for the rest of the real axis. Putting 



fc=,/^, (69) 



ip{x) should be zero outside the interval [0, a] and continuous at x = and 

X = a. 

For < X < a: 

V'(x) = Ae''''^ + A'e-''''' . (70) 
Since V'(O) = 0, one can infer that A' = —A, leading to: 

tp{x) = 2i^sin(A;x) . (71) 

Moreover, 'tp{a) =0 and therefore 

k = ^, (72) 
a 

where n is an arbitrary positive integer. If wc normalize the function (71), 
taking into account (72), then we obtain the stationary wavefunctions 

V'n(x) = y-sin(— -) (73) 

with the energies 

The quantization of the energy levels is extremely simple in this case. The 
stationary energies are proportional with the natural numbers squared. 



2P. Problems 

Problem 2.1: The attractive 5 potential 

Suppose we have a potential of the form: 



V{x) = -Vod{xy, Vb > 0; X G SR. 
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The corresponding wavefunction iIj{x) is assumed continuous. 

a) Obtain the bound states {E < 0), if they exist, locaUzed in this type 
of potential. 

b) Calculate the dispersion of a plane wave falling on the S potential and 
obtain the reflexion coefficient 



R 



'Jrefl\_ 

2 

tnc\ 



\x=0 



where iprefh fpinc arc the reflected and incoming waves, respectively. 
Suggestion: To determine the behavior of 'ip{x) in x=0, it is better to proceed 
by integrating the Schrodinger equation in the interval {—£,+£), and then 
to apply the limit e ^ 0. 

Solution, a) The Schrodinger eq. is: 



V'(^) + -^(^ + VoSix))ij{x) = . (75) 



Far from the origin we have a differential eq. of the form 



(P . 2mE , , , 

-^nx) = — ^nx)- (76) 



Consequently, the wavefunctions are of the form 

i;ix) = Ae-'i'' + Bel"" for x > and x < 0, (77) 



where q = y —2mE /fi^ G 5?. Since IV'P should be integrable , we cannot 
accept that a part of it grows exponentially. Moreover, the wavefunction 
should be continuous at the origin. With these conditions, we have 

V'(x) = ^e«^; (.T < 0), 

il^{x) = yle"^^; (x>0). (78) 
Integrating the Schrodinger eq. between —e and +£, we get 

- ^Wi^) - V''(-£)] - l^oV'(O) = E ip{x)dx « 2eE^{0) (79) 

Introducing now the result (78) and taking into account the limit e ^ 0,we 
have 

-—{-qA-qA)-VoA = 0, (80) 
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ov E = —m,{VQ /2h'^) [ — J- in units of Clearly, there is a single discrete 

energy. The normalization constant is found to be y4 = ^ mVo/h'^. The 

wavefunction of the bound state will be ipo = Ae^°l^l/^, where Vq is in ^ 
units. 



b) Take now the wavefunction of a plane wave 

2mE 

n 



It moves from the loft to the right and is reflected by the potential. If B and 
C are the amplitudes of the reflected and transmitted waves, respectively, 
then we have 

^{x) = ^e^^^ + Se-''^^; < 0), 

ip{x) = Ce^*^^; (x>0). (82) 

The joining conditions and the relationship ip'i^) — = cu 

/ = 2mVo/h'^ lead to 

A + B = C B = -—^—-A, 



ik{C-A + B) = -fC C = ^ , A (83) 

The reflection coefficient will be 



f + 2ik 
lik 
f + 2ik^ 



|V'mcP'"=° \A\^ m^Vi + h'k^- ^ ' 

If the potential is very strong (Vq — ^ oo), one can see that — > 1, i.e., the 
wave is totally reflected. 

The transmission coefficient, on the other hand, will be 

IV.ncP '"=° W mWS + n^k^- ^^^^ 

Again, if the potential is very strong (Vq oo) then T 0,i.e., the trans- 
mitted wave fades rapidly on the other side of the potential. 

In addition, R + T = 1 as expected, which is a check of the calculation. 
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Problem 2.2: Particle in a ID potential well of finite depth 

Solve the ID Schrodinger eq. for a finite depth potential well given by 

J — Vb daca \x\ < a 



Vix) 







daca \x\ > a . 



Consider only the bound spectrum (£' < 0) . 



+a 



E 



-Vo 



Fig. 2.5 

Solution. 

a) The wave function for |.t] < o and \x\ > a. 
The corresponding Schrodinger eq. is 

(x) + V(x)ijj(x) = Eip(x) . 

2m 

Defining 

2_ 2mE 2m(E + Vb) 

^ " " " ¥ ' ^^^^ 
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we get: 

1) for x<-a: 4j'i{x) - q'^^i = 0, = Aie'^'' + Bie''^''; 

2) for — a < X < a : ■i/'2 (2^) + ^^^2 = 0, ■02 = ^2 cos{kx) + B2 sm{kx); 

3) for X > a : 03(x) - q^^Pa = 0, ^3 = ^se''^ + ^ae"'"^- 

b) Formulation of the boundary conditions. 

The normahzation of the bound states requires solutions going to zero at 
infinity. This means Bi = = 0. Moreover, ijj^x) should be continuously 
differentiable. All particular solutions are fixed in such a way that ^ and 
are continuous for that value of x corresponding to the boundary between 
the interior and the outside regions. The second derivative tp" displays the 
discontinuity the 'box' potential imposes. Thus we are led to: 

ipii-a) = ip2{-a), ip2{a) = ipsia), 

V-K-a) = V2(-«), V2(«) = V'aC")- (88) 

c) The eigenvalue equations. 

From (88) we get four linear and homogeneous eqs for the coefficients 
Ai, A2, B2 and ^3: 

Aie~'^"' = A2 cos{ka) — B2 sin{ka) , 

qAier''"- = A2ksm{ka) + B2kcos{ka), 

B^e^'^"- = A2Cos{ka) + B2sm{ka), 

-gSae"*" = - A2k sm{ka) + B2k cos{ka) . (89) 

Adding and subtracting, one gets a system of eqs. which is easier to solve: 

(Ai + S3)e-«« = 2A2Cos{ka) 

g(^i + S3)e-«" = 2A2ksm{ka) 

(Ai-53)e"«" = -2S2sin(A;a) 

g(^i-B3)e-«" = 2B2kcos{ka). (90) 

Assuming + S3 7^ and ^2 7^ 0, the first two eqs give 

q = k tan(fea) , (91) 

which inserted in the last two eqs gives 

Ai = S3; B2 = 0. (92) 
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The result is the symmetric solution ip{x) = ip{—x), also called of positive 
parity. 

A similar calculation for Ai — ^ and 7^ leads to 



q = —kcot{ka) 



y 



^1 = -S.s 



^2=0. 



(93) 



The obtained wavefunction is antisymmetric, corresponding to a negative 
parity 

d) Quantitative solution of the eigenvalue problem. 
The equation connecting q and k, already obtained previously, gives the 
condition to get the eigenvalues. Using the notation 



from the definition (87) we get 



2 , 2 SmVba 2 



(94) 



(95) 



On the other hand, using (91) and (93) we get the equations 

?? = ^tan(0, r/=-Ccot(^). 

Thus, the sought energy eigenvalues can be obtained from the intersections 
of these two curves with the circle defined by (95) in the plane ^-r/ (see fig. 
2.6). 



11 = 5 tan 4 





Fig. 2.6 



There is at least one solution for arbitrary values of the parameter Vq, 
in the positive parity case, because the tangent function passes through the 
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origin. For the negative parity, the radius of the circle should be greater than 
a certain lower bound for the two curves to intersect. Thus, the potential 
should have a certain depth related to a given spatial scale a and a given 
mass scale m, to allow for negative parity solutions. The number of energy 
levels grows with Vq, a, and m. For the case in which mVo? oo, the 
intersections are obtained from 



tan(A;a) = oo 
-cot(A;a) = oo 



2n 



ka = 

2 

ka = mr, 



-TT, 



where n = 1, 2, 3, . . .; by combining the previous relations 

k{2a) = mr. 
For the energy spectrum this fact means that 

" 2m^2a^ " 



(96) 



(97) 



(98) 



Widening the well and/or the mass of the particle m, the diference between 
two neighbour eigenvalues will decrease. The lowest level (n = 1) is not 
localized at —Vq, but slightly upper. This 'small' difference is called zero 
point energy. 

e) The forms of the wavefunctions are shown in fig. 2.7. 





Fig. 2.7: Shapes of wave functions 
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Problem 2.3: Particle in ID rectangular well of infinite depth 

Solve the ID Schrodinger eq. for a particle in a potential well of infinite 
depth as given by: 



V{x) 



iov x' < X < x' + 2a 

oo for x' > X o a; > x' + 2a. 



The solution in its general form is 

ip(x) = Asm{kx) + Bcos{kx) , (99) 

where 



l2mE 

Since ijj should fulfill ip{x') = 'ip{x' + 2a) = 0, wc get: 



k = \l^. (100) 



A sm{kx') + B cos{kx') = (101) 
A sm[k{x' + 2a)] + B cos[k{x' + 2a)] = . (102) 

Multiplying (101) by sin[A;(x' + 2a)] and (102) by sin(A;x') and next subtract- 
ing the latter result from the first we get: 

B[ cos{kx') sm[k{x' + 2a)] - cos[A;(x' + 2a)] sin(A;x') ] = , (103) 

and by means of a trigonometric identity: 

Ssin(2aA;)=0 (104) 

Multiplying (101) by cos[A;(x' + 2a)] and subtracting (102) multiplied by 
cos(A;x') leads to: 

A[ sm{kx') cos[k{x' + 2a)] - sm[k{x' + 2a)] cos{kx') ] = , (105) 

and by means of the same trigonometric identity: 

^sin[A;(-2aA;)] = ^sin[A;(2aA;)] = . (106) 

Since we do not take into account the trivial solution V' = 0, using (104) 
and (106) one has sin(2aA;) = that takes place only if 2aA; = mr, with n an 
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integer. Accordingly, k = mr/2a and since = 2mE/h^ then it comes out 
that the eigenvalues are given by the following expression: 

E = . 107 

The energy is quantized because only for each fe„ = n7r/2a one gets a well- 
defined energy £'„ = [n? /2m][TTh/2a]'^ . 
The general form of the solution is: 

.rnrx. „ ,mrx. 
V'n = ^sm(^)+5cos(^), (108) 

and it can be normalized 

1= / 'il;il;*dx = a{A^ + B^), (109) 

Jx' 

wherefrom: 

A = ±^l/a-B2 . (110) 
Substituting this value of A in (101) one gets: 

1 . mrx' 

and plugging B in (110) we get 

, 1 .nvrx', /,,„x 
^ = ±^cos( ). (112) 

Using the upper signs for A and B, by substituting their values in (108) we 
obtain: 

V:.„ = ^sin(^)(x-x') . (113) 
Using the lower signs for A and B, one gets 
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3. ANGULAR MOMENTUM AND SPIN 



Introduction 

It is known from Classical Mechanics that the angular momentum 1 for 
macroscopic particles is given by 

l = rxp, (1) 

where r and p are the radius vector and the linear momentum, respectively. 

However, in Quantum Mechanics, one can find operators of angular mo- 
mentum type (OOAMT), which are not compulsory expressed only in terms 
of the coordinate xj and the momentum pk and acting only on the eigenfunc- 
tions in the x representation. Consequently, it is very important to settle 
first of all general commutation relations for the OOAMT components. 

In Quantum Mechanics 1 is expressed by the operator 

1 = -ihr X V, (2) 
whose components are operators satisfying the following commutation rules 

[^xj^j/] ~ i^zi ^yi^zl ~ i^xi ^zi^xl — ily (3) 

Moreover, each of the components commutes with the square of the angular 
momentum, i.e. 

f = ll + 1^ + 11 [kA=^, i = 1,2,3. (4) 

These relations, besides being correct for the angular momentum, are ful- 
filled for the important OOAMT class of spin operators, which miss exact 
analogs in classical mechanics. 

These commutation relations are fundamental for getting the spectra of the 
aforementioned operators as well as for their differential representations. 

The angular momentum 

For an arbitrary point of a fixed space (FS), one can introduce a function 
'ip{x, y, z), for which let's consider two cartesian systems E and E', where E' 
is obtained by the rotation of the z axis of E. 
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In the general case, a FS refers to a coordinate system, which is different 
of T, and S'. 

Now, let's compare the values of ip at two points of the FS with the same 
coordinates (x,y,z) in S and S', which is equivalent to the vectorial rotation 

i^(.x',y',z') = R'ilj{x,y,z) (5) 

where ii is a rotation matrix in 



(6) 



Then 

RiIj{x, y, z) = ip{x cos (j) — y smcp, x sin (/) + y cos (p, z). (7) 

On the other hand, it is important to recall that the wavefunctions are 
frame independent and that the transformation at rotations within the FS 
is achieved by means of unitary operators. Thus, to determine the form 
of the unitary operator U^{(l)) that passes ip to ip' , one usually considers 
an infinitesimal rotation d(f), keeping only the linear terms in d<j) when one 
expands ijj' in Taylor series in the neighborhood of x 

i'ix' ,y' ,z) PS + + 

~ i'{x,y,z) + d(j)(y^ - x^ 



dx dy , 

PS {l-id(j)l;,)Tp{x,y,z), (8) 

where we have used the notation^ 

lz = 'h'^{xpy-ypx). (9) 

As one will see later, this corresponds to the projection operator onto z of 
the angular momentum according to the definition (2) unless the factor . 
In this way, the rotations of finite angle (j) can be represented as exponentials 
of the form 

V'(x',y',z) =e*'-*^(x,y,z), (10) 

where 

= e''-<^. (11) 
^The proof of (8) is displayed as problem 3.1 
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In order to reassert the eoneept of rotation, we will consider it in a more 

general approach with the help of the vectorial operator A acting on ^, 
assuming that A^, Ay, A^ have the same form in S and S', that is, the 

mean values of A as calculated in S and S' should be equal when they are 
seen from the FS 



r{r){Aj + Ayj' + A,k')r{r) dr 
= j r{r){A:ci + Ayj + A,k)r{^ df, (12) 

where 

i' = icos (j) + jsincf), j' = z sin ^ + jcos 0, k' = k. (13) 
Thus, by combining (10), (12) and (13) we get 

e'^^'^A^e-'^^'^ = A^cos(t)- Ay sin (f>, 
e'^^'I'Aye-'^^'l' = A^sincl)- Ay cos 4>, 
e'^^'^A^e-'^^'^ = A^. (14) 

Again, considering infinitesimal rotations and expanding the left hand 
sides in (14), one can determine the commutation relations of A^, Ay and 
Az with Iz 

[lz,Ar^] = iAy, [l^jAy] = —lA^, [l^jA^] = 0, (15) 
and similarly for Ix and ly. 

The basic conditions to obtain these commutation relations are 
•k The eigenfunctions transform as in (7) when S — > S'. 

★ The components Ax, Ay, A^ have the same form in S and S'. 

★ The kcts corresponding to the mean values of in S and E' coincide 
(are the same) for a FS observer. 

One can also use another representation in which 'ip{x,y,z) does not 
change when S ^ S' and the vectorial operators transform as ordinary 
vectors. In order to pass to such a representation when we rotate by 
around z one makes use of the operator U {(f)) , that is 
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e'^^W{x,y,z) = ij{x,y,z), (16) 

and therefore 

^-ilz'l>Xe'^'='f' = A. (17) 
Using the relationships (14) we obtain 

= i^cos0 + ij/Sin(^ = e-^'-'^i^e*'-'^, 
A'y = -i^sin(/. + ij^cos(^ = e-^^-'^i^e^'-'^, 

= e-'^^'^A^e'^^^. (18) 

Since the transformations of the new representation are performed by 
means of unitary operators, the commutation relations do not change. 

Remarks 

★ The operator is invariant at rotations, that is 



e 



* It follows that 

[li,A^]=0. (20) 

★ If the Hamiltonian operator is of the form 

H=^p' + U{\r\), (21) 

then it remains invariant under rotations in any axis passing through 
the coordinate origin 

[ii,H]=0, (22) 
where li are integrals of the motion. 

Definition 

If Ai are the components of a vectorial operator acting on a wavcfunction 
depending only on the coordinates and if there are operators li that satisfy 
the following commutation relations 
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then li are known as the components of the angular momentum operator 
and we can infer from (20) and (23) that 

[k,P] = 0. (24) 

Consequently the three operatorial components associated to the com- 
ponents of a classical angular momentum satisfy commutation relations of 
the type (23), (24). Moreover, one can prove that these relations lead to 
specific geometric properties of the rotations in a 3D euclidean space. This 
takes place if we adopt a more general point of view by defining an angular 
momentum operator J (wc shall not use the hat symbol for simplicity of 
writing) as any set of three observables Jx, Jy §i Jz which fulfill the com- 
mutation relations 

[Ji, Jj] = iSijkJk- (25) 
Moreover, let us introduce the operator 

J2 = jj + J2 + Jf, (26) 

the scalar square of the angular momentum J. This operator is hermitic 
because Jx, Jy and Jz are hermitic and it is easy to show that commutes 
with the three components of J 

[j2,Ji] = 0. (27) 

Since commutes with each of the components it follows that there is 
a complete system of eigenfunctions, i.e. 

J^V'tm = /(7^)V'7M' Ji'^'Y^^ = 9{l^)'^7ti, (28) 

where, as it will be shown in the following, the eigenfunctions depend on two 
subindices, which will be determined together with the form of the functions 
7(7) and g{n). The operators Jj and Jfe {i 7^ k) do not commute, i.e. they do 
not have common eigenfunctions. For physical and mathematical reasons, 
we are interested to determine the common eigenfunctions of and Jz, 
that is, we shall take i = z in (28). 

Instead of using the components Jx and Jy of the angular momentum J, 
it is more convenient to work with the following linear combinations 

J^ — Jx ~t~ ^*^y, J— ~ Jx iJy- (29) 
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Contrary to the operators a and a' of the harmonic oscillator (see chapter 
5), these operators are not hermitic, they are only adjunct to each other. 



The following properties are easy to prove 

[J„J±]=±J±, [J+,J_] = 2J,, (30) 

[J^J+] = [J^J_] = [J^J,]=o. (3i) 

JziJ±i^'yn) = {J±Jz + [Jz, J±]}^J^^ = ± l){J±il;j^). (32) 



Therefore J^ijj^^^ are eigenfunctions of Jz corresponding to the eigenval- 
ues ;U ± 1 , that is these functions are identical up to the constant factors 
and j3n (to be determined) 

J-tpji^ = (33) 

On the other hand 

Therefore, taking a phase of the type e*'* (where a is real) for the function 
tp^n one can put real and equal to /3^, which means 

J+ipj,ij.~i = ajji'^^^, J-ipjij. = atJ,ipj,ij,-i, (35) 

and therefore 

7 = (V'7M' [Jx +Jy+ -^zIV'tm) = + a + b, 

b = {tp'yn,Jytpju.) = {Jy'ipju.,Jy'ipju.) >0. (36) 

The normalization constant cannot be negative. This implies 

7 > m', (37) 

for a fixed 7; thus, n has both superior and inferior limits (it takes values 
in a finite interval). 

Let A and A be these limits, respectively, for a given 7 

J+V7A = 0, J-V7A = 0. (38) 
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Using the following operatorial identities 

J_J+ = - + J, = J2 - J,(J, - 1), 

J+J_ = - + J, = J2 - J^(J, + 1), (39) 

acting on as well as on one gets 

7 - - A 
7 - + A 
(A-A + 1)(A + A) 

In addition, 

A> A^A = -A = J^7 = J(J + 1). (41) 

For a given 7 the projection ji of the momentum takes 2J + 1 values that 
differ by unity, from J to — J. Therefore, the difference A — A = 2 J should 
be an integer and consequently the eigenvalues of that are labelled by m 
are integer 

m = k, k = 0,±l,±2, ... , (42) 

or half- integer 

m = k + ^, k = 0,±l,±2, ... . (43) 

A state having a given 7 = J( J + 1) presents a degeneration of order 
g = 2J + 1 with regard to the eigenvalues m (this is so because Jj, 
commute with but do not commute between themselves. 

By a "state of angular momentum J" one usually understands a state 
of 7 = J( J + 1) having the maximum projection of its momentum, i.e. J. 
Quite used notations for angular momentum states are ipjm and the Dirac 
ket one \jm). 

Let us now obtain the matrix elements of Jx, Jy in the representation in 
which and are diagonal. In this case, one obtains from (35) and (39) 
the following relations 

J(J + l)-(m- l)"-(m-l) = a^, 

am = ^(J + m)(J-m + l). (44) 



= 0, 
= 0, 

= 0. (40) 
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Combining (44) and (35) leads to 

J+i^jm-i = ^J{J + m){J -m + l)V'jm • (45) 
It follows that the matrix element of J+ is 



- 1) = \J{J + m){J -m + l)5„m, (46) 

and analogously 



{in\J^\jm) = {J + m){J - m + l)6nm-i ■ (47) 
Finally, from the definitions (29) for J+, J_ one easily gets 

{jm\Jx\jm-l) = ^\/iJ + m){J -m + 1), 

{jm\Jy\jm - I) = ^^{J + m){J -m + l) . (48) 

Partial conclusions 

a Properties of the eigenvalues of J and 

If j{j + 1)%'^ and mh are eigenvalues of J and Jz associated to the 
eigenvectors \kjm), then j and m satisfy the inequality 

—j<m< j. 

P Properties of the vector J-\kjm) 

Let \kjm) be an eigenvector of and Jz with the eigenvalues j{j+\)h? 
and mh 

— (i) If m, = — j, then J-\kj — j) = 0. 

— (ii) If m > —j, then J^\kjm) is a nonzero eigenvector of and 
Jz with the eigenvalues j(j + l)?i^ and (m — l)h. 

7 Properties of the vector J+ 1 A;jm) 

Let \kjm) be a (ket) eigenvector of and J2 for the eigenvalues 
j{j + 1)^^ and mh 

★ If m = J, then Jj^\kjm) = 0. 

* If m < J, then Jj^\kjm) is a nonzero eigenvector of and Jz 
with the eigenvalues j{j + 1)^^ and {m + l)h 
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S Consequences of the previous properties 



Jz\kjm) = mhlkjm), 

J+\kjm) = mfi^ j{j + 1) — m{m + l)\kjm, + 1), 
J-\kjm) = mh\J j{j + 1) — m{m — l)\kjm + 1). 

Applications of the orbital angular momentum 

Until now we have considered those properties of the angular momentum 
that could be derived only from the commutation relations. Let us go back 
to the orbital momentum 1 of a particle without intrinsic rotation and let 
us examine how one can apply the theory of the previous section in the 
important particular case 

[ii,Pj] = ieijkPk- (49) 

First, Iz and pj have a common system of eigenfunctions. On the other 
hand, the Hamiltonian of a free particle 




being the square of a vectorial operator has a complete system of eigenfunc- 
tions with and Iz- In addition, this implies that the free particle can be 
found in a state of well-defined E, I, and m. 



Eigenvalues and eigenfunctions of P and Iz 

It is more convenient to work in spherical coordinates because, as we will 
see, various angular momentum operators act only on the angle variables 
9, (p and not on r. Thus, instead of describing r by its cartesian compo- 
nents X, y, z we determine the arbitrary point M of vector radius r by the 
spherical 3D coordinates 

x = rcos0sin^, y = r sm(j)sm9, z = rcos6, (50) 

where 

r > 0, 0<e <7T, 0<(p<2'K. 
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Let $(r, 6, (p) and $'(r, 6, (p) be the wavefunctions of a particle in S and 
S', respectively, in which the infinitesimal rotation is given by 6a around 
the z axis 



^{r,e,(P) + 6a—, (51) 



or 

$'(r, 9, = (1 + ilja)^r, 6, 4>). (52) 

It follows that 

For an inifinitesimal rotation in x 

= {l + iya)Hr,9,(l>), (54) 

but in this rotation 

z' = z + ySa; z' = z + ySa; x' = x (55) 

and from (50) one gets 

r cos (9 + d9) = r cos ^ + r sin sin ^(5a, 
r sin sin(^ + d^) = rsin^sin^ + rsin0sin<^ — rcos^fe, (56) 

i.e. 

d9 

sin 9d9 = sin 9 sin (p6a ^ — = — sin (f), (57) 

da 

and 

cos sin d0 + sin cos = — cos^^a, 

cos sin = — cos ^ — cos ^ sin 0— . (58) 

da da 

Substituting (57) in (56) leads to 

^ = -cot6'cos^ . (59) 
da 
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With (56) and (58) substituted in (51) and comparing the right hand sides 
of (51) one gets 

4 = i [sin (p— + cot 9 cos j ■ (60) 



For the rotation in y, the result is similar 



d d 
L = z ( - cos 0— + cot 9 sin (j)— 
^ o9 0(p 



(61) 



Using Ix, ly one can also obtain l±, r 



l± 



n 



exp 



+ V + 

1 ^2 1 d / d 
sin^ 9d^'^'^'^d9V''^ 89 



(62) 



Prom (62) one can see that I is identical up to a constant to the angular 
part of the Laplace operator at a fixed radius 



1 d f ^df 



dr V dr 



1 d / df\ 1 
^89 V'^ 89) ^ sin^ 9 d(j)'^ 



(63) 



The eigenfunctions of Iz 



d(t> ' 



27r 



(64) 



Hermiticity conditions of 

/•27r ^ / /•27r \ * 

rhgd<i>=(^j^ gnjdcpj +/*5(27r)- 7*3(0). 
It follows that L is hermitic in the class of functions for which 



(65) 



/*5(27r) = rg{0). 



(66) 
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The eigenfunctions ^rn of Iz belong to the integrable class (0, 27r) and 
they fulfill (66), as it happens for any function that can be expanded in 

k 

Ficf>) = ^afee^^^ k = 0,±l,±2, ... , 

k 

G{cj>) = Y.bke''"^, A: = ±l/2,±3/2,±5/2 ... , (67) 

with k only integers or half- integers, but not for combinations of F{(/)) and 
G((p). The correct choice of m is based on the common eigenfunctions of Iz 
and P. 

Spherical hcirmonics 

In the {?} representation, the eigenfunctions associated to the eigenvalues 
Z(Z + 1)^^ of 1^ and mh of Iz are solutions of the partial differential equations 

d 

—i-^'il){r., 9 .,(!)) = mfi'il^{r,6,(j)). (68) 

Taking into account that the general results presented above can be 
applied to the orbital momentum, we infer that / can be an integer or half- 
integer and that, for fixed I, m can only take the values 

-/ -I- 1, ... ,1-1,1. 

In (68), r is not present in the differential operator, so that it can be 
considered as a parameter. Thus, considering only the dependence on 6, (f) 
of Tp, one uses the notation Yim{d,(t>) for these common eigenfunctions of 1^ 
and Iz, corresponding to the eigenvalues l{l + 1)11^, mfi. They are known as 
spherical harmonics 

\^YUe,(i)) = i{i + i)n^Yi^{e,4>), 

lzYim{9,cp) = mnYi^{e,<P). (69) 

For more rigorousness, one should introduce one more index in order to 
distinguish among the various solutions of (69) corresponding to the same 
{I, m) pairs for particles with spin. If the spin is not taken into account, 
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these equations have a unique solution (up to a constant factor) for each 
ahowed pair of {l,m); this is so because the subindices l,m are sufficient 
in this context. The solutions YimiO, </>) have been found by the method of 
the separation of variables in spherical variables (see also the chapter The 
hydrogen atom) 

^lm{r,e,4>) = f{r)i;im{e,4>), (70) 

where /(r) is a function of r, which looks as an integration constant from 
the viewpoint of the partial differential equations in (68). The fact that /(r) 
is arbitrary proves that L"^ and Iz do not form a complete set of observables^ 
in the space e^.^ of functions of r {r,9,(j)). 

In order to normalize 'i/'/m(^) 4'): it is convenient to normalize ¥[^{0, (j)) 
and /(r) separately 

del) r sin e\ijimi9,<l))\^de = 1, (71) 

oo 

rV(r)|'dr = 1. (72) 



^0 



The values of the pair (/, m) 

(a): l,m should be integers 

Using Iz = ^-ir, we can write (69) as follows 



-^YUd,cl,)=mnYUB,ct>). (73) 

I 0(p 

Thus, 

y,„(0,0) =F;„(0,0)e^™<^. (74) 

If < (f) < 2tt, then we should tackle the condition of covering all space 
according to the requirement of dealing with a function continuous in any 
angular zone, i.e. ca 

Yim{9,(P = 0)=Yimie,(P = 27r), (75) 

implying 

e^'"" = 1. (76) 



^By definition, the hermitic operator A is an observable if the orthogonal system of 
eigenvectors form a base in the space of states. 

^Each quantum state of a particle is characterized by a vectorial state belonging to an 
abstract vectorial space Sr- 
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As has been seen, m is either an integer or a half- integer; for the apph- 
cation to the orbital momentum, m should be an integer, (e^*"*"^ would be 
— 1 if m is a half-integer) . 

(/3): For a given value of /, all the corresponding Yi^ can be obtained by 
algebraic means using 

l+Yu{e,^) = 0, (77) 
which combined with eq. (62) for leads to 

(^-^-lcote^Fn{e) = 0. (78) 

This equation can be immediately integrated if we notice the relationship 

d(sin 6) , , 

cot6'd6'=-^ — . 79 
smfc* 

Its general solution is 

F,; = Q(sin0)', (80) 

where q is a normalization constant. 

It follows that for any positive or zero value of Z, there is a function 
YlliQ., (j)), which up to a constant factor is 

Y^\9,,P) =ci{smeye'^'f'. (81) 

Using repeatedly the action of Z_ , one can build the whole set of functions 
Yii-i{9, 4>), . . . , Yio{6, (p), ... , Yi-i{6, (p). Next, we look at the way in which 
these functions can be put into correspondence with the eigenvalue pair 
l{l + l)h,m?i (where I is an arbitrary positive integer such that I < m < I 
). Using (78), we can make the conclusion that any other eigenfunction 
^lm{G-,4>) can be unambigously obtained from Yu. 

Properties of spherical hcirmonics 

a Iterative relationships 

Prom the general results of this chapter, we have 



l±Yim{e,ct>) = hJl{l + l) - m{m±l)Yim±i{d,4>)- (82) 
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Using (62) for l± and the fact that Yim{0, 4>) is the product of a 6'-dependent 
function and 6=*="^, one gets 



d_ 
89 



- m cot Yim{e, <i>) = ^Jl{l + l)-m{m±l)Yl^±l{9, cf) (83) 



(3 Orthonormalization and completeness relationships 

Equation (68) determines the spherical harmonics only up to a constant 
factor. We shall now choose this factor such that to have the orthonormal- 
ization of Yijn{d, <p) (as functions of the angular variables 9, (p) 



rZTV riv 

/ #/ sm9d9Yi*^{9,mm{0,cP) (84) 
Jo Jo 

In addition, any continuous function of 9, (j) can be expressed by means of 
the spherical harmonics as follows 

oo I 

f{9,<l>)=J2 12 ^lrnYlm{0,<l>), (85) 
;=o m=-l 

where ^ 

cirn = 1^ d(t> (\m9d9YilX0,<t>)f{0A)- (86) 
Jo Jo 

The results (85), (86) are consequences of defining the spherical harmon- 
ics as an orthonormalized and complete base in the space en of functions of 
9, (f). The completeness relationship is 

oo I 

J2J2^imi0^4>)Yi*m{e',(l>) = S{cos9-cos9')S{4>,<t^), 

1=0 m=l 

= ^s{o-m<i>,<i>)- (87) 

sm (7 

The 'function' d{cos 9 — cos 9') occurs because the integral over the variable 
9 is performed by using the differential element sin^d^ = —d{cos9). 

Pcirity operator V for spherical hcirmonics 

The behavior of V in 3D is rather close to that in ID. When it is applied to 
a function of cartesian coordinates tp{x, y, z) changes the sign of each of the 
coordinates 

Vi/jix, y, z) = ip{-x, -y, -z). (88) 
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V has the properties of a hermitic operator, being also a unitary operator, 
as well as a projector since "P^ is an identity operator 



(V(r)lPlV'(r)) = (V(r)lV(-r)) = (V(-r)|^(r')), 

V'^\r)=r{r\r))=r\-v) = \r). (89) 

Therefore 

= 1, (90) 

for which the eigenvalues are P = ±1. The eigenfunctions are called even 
if P = 1 and odd if P = — 1. In nonrelativistic quantum mechanics, the 
operator H for a conservative system is invariant with regard to discrete 
unitary transformations, i.e. 

VHV = V'^HV = H. (91) 

Thus, H commutes with V and the parity of the state is a constant of the 
motion. In addition, V commutes with the operators I and l± 

[P,y=0, [P,/±]=0. (92) 

Because of all these properties, one can have the important class of wave 
functions which are common eigenfunctions of the triplet P, P and f^. It 
follows from (92) that the parities of the states which difer only in 1^ coincide. 
In this way, one can identify the parity of a particle of definite orbital angular 
momentum I. 

In spherical coordinates, we shall consider the following change of vari- 
ables 

r^r, e ^-K -6 (p^7: + (/). (93) 

Thus, using a standard base in the space of wavcfunctions of a particle 
without 'intrinsic rotation', the radial part of the base functions ipklmi^ 
is not changed by the parity operator. Only the spherical harmonics will 
change. Prom the trigonometric standpoint, the transformations (93) are as 
follows 

sin(7r-0) ^ sine, cos(7r-e) ^ - cos e'"*('^+'^ ^ (_i)mgim<^ (94) 
leading to the following transformation of the function Yii{9, 0) 

Yii{ct> - e,Tr + cP) = {-1)%{9,^) . (95) 
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From (95) it follows that the parity of Yu is (— 1)^ On the other hand, l- 
(as well as Z+ is invariant to the transformations 



d{7r-e) de' a(7r + 0) 

(96) 

In other words, l± are even. Therefore, we infer that the parity of any 
spherical harmonics is (— 1)^ that is it is invariant under azimuthal changes 

YU(l)-9,7r + ct>) = {-iyYim{e,(t>). (97) 

In conclusion, the spherical harmonics are functions of well-defined parity, 
which is independent of m, even if I is even and odd if / is odd. 



The spin operator 

Some particles have not only orbital angular momentum with regard to 
external axes but also a proper momentum , which is known as spin denoted 
here by S. This operator is not related to normal rotation with respect to 
'real' axes in space, although it fulfills commutation relations of the same 
type as those of the orbital angular momentum, i.e. 

[Si, Sj] = ieijkSk , (98) 

together with the following properties 

(1) . For the spin operator all the formulas of the orbital angular momentum 

from (23) till (48) are satisfied. 

(2) . The spectrum of the spin projections is a sequence of either integer or 

half-integer numbers differing by unity. 

(3) . The eigenvalues of S^ are the following 

SV. = 5(S+1)V'.. (99) 



(4) . For a given S, the components Sz can take only 25" + 1 values, from 

-S to +S. 

(5) . Besides the usual dependence on r and/or p, the eigenfunctions of the 

particles with spin depend also on a discrete variable, (characteristic 
for the spin) a denoting the projection of the spin on the z axis. 
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(6). The wavefunctions ilj{f,a) of a particle with spin can be expanded in 
eigenfunctions of given spin projection Sz, i.e. 



s 



"0(^,0") = Yl V'a(r)x(<^) 



(100) 



a=-S 



where ^'^■(r) the orbital part and x{(^) is the spinorial part. 

(7) . The spin functions (the spinors) xi^i) ^-re orhtogonal for any pair ai ^ 

(Tfe. The functions ijjcr{r)x{<^) the sum of (100) are the components 
of a wavefunction of a particle with spin. 

(8) . The function ipcr{r) is called the orbital part of the spinor, or shortly 

orbital. 

(9) The normalization of the spinors is done as follows 



The commutation relations allow to determine the explicit form of the 
spin operators (spin matrices) acting in the space of the eigenfunctions of 
definite spin projections. 

Many 'elementary' particles, such as the electron, the neutron, the pro- 
ton, etc. have a spin of 1/2 (in units of h) and therefore the projection 
of their spin takes only two values, {Sz = ±1/2 (in h units), respectively. 
They belong to the fermion class because of their statistics when they form 
many-body systems. 

On the other hand, the matrices Sx, Sy, Sz in the space of S^, Sz are 



s 




(101) 



a=-S 







(102) 



Definition of the Pauli matrices 



The matrices 



= ^Si 



(103) 
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are called the Pauli matrices. They are hermitic and have the same charac- 
teristic eq. 

- 1 = 0. (104) 
Therefore, the eigenvalues of ax, CFy and cr^ are 

A = ±1. (105) 

The algebra of these matrices is the following 

(^i = I, (^k(^j = -<^j(^k = if^z, (^j(^k = i £jkm- + Sjkl ■ (106) 

I 

In the case in which the spin system has spherical symmetry 

Mr,+k), Mr,-k), (107) 

are different solutions because of the different projections Sz- The value 
of the probability of one or another projection is determined by the square 
moduli IIV'ilP or ||V'2|P, respectively, such that 

||V'i||' + ||V'2||' = l. (108) 
Since the eigenfunctions of Sz have two components, then 




(109) 



so that the eigenfunction of a spin one-half particle can be written as a 
column matrix 




(110) 



In the following, the orbitals will be replaced by numbers because we are 
interested only in the spin part. 



Transformations of spinors to rotations 

Let ijj be the wavefunction of a spin system in S. We want to determine 

the probability of the spin projection in a arbitrary direction in 3D space, 
which one can always chose as the z' of S'. As we have already seen in the 
case of the angular momentum there are two viewpoints in trying to solve 
this problem 
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(a) V does not change when S ^ S' and the operator A transforms as 
a vector. We have to find the eigenfunctions of the projections S'^ 
and to expand in these eigenfunctions. The square moduh of the 
coefficients give the result 

5; = 5^cos(/) + SySin0 = e-'^'f'S^e'^'t', 
S'y = - sin (t) + Sy cos (I) = e'^^'f' Syc'^'l' , 

= -S, = e'^'t'S,, (111) 

for infinitesimal rotations. Then, from the commutation relations for 
spin one can find 

i=S„ (112) 
where I is the infinitesimal generator. 

(/3) The second representation is: 

S does not change when T, ^ H' and the components of ip does change. 
The transformation to this representation can be performed through 
a unitary transformation of the form 

V^S'V = A, 
Using (111) and (113) one gets 

Y\ = ^iSA^ (114) 

and from (114) we are led to 

Using the explicit form of Sz and the properties of the Pauli matrices 
one can find the explicit form of Vj, such that 

m) = ( J, ) . (116) 
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A result of Euler 



One can reach any reference frame S' of arbitrary orientation with regard 
to E through only three rotations; the first of angle (j) around z, the next of 
angle 6 around x' and the last of angle jpa around z' , i.e. This important 
result belongs to Euler. The parameters {(p,6,i/ja) are called Euler's angles. 
Thus 

v\^,e,iPa) = vU^a)vii9)vH^)- (117) 

The matrices are of the form (116), whereas is of the form 

y^i^)=(-Tje TIV (118) 

so that 



t sin I cos ^ 



e 2 cos o «e 2 sm ■ 



yt(^,,,^„)= Z^'^'l . (119) 



, ., -7:1 

le 2 sm ^62 cos 2 



It comes out in this way that by the rotation of S, the components of 
the spinorial function transforms as follows 

ipi = tpie 2 cQs-+iip2e 2 sm-, 

. ■ ip — tpa ■ '^~\~'4'a 

V'2 = ^V'le'^^ sin - + ^2e~*^~ cos -. (120) 

Prom (120) one can infer that there is a one-to-one mapping between any 
rotation in £'3 and a linear transformation of E2, the two-dimensional Eu- 
clidean space. This mapping is related to the two components of the spinorial 
wavefunction. The rotation in does not imply a rotation in E2, which 
means that 

($V'> = (^IV'> = + ^1^2. (121) 

Prom (119) one finds that (121) does not hold; nevertheless there is an 
invariance in the transformations (119) in the space E2 of spinorial wave- 
functions 

= V'i$2-V'2$i- (122) 

The linear transformations that preserve invariant bilinear forms invari- 
ant arc called binary transformations. 

A physical quantity with two components for which a rotation of the 
coordinate system is a binary transformation is know as a spin of first order 
or shortly spin. 
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The spinors of a system of two fermions 

The eigenfunctions of jS^, with i = 1,2 have the following form 

A very used operator in a two-fermion system is the total spin 

S=iS+2S . (124) 
The spinors of Sz are kets \S,a), which are linear combinations of 




(125) 



The spinorial functions in (125) are assumed orthonormalized. In En 
the ket | + +) has Sz = 1 and at the same time it is an eigenfunction of the 
operator 

S =is^ + 2i,s){2s)+2s\ (126) 

as one can see from 

= \ + +) = ^\ + +) + 2{iS^ -2 3^+1 Sy -2 Sy+iSz -2 Sz)\ + +)(127) 

S'^ = I + +) = 2| + +) = 1(1 + 1)1 + +). (128) 
If we introduce the operator 

S_ =1 s_ +2 s_, (129) 

one gets 

[5_,52]=0. (130) 

Then (S'_)'^|l, 1) can be written in terms of the eigenfunctions of the operator 
i.e. 

5_|1, 1) = S-\ + +) = V2I + -) + V2I - +). (131) 
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Thus, Sz = in the state >S'_|1, 1). On the other hand, from the normahza- 
tion condition, we have 



M = j=,{\ + -) + \-+)) (132) 
5_|1,0) = I --) + !--) = a|l,-l). (133) 

In addition, the normahzation condition gives 

|1,-1) = |-,-). (134) 

There is only one other Uncar-independent combination of functions of 
the type (125), which is different of |1, 1), |1,0) and |1, —1), which is 

^4 = ^(1 + -)-!-+)), (135) 

Szip4 = 0, S'^ipi. (136) 

Therefore 

V'4 = |0,0). (137) 

■04 describes the state of a system of two fermions having the total spin equal 
to zero. The latter type of state is called singlet. On the other hand, the 
state of two fermions of total spin one can be called triplet having a degree 
of degeneration g = 3. 



Total angular momentum 

The total angular momentum is an operator defined as the sum of the an- 
gular and spin momenta, i.e. 

J = i+S, (138) 

where I and S, as we have seen, act in different spaces, though the square 
of I and S commute with J 

[JiJj] = ^e^jkJk, [JuP] = 0, [Ji,S^] = 0, (139) 

Prom (139) one finds that and S'^ have a common eigenfunction system 
with and J^. 
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Let us determine the spectrum of the projections for a fermion. The 
state of maximum projection can be written 

= ^ii(^l^=\l,l,+) (140) 

J.V' = il+'^)^P,^j = l + i (141) 
We introduce the operator J_ defined as 

j_=L+5_=L+(^ J (142) 

On account of the normalization a = -^(J + M){J — M + 1), one gets 



J_# J j =V2l\l,l-l,+) + \l,l-l,-), 



(143) 



so that the value of the projection of j_ in j^ip will be 

j, = (Z-l) + i = /-i . (144) 

It follows that j- lowers by one unit the action of J^. 
In the general case we have 

t = t + kt-^S- . (145) 

One can sec that (145) is obtained from the binomial expansion considering 
that si and all higher-order powers of s arc zero. 

=t\l,l,+) +ki1-%l,-). (146) 

Using 
we get 



I, +) = J^\i, i-k,+)+ jmrnikii, i-k+i,-). (147) 



(2i-fe)!l''''' ^/ ^ y (2i-fe+l)! 

Now noticing that m = 1 — k 



1-|M,+) = ^|^^^\hm,+) + ,|^^^{l-m)\l,m+l,-). (148) 
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The eigenvalues of the projections of the total angular momentum are given 
by the sequence of numbers differing by one unit from j = I + ^ pina to 
j = 1 — \. All these states belong to the same eigenfunction of J as |Z, Z, +) 
because [J_, J^] = 0: 

= {l^ + 2iS + S'^)\Ul,+), 

= [Z(Z + l) + 2Zi + |]|Z,/,+) (149) 

where + = + + 

In the left hand side of (149) a contribution different of zero gives only 
j = IzSz- Thus, the obtained cigenfunctions correspond to the pair j = 1 + ^, 
rrij = m + ^ ; they are of the form 

\ \ I ui \ I — jYl 

K + 2 , - + 2) = V ^[Tl-'^' + V + ^'''^ 

The total number of linearly independent states is 

iV = (2Z + l)(2s + l) = 4/ + 2, (151) 

of which in (150) only (2j+l)=21+3 have been built. The rest of 2Z - 1 
cigenfunctions can be obtained from the orthonormalization condition: 

K - i,m - i) = y^K,m, +) - yi±^|i,m + 1, -). (152) 

If two subsystems are in interaction in such a way that each of the 
angular momenta is conserved, then the cigenfunctions of the total angular 
momentum 

J = h+32, (153) 
can be obtained by a procedure similar to the previous one. For fixed eigen- 
values of ji and j2 there are (2ji + l)(2j2 + l) orthonormalized cigenfunctions 
of the projection of the total angular momentum J^; the one corresponding 
to the maximum value of the projection J^, i.e. Mj = ji + j2, can be built 
in a unique way and therefore J = ji + J2 is the maximum value of the 
total angular momentum of the system. Applying the operator J = ji + j2 
repeatingly to the function 

|ji +i2,ii + j2,ii +j2) = • |j2,i2), (154) 

one can obtain all the 2(ji + ^2) + 1 cigenfunctions of J = ji + j2 with 
different Ms: 

-(ii+i2)<M< (ii+ia). 
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For example, the eigenfunction of M = ji + j2 ~ 1 is 



(155) 

Applying iteratively the operator J_, all the 2{ji +j2 — l) — l eigenfunc- 
tions of J = jl + ^2 — 1 can be obtained. 
One can prove that 

\ji - 32\ < J <ji +32 , 

so that 

max ./ 

^ (2 J + 1) = (2Ji + 1)(2J2 + 1) . (156) 

min J 

Thus 

\J,M,ji,j2)= X! Ui'mij2m2\JM)\ji,mi,j2,m2) , (157) 

mi+m2=M 

where the coefficients {jiniij2m2\JM) determine the contribution of the 
various kets \ ji,mi, j2,m2) to the eigenfunctions of J^, having the eigen- 
values J{J + 1), M. They are called Clebsch-Gordan coefficients. 
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3N. Note: 1. The operator corresponding to the Runge-Lenz vector of the 
classical Kepler problem is written as 

(f xp)-{px i) , 

where atomic units have been used and the case Z = 1 (hydrogen atom) 
was assumed. This operator commutes with the Hamiltonian of the atomic 
hydrogen H = ^ — f > that is it is an integral of the atomic quantum motion. 
Its components have commutators of the type [^j, Aj] = —2ieijklk ■ H] the 
commutators of the Runge-Lenz components with the components of the 
angular momentum are of the type = ieij^Ak- Thus, they respect 

the conditions (23). Proving that can be a useful exercise. 



> f 1 

A = - + - 
r 2 
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3P. Problems 



Problem 3.1 

Show that any translation operator, for which il){y + a) = Taip{y), can be 
written as an exponential operator. Apply the result for y = f and for a the 
finite rotation a around z. 
Solution 

The proof can be obtained expanding tp{y + a)) in Taylor series in the in- 
finitesimal neighborhood around x, that is in powers of a 



'4){y + a) = Y^ 



oo „n jn 



n=Q 



We notice that 



E 

n=0 



and therefore one has Ta = e"'dx in the ID case. In 3D, y = r and a a. 
The result is Tg = e^'^ . 

For the finite rotation a around z we has y = (f) and a = a. It follows 

Ta = Ra = e"^ . 

Another exponential form of the rotation around z is that in terms of 
the angular momentum operator as was already commented in this chapter. 
Let x' = X + dx and consider only the first order of the Taylor series 



i^{x' ,y' ,z!) = ilj{x,y,z) + {x' - x)—'tp{x',y',z' 

d 

+{y' -y)^^{x',y',z') 
+{z'-z)-^^P{x',y',z') 



Taking into account 



x' = X — yd(f), 



y + xd(p, 
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one can reduce the series from three to two dimensions 

^(^) = ^{f) + {x-yd(l)-x)^^ + {y + xd(l)-yf'^^''^ 



dx 

ip[r) — ya<p — h xacpx 



dy' 



l-d(f) 



dx 
d 



d 



oy ox 



Since iL 



x-^ — y-S-) it follows that R 



dy 

1 - 



In the 



second order one can get ^{ilzdcp)^, and so forth. Thus, R can be written 
as an exponential 



R 



Problem 3.2 

Based on the expressions given in (14) show that one can get (15). 



A, 
A, 



Solution 

Let us consider only linear terms in the Taylor expansion (infinitesimal ro- 
tations) 

^iLd^ = 1 + i4# + ^{iUcpf + ... , 

so that 

(l + 4#)ix(l-4d<^) = Ax 

{Ax + ilzd(t)Ax)il - ilzd(j)) = Ax 

Ax — Axilzdcj) + ilzd(pAx + lzd(pAxlzd(j) = Ax 

i{lzAx - Axlz)d4) = -Aydcj). 

We easily arrive at the conclusion 

^zi Ax\ = iAy . 

In addition, [/2,^j/] = iAx can be obtained from 

(1 + ilzd(f))Ay{l — ilzd(f)) = Axd(f) — Ay, 

{Ay + ilzd(f)Ay){l — ilzdcf)) = A^dcf) — Ay, 

Ay — Ayilzdcj) + ilzd(l)Ay + lzd(l)Ay lzd(j) = Axdcf) — Ay , 

i{lzAy — Aylz)d(j) = —Axdcj). 
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Problem 3.3 

Determine the operator based on the Hamiltonian of an electron with 
spin in a magnetic field of induction B. 

Solution 

The Hamiltonian in this case is H{p,r,&) = H{p,r) + o" • B, where the 
latter term is the Zeeman Hamiltonian of the electron. Since cix commutes 
with the momenta and the coordinates, applying the Heisenberg equation 
of motion leads to 

Using [(7j;,cry] = ■icr^, one gets : 



—- = [PyBz - GzBy) = (CT X B)^ 

at nie rup. 
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4. THE WKB METHOD 

In order to study more realistic potentials with regard to rectangular 
barriers and wells, it is necessary to employ approximate methods allowing 
to solve the Schrodinger equation for more general classes of potentials and 
at the same time to give very good approximations of the exact solutions. 

The aim of the various approximative methods is to offer solutions of 
acceptable precision and simplicity that can be used for understanding the 
behaviour of the system in quasianalytic terms. 

Within quantum mechanics, one of the oldest and efficient approximate 
method for getting rather good Schrodinger solutions was developed almost 
simulataneously by G. Wentzel, H. A. Kramers and L. Brillouin in 1926, 
hence the acronym WKB under which this method is known (or JWKB as 
is more correctly used by many authors, see note 4N). 

It is worth mentioning that the WKB method applies to ID Schrodinger 
equations and that there are serious difficulties when trying to generalize it 
to more dimensions. 

In order to solve the Schrodinger equation 

with a potential of the form 

u{y) = u,f(y-) , (2) 
we first perform the changes of notations and of variable 

2 

""^ c' 

x = ?. (5) 



Prom eq. (5) we get 



dx dx dy dy 



(P d / d \ / d \/ d \ 2 rrj\ 
dx"^ dx^ dy' ^ dx'^ dx' dy"^ 
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and the Schrodinger eq. reads 



(8) 



Multiplying by —1/^^ and defining r{x) = rj — f{x), it is possible to write it 
as folows 



+ ^r(.)V^ = . 
To solve (9), the following form of the solution is proposed 



■0(a;) = exp 



- I q{x)da 



(9) 



(10) 



Therefore 



dx^ dx\x 
di^i/j i { i 2, ^ 

Factorizing V^, we have 

_ 

dx"^ 



rx 

I q{x)da 

J a 



- j q{x)da 



+ 



dq(x) 
dx 



exp 



- j q{x)da 

S Ja 



(11) 



Discarding for the time being the dependence of x, the Schrodinger eq. can 
be written 



1 



i dq 



1 



tp = 



and since in general ^ 0, we get: 



i^— + r - q = , 
ax 



(12) 



(13) 



which is a nonlinear differential eq. of the Riccati type whose solutions are 
sought in the form of expansions in powers of ^ under the assumption that 
^ is very small. 

More precisely, the series is taken of the form 



q{x) = Y,{-iCrQn{x) . 
n=Q 



(14) 
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Plugging it into the Riccati eq., we get 



/i=0 



fl5) 



n=0 /i=0 1^=0 

By a rearrangement of the terms one is led to 



n=0 



n+1 dqn 

dx 



;^i=0 i/=0 



(16) 



Double series have the following important property 



fi=0 u=0 



n=Q m=0 



where fi = n — m ,i' = m. 
Thus 



E(-i)"(^o 



n=0 



n+l dqn 

dx 



+ r{x) -EE (-^0""'"+"^9mgn-m = . (17) 
n=0 m=0 



Let us see explicitly the first several terms in each of the series in eq. 
(17): 

±i-irm-^'^ = + - + ... m 

• * riif rinc riry rim 

n=0 



dx dx dx dx 

oo n 

E E {-iO'^Qmqn-m = Qo - «2^«0gi 



(19) 



n=0 m=0 

Asking that the first terms in both series contain i^, one should write them 
as 



n=l 

which leads to 



n-1 



dx 



+ r{x) - (?o - E E i-iO^'Qmqn ■ 

n=l m=0 



m = , 



E 

n=l 



-HO"^-EHO' 



QmQn 



r{x) - ql 



. (20) 



m=0 
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In order that this equation be right the following conditions should be 
satisfied 

r(x)-g2 = qo = ±yf^ (21) 



m=0 

dqn 



dx 

- = - E IrnQn-m n > 1 . (22) 
m=0 

The latter is a recurrence relatioship, which occurs naturally in the WKB 
method. Recalling that we have defined r(a;) = ?7—/(x), "H = ^ ^ f{^) = 
by means of eq. (21) we get 



go 



This clearly indicates the classical nature of the WKB momentum of the 
particle of energy E in the potential u and units of ^/2rnuo. Thus 



go = Pix) = ^Jl^-f{x) 
is not an operator. If we approximate till the second order, we get 

q{x) =qo- iiqi - ^^q2 
and using the WKB recurrence relationship (22) we calculate qi and q2 

dx 2 go 2 dx 

«i = -^^(ln|M^)|) (24) 

7 dqi 2 

A glance to eq. (24), affords us to consider qi as the slope, up to a change 
of sign, of In I go I; when go is very small, then gi^O =^ — ^t/iS>0 and 
therefore the series diverges. To avoid this the following WKB condition 
is imposed 

\qo\ » I -^qi\ = ^\qi\ ■ 
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It is worth noting that this WKB condition WKB is not fulfilled at those 
points Xk where 

(li){xk) =P{xk) = . 



Since qo = P 



"^"Imuo^^ previous equation leads us to 
E = u{xk) . 



(26) 



In classical mechanics the points Xk that satisfies (26) are called turning 
points because the change of the sense of the motion of a macroscopic 
particle takes place there. 

By means of these arguments, we can say that qq is a classical solution of 
the problem under examination; also that the quantities qi k. q2 are the first 
and the second quantum corrections, respectively, in the WKB problem. 

To obtain the WKB wavefunctions we shall consider only the classical 
solution and the first quantum correction that we plug in the WKB form of 





i 




i/j = exp 


- / q{x)dx 

S Ja 


= exp 



t/j = exp I ^ y qodx j • exp 



rx 
J a 

qidx 



For the second factor, we get 

qidx , 



exp 



exp 



exp 



A 



Vp(x) 



where A is a constant, whereas for the first factor we get 



exp 



Ua 



qodx 



exp 



Ua 



p{x)dx 



Thus, we can write tp in the following form 

1 



4^ 



± 



Vp(x) 



exp 



/ p{x)da 

Ja 



(27) 
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The latter are known as the WKB solutions of the ID Schrddinger equation. 
The general WKB solution in the region in which the WKB condition is 
satisfied is written down as 

tjj = a+tjA + a^ijr . (28) 

As already mentioned there is no WKB solution at the turning points. This 
raises the question of the manner in which one has to do the passing from 
iIj{x < Xk) to ip{x > Xk). The solution of this difficulty is achieved by 
introducing the WKB connection formulas. 



The connection formulas 

We have already seen that the WKB solutions are singular at the clas- 
sical turning points; however, these solutions are correct both on the left 
and right side of these turning points x^. A natural question is how do we 
change ip{x < x^) in ip{x > x^) when passing through the turning points. 
The explicit answer is given by the connection formulas. 

Prom the theory of differential equations of complex variable it can be 
proved that really there are such connection formulas and that they are the 
following 

i^iix) = jexp ( - / J-r{x)dx) ->■ 

[-r(x)]4 \ Jx ' J 

^^L_cos( r J^)dx-^) , (29) 

where V'i(a^) tias only an attenuated exponential behavior for x < x^- The 
first connection formula shows that the function ipix), which at the left of 
the turning point behaves exponentially decaying, turns at the right of x^ 
into a cosinusoide of phase (p = j double amplitude with regard to the 
amplitude of the exponential. 

In the case of a more general function ipi^), such as a function with both 
rising and decaying exponential behavior, the connection formula is 

sin ((p + — Y- ( ( \/-r(x)dx] <— 

< cos ( I Jr{x)dx + (p\ , (30) 
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under the condition that (j) sa do not take a value that is too close to — J. 
The reason is that if ^ = — ^, then the sinus function is zero . The latter 
connection formula means that a function whose behavior is of the cosinusoid 
type at the right of a turning point changes into a growing exponential with 
sinusoid-modulated amplitude at the right of that point. 

In order to study the details of the procedure of getting the connection 
formulas we recommend the book Mathematical Methods of Physics by J. 
Mathews & R.L. Walker. 



Estimation of the WKB error 

We have found the solution of the Schrodingcr equation in the regions 
where the WKB condition is satisfied. However, the WKB solutions are 
divergent at the turning points. We thus briefly analyze the error introduced 
by using the WKB approximation and tackling the connection formulas in 
a close neighbourhood of the turning points. 

Considering x = as a turning point, we have qoix/.) = p^x/.) =0 ^ 
E = u{xk)- At the left of Xk, that is on the 'half-line' x < Xk, we shall 
assume E < u{x) leading to the WKB solution 



tp{x) = 




+ 



(31) 



Similarly, at the right of Xk (on the 'half-line x > Xk) we assume E > u{x); 
therefore the WKB solution in the latter region will be 




(32) 
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If 'ip{x) is a real function, it will have this property both at the right and 
the left of Xfe. It is usually called the "reality condition". It means that if 
a,b e^, then c = d*. 

Our problem consists in connecting the approximations on the two sides 
of Xk such that they refer to the same solution. This means to find c and d if 
one knows a and b, as well as viceversa. To achieve this connection, we have 
to use an approximate solution, which should be correct along a contour 
connecting the regions on the two sides of x/. , where the WKB solutions are 
also correct. A method proposed by Zwann and Kemhle is very useful in 
this case. It consists in going out from the real axis in the neighbourhood 
of Xk on a contour around x^ in the complex plane. It is assumed that on 
this contour the WKB solutions are still correct. Here, we shall use this 
method as a means of getting the estimation of the error produced by the 
WKB method. 

The estimation of the error is always an important matter for any ap- 
proximate solutions. In the case of the WKB procedure, it is more significant 
because it is an approximation on large intervals of the real axis that can 
lead to the accuulation of the errors as well as to some artefacts due to the 
phase shifts that can be introduced in this way. 

Let us define the associated WKB functions as follows 



W± 




(33) 



that we consider as functions of complex variable. We shall use cuts in 



order to avoid the discontinuities in the zeros of r{x) 



E-u{x) 



These 



functions satisfy the differential equation that is obtained by differentiating 
with respect to x, leading to 



Let us notice that 



Ir' 



4 r , 
/X 2 



r 1 r" 5 / r' 
^ aV ~ 16 



, , Ir" 5 fr' 



w± 

w± = o . 



(34) 



(35) 
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then W± are exact solutions of the equation 



1 



^r(x) + s{x) 



W± = 0, 



(36) 



although they satisfy only approximately the Schrodinger equation, which 
is a regular equation in x = x^, whereas the same equation for the associate 
WKB functions is singular at that point. 

We shall now define the functions a±{x) satisfying the following two 
relationships 

V'(x) = a+{x)W+{x) + a-{x)W-{x) (37) 



ijj'ix) = a+{x)W'^{x) + a-{x)WL{x) 



(38) 



where tp{x) is a solution of the Schrodinger equation. Solving the previous 
equations for a±, we get 



TpWL - Tp'W- 

W+WL - WLW- 



a- = 



'W+WL - WLW- ' 



where the numerator is just the Wronskian of tV+ and It is not difficult 
to prove that this takes the value —^i, so that a± simplifies to the following 
form 

a+ = |i {ipWL - V-'W-) 



Doing the derivative in x in the eqs. (39) and (40), we have 

In the brackets, the first and the fourth terms are zero; recalling that 



(39) 
(40) 



(41) 



< + ^r(x)V' = k W'l + 



1 



^r{x) + s{x) 



we can write eq. (41) in the form 
da± ^ 



dx 2 

da± 



-V' 



(42) 
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which based on eqs. (33) and (37) becomes 



da± ^ . s{x) 
dx 2 [r(x)] 



Oi± + a=p exp ^T^^ j \l r{x)dx^ 



(43) 



Eqs. (42) and (43) are useful for estimating the WKB error in the ID 
case. 

The reason for which can be considered as a measure of the WKB er- 
rors is that in the eqs. (31) and (32) the constants a, b and c, d, respectively, 
give only approximate solutions ijj^ while the functions a± when introduced 
in the eqs. (37) and (38) produce exact solutions. From the geometrical 
viewpoint the derivative gives the slope of the tangent to these functions 
and indicates the measure in which a± deviates from the constants a, 6, c 
and d. 

4N. Note : The original (J)WKB papers are the following: 

G. Wentzel, "Eine Verallgemeinerung der Wellenmechanik" , ["A generaliza- 
tion of wave mechanics"], 

Zeitschrift fiir Physik 38, 518-529 (1926) [received on 18 June 1926] 

L. Brillouin, "La mecanique ondulatoire de Schrodinger: une methode generale 
de resolution par approximations successives" , [ "Schrodinger 's wave mechan- 
ics: a general method of solving by succesive approximations"], 
Comptes Rendus Acad. Sci. Paris 183, 24-26 (1926) [received on 5 July 
1926] 

H. A. Kramers, "Wellenmechanik und halbzahlige Quantisierung" , ["Wave 
mechanics and half-integer quantization"], 

Zf. Physik 39, 828-840 (1926) [received on 9 Sept. 1926] 

H. Jeffreys, "On certain approx. solutions of linear diff. eqs. of the second 

order" , 

Proc. Lond. Math. Soc. 23, 428-436 (1925) 

4P. Problems 

Problem 4-1 

Employ the WKB method for a particle of energy E moving in a potential 
u{x) of the form shown in fig. 4.1. 
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Fig . 4.1 



Solution 

The Schrodinger equation is 

^ + ^[E-uix)]i. = 0. (44) 

As one can see, we have 

, > 2m , I is positive for a < x < b 

r(x) = [E — uix)] S ■ X- r , 
t)/' 1 IS negative lor x < a, a; > b. 

If V'(^) corresponds to the region x < o, when passing to the interval 
a < X <b, the connection formula is given by eq. (29) telling us that 



■i/;(x) J- COS 



where A is an arbitrary constant. 

When ^(x) corresponds to the region x > b, when passing to the segment 
a < X < 6, we have in a similar way 



\E-u\ 



4 
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where B is another arbitrary constant. The reason why the connection 
formula is again given by eq. (29) is easily understood examinining what 
happens when the particle reaches the second classical turning point at x = 
h. This produces the inversion of the direction of motion. Thus, the particle 
appears to come from the right toward the left. In other words, we are in 
the first case (from the left to the right), only that as seen in a mirror placed 
at the point x = a. 

These two expressions should be the same, independently of the con- 
stants A and B, so that 



cos 



cos I I J^^{E — u)dx — ^ I + cos 



Recalling that 



:A + B\ fA-B 
cos A + cos B = 2 cos ( — - — | cos 



= . 

(47) 



eq. (47) can be written 
2 cos 



•cos 




{E-u)dx-- + J — {E-u)dx-- 



n^^ ' 4 



^{E -u)dx + ^ 



0, (48) 



which implies that the arguments of the cosinusoids are integer multiples of 
^. On the other hand, the argument of the second cosinusoid do not lead to 
a nontrivial result. Therefore, we pay attention only to the argument of the 
first cosinusoid, which prove to be essential for getting an important result 



' 2m 



TT 



2m 



TT 



n 



—TT for n odd 
2 



/ 



1 2m 



TT 

2 {E — u)dx — — = riTT 

a \ h 2 
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- u)dx = (n + -)7r 
^ / j2m{E-u)dx = (n + hirh . (49) 

Ja ^ 

This result is very similar to the Bohr - Sommerfeld quantization rules. 

We recall that Bohr's postulate says that the orbital angular momen- 
tum of an electron moving on an "allowed atomic orbit" is quantized as 
L = nh, n = 1,2,3, .. .. We also recall that the Wilson - Sommerfeld quan- 
tization rules assert that any coordinate of a system that varies periodically 
in time should satisfy the 'quantum' condition: § Pqdq = riqh, where g is a 
periodic coordinate, pq is the associated momentum, riq is an integer, and h 
is Planck's constant. One can see that the WKB result is indeed very similar. 



Problem 4-2 

Estimate the error of the WKB solution WKB at a point xi ^ Xk, where 
Xk is a classical turning point for the differential equation y" + xy = 0. The 

solution of this problem is of importance in the study of uniform fields, such 
as the gravitational and electric fields generated by large planes. 



Solution: 

For this differential equation we have 

5 

^ = 1, r{x)=x k s{x) = -—x~'^. 

r{x) = x has a single zero at = 0, therefore for x ^ 0: 

W±=x~4exp^±i J ^/xdx^ = x~i exp ^±^ix^^ . (50) 

Derivating W± up to the second derivative in a;, we realize that the following 
differential equation is satisfied 

+ (x - ^x-^)W± = . (51) 

The exact solution y{x) of the latter differential equation can be written 
as a linear combination of W±, as it has been shown in the corresponding 
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section where the WKB error has been tackled; recall that the following 
form of the linear combination was proposed therein 

y{x) = a+(x)W+(a;) + a-{x)W-{x) 



For large x, the general solution of the differential equation can be writ- 
ten in the WKB approximation as follows 



y{x) 



Ax 4 cos 



2 3 , 
-X2 + d 



for X ^ oo 



(52) 



Thus, a+ —>■ ^e^^ and a- — > ^e~** for x —>■ oo. We want to calculate the 
error due to these WKB solutions. A simple measure of this error is the 
deviation of a+ and of a_ relative to the constants A. Using the equation 



da 



dx 



six) 
2*Vf(^ 



a± + azf exp (^2i J \Jr{x)dx 



and doling the corresponding substitutions, one gets 



dx 



^ -2 

16" 



_ i_ 

X 2 



2 ^ 2 



2 3 

exp ( ^2i-x2 



(53) 



Taking Aa-i- as the changes displayed by a± when x varies between xi and 
oo, we can do the required calculation by means of 



Aq;± 
^472 



2 r'^ 
A Jxi 



da± 
dx 



dx 



32 



±i5 



2 -S 
-X-, 



+ e 



J x\ 



4 3 

X 2 exp ( ^z-x2 



dx 



(54) 



The second term in the parentheses is less important than the first one 
because the complex exponential oscillates between 1 and —1 and therefore 
x~2 < a;~2. Consequently 



A/2 



, 5 . 

± — le 
48 



±ii5„ 2 
Xi 



(55) 



and as we can see the error introduced by the WKB approximation is in- 
deed small if we take into account that the complex exponential oscillates 

_3 

between —1 and 1, while X]^ ^ is also small. 
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5. THE HARMONIC OSCILLATOR (HO) 

The solution of the Schrodinger eq. for HO 

The HO can be considered as a paradigm of Physics. Its utihty is manifest 
in many areas from classical physics until quantum electrodynamics and 
theories of gravitational collapse. 

Prom classical mechanics we know that many complicated potentials can be 
well approximated near their equilibrium positions by HO potentials 

V{x)^\v"{a){x-af . (1) 

This is a ID case. For this case, the classical Hamiltonian function of a 
particle of mass m, oscillating at the frequency ui has the following form: 

and the quantum Hamiltonian corresponding to the space of configurations 
is given by 

2m^ dx' 2 ^ ' 

1 2 2 

Since we consider a time-independent potential, the eigenfunctions 
and the eigenvalues En are obtained by means of the time-independent 
Schrodinger equation 

H^n = En^n ■ (5) 

For the HO Hamiltonian, the Schrodinger eq. is 



dx'^ 



2mE vr?ijJ^ ^ 

~1? tT'' 



^ = {) . (6) 



We cancealed the subindices of E and because they are not of any 
importance here. Defining 

e^'^ (7) 
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the Schrodinger eq. becomes 



+ [r - Yx^]^ = , (9) 



which is known as Weber's differential equation in mathematics. 
We shall make now the transformation 

y = Ax^ . (10) 

In general, by changing the variable from x to y, the differential operators 
take the form 

dx dx dy 

= A(^iL) = + (!^)2_^f_ (-12) 

dx'^ dx dx dy dx^ dy dx dy"^ 

Applying this obvious rule to the proposed transformation we obtain the 
following differential eq. in the y variable 

d^^r id^ 1 ^ 



and, by definind : 
we get 



e p E 

K= TTT = t: = 1— , (14) 

2A 2muj ^ ' 



Let us try to solve this equation by first doing its asymptotic analysis in 
the limit y — ^ 00. We first rewrite the previous equation in the form 

(f-^ Id^ ,K l,^ ^ 
We notice that in the limit y ^ 00 the equation behaves as follows 



dy'^ 4 
This equation has as solution 



- 7*00 = . (17) 



*oo(y) =^exp| + Bexp^ . (18) 
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Taking ^ = 0, we eliminate exp | since it diverges in the limit y — > oo, 
keeping only the attenuated exponential. We can now suggest that ^ has 
the following form 

= exp ^V(2/) • (19) 
Plugging it in the differential equation for y ( eq. 15) one gets: 



The latter is a confluent hyper geometric equation Q : 

The general solution of this equation is 

y{z) = A iFi (a; c,z)+B z^"^ (a - c + 1; 2 - c, z) , (22) 
where the confluent hypergeometric equation is defined by 

.^.c;-)-^!^ . (23) 

n=0 ^ ' 

Comparing now our equation with the standard confluent hypergeomet- 
ric equation, one can see that the general solution of the first one is 

1 1 13 

iP{y) = A ^Fi{a;-,y) + B y^ ^Fi{a + -■ -,y) , (24) 

where 

« = -(2-4)- (25) 

If we keep these solutions in their present form, the normalization con- 
dition is not satisfied for the wavefunction because from the asymptotic 
behaviour of the confluent hypergeometric function ^ it follows ( taking into 
account ony the dominant exponential behavior ) : 

— y y 1 

^{y) = e~ilj{y) const. e2y"~2 . (26) 



^It is also known as Kummer's differential equation. 
^ The asymptotic behavior for | a:: |— > oo is 



iFi(a; c, z) jilsLe— x"" + f^e^x-^"^ 



91 



The latter approximation leads to a divergence in the normalization in- 
tegral, which physically is not acceptable. What one does in this case is to 
impose the termination condition for the series ^ , that is , the series has 
only a finite number of terms and therefore it is a polynomial of n order. 
We thus notice that asking for a finite normalization constant (as already 
known, a necessary condition for the physical interpretation in terms of prob- 
abilities), leads us to the truncation of the series, which simultaneously gen- 
erates the quantization of energy. 
In the following we consider the two possible cases 

1) a = —n and B = 

K 1 , , 

2-4=^- ^''^ 
The eigenfunctions are given by 

— Ax^ 1 
= D„exp^— iFi(-n;-,Ax2) (28) 

and the energy is: 

En = huj{2n + ^) . (29) 

2) a + i = -n and ^ = 



K 1 1 , , 

2-4=^+2- 



The eigenfunctions are now 



-Xx^ 3 
^n(x) = Dnexp ^ X iFi(-n; -, Ax^) , (31) 

whereas the stationary energies are 

En = nio[{2n + l) + ^] . (32) 

The polynomials obtained by this truncation of the confluent hypergeo- 
metric series are called Hermite polynomials and in hypergeometric notation 
they are 

H2n{r]) = (-1)"^ iFi(-n;i r/2) (33) 



®The truncation condition of the confluent hypergeometric series iFi(a; c, z) is a — —n, 
where n is a nonnegative integer ( i.e., zero included). 
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We can now combine the obtained results ( because some of them give 
us the even cases and the others the odd ones ) in a single expression for 
the eigenvalues and eigenfunctions 

= D„ exp -^HniVXx) (35) 

En = {n+^)hLO n = 0,l,2 ... (36) 

The HO energy spectrum is equidistant, i.e., there is the same energy 
difference Tilo ibetween any consequitive neighbour levels. Another remark 
refers to the minimum value of the energy of the oscillator; somewhat sur- 
prisingly it is not zero. This is considered by many people to be a pure 
quantum result because it is zero when ^ — 0. It is known as the zero point 
energy and the fact that it is different of zero is the main characteristic of 
all confining potentials. 



The normalization constant is easy to calculate 

^ 1 

l\ 1 

Dr,. = 



(37) 



tor 



7r2"n! 

Thus, one gets the following normalized eigenfunctions of the ID opera- 

A 1 



7r2'*n! 



-Xx'^ 

exp(^— ) HniVXx) . 



(38) 



Creation and anihilation operators: a) and a 

There is another approach to deal with the HO besides the conventional one 
of solving the Schrodinger equation. It is the algebraic method, also known 
as the method of creation and annihilation (ladder) operators. This is a very 
efficient procedure, which can be successfully applied to many quantum- 
mechanical problems, especially when dealing with discrete spectra. 
Let us define two nonhermitic operators a and : 

Imu! , ip , 
V 2ft mu 
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V 2n ^ 



These operators are known as annihilation operator and creation 
operator, respectively (the reason of this terminology will be seen in the 
following, though one can claim that it comes from quantum field theories) . 
Let us calculate the commutator of these operators 

, 4., mu! , ip ip , 1 / -r 1 -r i\ 1 

[a,a^\ = i^b^H ,x = ^{-i[x,p\ +i\p,x\) = 1 , (41) 

where we have used the commutator 

[x,p]=in. (42) 

Therefore the annihilation and creation operators do not commute, since 
we have 

[a,a)] = l. (43) 
Let us also introduce the very important number operator N: 

N = a^a . (44) 

This operator is hermitic as one can readily prove using {AB)^ = B^A^ 

= {a^a)^ = a\a^y = a^a = N . (45) 
Considering now that 

a^a = — + + — = — - - 46 

2n m'^u)'^ zh huj 2 

we notice that the Hamiltonian can be written in a quite simple form as a 
function of the number operator 

H = hu{N + ^) . (47) 

The number operator bear this name because its eigenvalues are precisely 
the subindices of the eigenfunctions on which it acts 

N \ n) = n\ n) , (48) 

where we have used the notation 

Psin = \n) . (49) 
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Applying this fact to (47), we get 

H I n >= ?iu{n + ^) I n > . (50) 

On the other hand, from the Schrodinger equation we know that H \ 
n >= E \ n >. In this way, it comes out that the energy eigenvalues are 
given by 

En = nu{n + ^) . (51) 

This result is identical (as it should be) to the result (36). 
We go ahead and show why the operators a and a) bear the names they 
have. For this, we calculate the commutators 

[N,a] = [a) a, a] = a''" [a, a] + [a''",a]a = —a , (52) 

which can be obtained from [a, a] = and (43). Similarly, let us calculate 

[N,a^] = [a^a,a^] = a^[a,a^] + [a\a^]a = J . (53) 

Using these two commutators, wc can write 

n > 

= {a) + a)N) I n > (54) 
= a\l + n) I n >= (ra + l)a^ \ n> . 

By a similar procedure, one can also obtain 

N[a \n>) = ([iV, a] + aN) \ n >= {n - l)a \ n > . (55) 

The expression (54) implies that one can consider the ket a''' | n > as an 
eigenket of that number operator for which the eigenvalue is raised by one 
unit. In physical terms, this means that an energy quanta has been produced 
by the action of on the ket. This already expains the name of creation 
operator. Similar comments with corresponding conclusion can be infered 
for the operator a, originating the name of annihilation operator (an energy 
quanta is eliminated from the system when this operator is put in action). 
Moreover, eq. (54) implies the proportionality of the kets \ n > and 
I n + 1 >: 

a} I n >= c I 77, + 1 > , (56) 
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where c is a constant that should be determined. Considering in addition 
(at I n >)t =< n\a = c*<n + l\ , (57) 
one can perform the fohowing calculation 

< n I a(at \n>) = c* <n + l\ (c | n + 1 >) (58) 

< n I aat | n >= c*c <n + l|?7, + l> (59) 
< n I aflt I n >=| c p . (60) 
But from the commutation relation for the operators a and 

[a, at] = aat — ata = aat — N = 1 , (61) 

we have 

aat = AT + 1 . (62) 
Substituting in (60), we get 

<n \ N + l\ n >=< n \ n > + < n \ N \ n >= n + I =\ c l"^ . (63) 

Asking conventionally for a positive and real c, the following value is 
obtained 

c = Vn + 1 . (64) 

Consequently, we have 

at I n >= Vn + 1 \n + l> . (65) 

For the annihilation operator, following the same procedure one can get 
the following relation 

a I n >= \fn \ n — \> . (66) 

Let us show now that the values of n should be nonnegative integers. 
For this, we employ the positivity requirement for the norm, applying it 
to the state vector a\n >. The latter condition tells us that the interior 
product of the vector with its adjunct ((a | n >)t =< n \ at) should always 
be nonnegative 

(< Ti I at) • (a I n >) > . (67) 
This relationship is nothing else but 

<n\a^a\n >=< n | TV | n >= n > . (68) 
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Thus, n cannot be negativ. It should be an integer since were it not by 
applying iteratively the annihilation operator we would be lead to negative 
values of n, which would be a contradiction to the previous statement. 
It is possible to express the state n (| n)) directly as a function of the 
ground state (| 0)) using the creation operator. Let us see how proceeds this 
important iteration 



2) 
3) 



r , , , 



2) 



1) =aT 
^ V2!^ 



0) 
0) 

0) 



(69) 
(70) 

(71) 



n) = [^] I 0) . (72) 
V n! 



One can also apply this method to get the eigenfunctions in the config- 
uration space. To achieve this, we start with the ground state 

a I 0) = . (73) 

In the X representation, we have 



»oW = ^^(x+£;)>l.„(x) = 0. (74) 

Recalling the form of the momentum operator in the x representation, we 
can obtain a differential equation for the wavefunction of the ground state. 
Moreover, introducing the definition xq = ^ we have 

{x + xl^)^o = 0- (75) 

The latter equation can be readily solved, and normalizing (its integral from 
— oo to oo should be equal to unity), we obtain the wavefunction of the 
ground state 

^o{x) = ( )e"^^^^' . (76) 
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The rest of the cigcnfunctions, which describe the HO excited states, can be 
obtained employing the creation operator. The procedure is the following 

*- = "'*° = <7k"""^°^'*° '''' 

By mathematical induction, one can show that 

1 1 r, d 1/ x ^2 

Time evolution of the oscillator 

In this section wc shall ilhistratc on the HO example the way of working with 
the Hciscnbcrg representation in which the states are fixed in time and only 
the operators evolve. Thus, we shall consider the operators as functions 
of time and obtain explicitly the time evolution of the HO position and 
momentum operators, a and a\ respectively. The Heisenberg equations of 
the motion for p and x are 

Hence the equations of the motion for x and p in the HO case are the 
following 

^ = -muj^x (82) 

f = I. (83) 
dt m ^ ^ 

These are a pair of coupled equations, which are equivalent to a pair of 

uncoupled equations for the creation and annihilation operators. Explicitly, 

we have 



da Imio d , ^ ip . ,„ 

di = l-mJt^^^^^ ^''^ 

da jmuj dx i dp 

dt ^ y ^^~di ' ^ ^ 
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Substituting (82) and (83) in (85), we get 



da jmuj p . . 

— = \ — :- iLOX) = —iLoa . (oo) 

Similarly, one can obtain a differential equation for the creation operator 

da} + , 

— = lua^ . (87) 

The differential evolution equations for the creation and annihilation opera- 
tors can be immediately integrated leading to the explicit evolution of these 
operators as follows 

a{t) = a(0)e-^'^* (88) 
a\t) = a\0)e^^ . (89) 

It is worth noting based on these results and eqs. (44) and (47) that 
both the Hamiltonian and the number operator are not time dependent. 

Using the latter two results, we can obtain the position and momentum 
operators as functions of time as far as they are expressed in terms of the 
creation and annihilation operators 



X = J^{a + a^) (90) 
V zniuj 



p ^ J"J^(a'-a). (91) 



Substituting them, one gets 



x(t) = x(0)cosa;t + ffisina;t (92) 

mu> 



p{t) = —mujx{0)smu}t+p{0)cosu}t. (93) 

The time evolution of these operators is the same as for the classical 
equations of the motion. 

Thus, we have shown here the explicit evolution form of the four HO ba- 
sic operators, and also we illustrated the effective way of working in the 
Heisenberg representation. 
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The 3D HO 



We commented on the importance in physics of the HO at the very beginning 
of our analysis of the quantum HO. If we will consider a 3D analog, we would 
be led to study a Taylor expansion in three variables[[] retaining the terms 
up to the second order, we get a quadratic form in the most general case. 
The problem at hand in this approximation is not as simple as it might look 
from the examination of the corresponding potential 

V{x, y, z) = ax^ + by^ + cz'^ + dxy + exz + fyz . (94) 

There are however many systems with spherical symmetry or for which 
this symmetry is sufficiently exact. In acest caz: 

V{x,y,z)=K{x^ + y^ + z^) , (95) 

which is equivalent to saying that the second unmixed partial derivatives 

have all the same value, denoted by K in our case). We can add that this is 

a good approximation in the case in which the values of the mixed second 

partial derivatves are small in comparison to the unmixed ones. 

When these conditions are satisfied and the potential is given by (95), we 

say that the system is a 3D spherically symmetric HO. 

The Hamiltonian in this case is of the form 

where the Laplace operator is given in spherical coordinates and r is the 
spherical radial coordinate. 

Since the potential is time independent the energy is conserved. In addition, 
because of the spherical symmetry the orbital momentum is also conserved, 
having two conserved quantities, we may say that to each of it one can 
associate a quantum number. Thus, we can assume that the eigenfunctions 
depend on two quantum numbers (even though for this case we shall see 
that another one will occur). Taking care of these comments, the equation 
of interest is 



^It is possible to express the Taylor series in the neighbourhood of ro as an exponential 
operator 
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H^nl = Enl^nl ■ (97) 

The Laplace operator in spherical coordinates reads 

and can be also inferred from the known fact 

= -^'[-^^(^i^^^) + -^2^2] ■ (99) 
sm 06 06 sm 6^ oip^ 

The eigenf unctions of are the spherical harmonics, i.e. 

L'^Yi^, {9, <p) = + l)y^n,, {9, <f) (100) 

The fact that the spherical harmonics 'wear' the quantum number mi 
introduces it in the total wavefunction ^nimi ■ 

In order to achieve the separation of the variables and functions, the follow- 
ing substitution is proposed 

*n;m, {r, 9, if) = ^^Yim, {9, if) . (101) 

Once this is plugged in the Schrodinger equation, the spatial part is 
separated from the angular one; the latter is identified with an operator 
that is proportional to the square of the orbital momentum, for which the 
eigenfunctions are the spherical harmonics, whereas for the spatial part the 
following equation is obtained 



<, + (i^-^.^-^)«»<M = o. (102) 



,2mEni rn^uJ^ 2 + 

Using the definitions (7) and (8), the previous equation is precisely of the 
form (9), unless the angular momentum term, which is commonly known as 
the unghiular, care in mod comun se cunoa§te ca angular momentum harrier 

+ {k' - - ^-^^)Rni = . (103) 

To solve this equation, we shall start with its asymptotic analysis. If we 
shall consider first r ^ 00, we notice that the orbital momentum term is 
negligible, so that in this limit the asymptotic behavior is similar to that of 
(9), leading to 
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-Ar2 

Rni{r) ~ exp — ^ — for lim r — > oo . (104) 

If now we pass to the behavior close to zero, we can see that the dominant 
term is that of the orbital momentum, i.e., the differential equation (102) in 
this limit turns into 

- ^-^^Rnl = . (105) 

This is a differential equation of the Euler type|^ , whose two independent 
solutions are 

Rnl{r) ~ r'+^ or for lim r ^ . (106) 

The previous arguments lead to proposing the substitution 

Rni{r)=r'+^e^V^^<P{r) ■ (107) 

One can also use another substitution 

_ w2 

Rni (r ) = r-' exp -^v{r) , (108) 

which, however, produces the same solutions as (107) (showing this is a help- 
ful exercise). Substituing (107) in (103), the following differential equation 
for (j) is obtained 

f + 2(^-^ - \r)4>' - [\{2l + 3) - k^]4) = . (109) 

Using now the change of variable w = Ar^, one gets 

w<t^" + + w)ct>' - + ^) - |]</> = , (110) 



*An equation of the Euler type has the form 

a;"?/'"' (a;) + x^~^y^"~^\x) + ■ • • + xy'{x) + y(x) = . 

Its solutions are of the type x" that are plugged in the equation obtaining a polynomial 
in a. 
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where = |x ~ ^ been introduced. We see that we found again a dif- 
ferential equation of the confluent hypergeometric type having the solutions 
(see (21) and (22)) 

^{r) = A ,F,[\{l+^-ny,l+l,Xr']+B r-(^'+i) ,F,[\{~l+^-Ky, Xr^] . 

(Ill) 

The second particular solution cannot be normalized because diverges 
strongly in zero. This forces one to take B = 0, therefore 

<P{r) = A ,F,[^(l + ^-ny,l + l,Xr'] • (112) 

Using the same arguments as in the ID HO case, that is, imposing a regular 
solution at infinity, leads to the truncation of the series, which implies the 
quantization of the energy. The truncation is explicitly 

l(l + l-'^) = -n, (113) 
where introducing k we get the energy spectrum 

Eni = nco{2n + 1 + 1)- (114) 

One can notice that for the 3D spherically symmetric HO there is a zero 
point energy ^hoo. 

The unnormalized eigenfunctions are 

^nlm{r,e,ip)=r^e^ iFi{-n;l + ^, Xr^) Yi^{e,ip) . (115) 

5P. Problems 

Problem 5.1 

Determine the eigenvalues and eigenfunctions of the HO in the 
momentum space. 

The quantum HO Hamiltonian reads 

n = h -muj X . 

2m 2 

In the momentum space, the operators x and p have the following form 

P 
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dp 

Thus, the HO quantum Hamiltonian in the momentum representation is 

H = -mw h -r-^ . 

2m 2 dp^ 

We have to solve the eigenvalue problem (i.e., to get the eigenfunctions and 
the eigenvalues) given by (5), which, with the previous Hamiltonian, turns 
into the following differential equation 

+ ^2-^ - * W = . 116 

dp^ mWu^ mP'Wu^ 

One can see that this equation is identical, up to some constants, with the 
differential equation in the space of configurations (eq. (6) ). Just to show 
another way of solving it, we define two parameters, which are analogous to 
those in (7) and (8) 

k' = ^^ \=^. (117) 

With these definitions, we get the differential cq. (9) and therefore the 
solution sought for (after performing the asymptotic analysis) is of the form 

^r(y) = e-^^0(y) , (118) 

where y = Xp^ and A is defined in (117). Substitute (118) in (116) taking 
care to put (118) in the variable p. One gets a differential equation in the 
variable ^ 

^-2Ap^?*>+(*^-A)«p) = 0. (119) 

We shall now make the change of variable u = \fXp that finally leads us to 
the Hermite equation 

^_2.^ + 2n^(n)=0, (120) 
du^ du 

where n is a nonnegative integer and where we have put 

^ - 1 = 2n . 
A 
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From here and the definitions given in (117) one can easily conclude that 
the eigenvalues are given by 

En = huj{n + ^) . 

The solutions for (120) are the Hermite polynomials (j){u) = Hn{u) and the 
unnormalized eigenfunctions are 



Problem 5.2 

Prove that the Hermite polynomials can be expressed in the fol- 
lowing integral representation 

Hn{x) = ^ / (x + iyre-y dy . (121) 

This representation of Hermite polynomials is not really usual, though 
it can prove useful in many cases. In order to accomplish the proof, we shall 
expand expand the integral and next prove that what we've got is identical 
to the series expansion of the Hermite polynomials that reads 

[HI 

where the symbol [c], indicating where the series terminates, denotes the 
greatest integer less or equal to c. 

The first thing we shall do is to expand the binomial in the integral by using 
the well-known binomial theorem 

(n — my.mi 

m=0 ^ ' 

Thus 

{x + tyT = E 7 , (123) 

[n — m)\m\ 

m=0 ^ ' 
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which plugged in the integral leads to 

Inspecting of the integrand we realize that the integral is not zero when 
m is even, whereas it is zero when m is odd. Using the even notation m = 2k, 

we get 

on '•2 J I ,-oo 

Under the change of variable u = y^, the integral turns into a gamma 
function 

f-1 

9" ^ n\ r°° 1 

±= y !h^^i2k^n-2k / ^fe-ig-^rf^ (126) 

V^fr'o{n-2ky.{2k)\ Jo ^ \ J 

more precisely r{k + ^), which can be expressed in terms of factorials ( of 
course for k an integer) 

^„ 1, (2fc)! ^ 
Plugging this expression in the sum and using i'^^ = (—1)'^, one gets 

which is identical to (122), hence completing the proof. 
Problem 5.3 

Show that Heisenberg's uncertainty relation is satisfied by doing 
the calculation using the HO eigenfunctions . 

We have to show that for any we have 

<{Ap)\Axf> >^, (128) 
where the notation <> means the mean value. 

We shall separately calculate < (Ap)^ > and < (Ax)^ >, where each of 
these expressions is 
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< {Apf >=< {p- < p >f >=< p^-2p <p> + <p >^>=< p^ > - <p>'^ , 



< (Ax)^ >=< [x- < X >)^ >=< x'^-2x < X > + < X >^>=< x'^ > - < X 

First of all, we shall prove that both the mean of x as well as of p are 
zero. For the mean of x, we have 



< X >= / x[^nix)] dx . 

J —oo 

This integral is zero because the integrand is odd. Thus 

< X >= . (129) 

The same argument holds for the mean of p, if we do the calculation in 
the momentum space, employing the functions obtained in problem 1. It is 
sufficient to notice that the functional form is the same (only the symbol 
does change). Thus 

<p>=0 . (130) 

Let us now calculate the mean of x^. We shall use the virial theorem 0. 
We first notice that 

< V >= -muj"^ < x^ > . 
2 

Therefore, it is possible to relate the mean of x^ directly to the mean of the 
potential for this case (implying the usage of the virial theorem). 



2 

We also need the total energy 



< x^ >= —2 <V > . (131) 



< H >=< T> + <V > 



®We recall that the virial theorem in quantum mechanics asserts that 

2 < T >=< r ■ v''^(r) > . 
For a potential of the form V = \x", the virial theorem gives 

2<T>=n<V > , 
where T is the kinetic energy and V is the potential energy. 
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for which again one can make use of the virial theorem (for n = 2) 

<H >=2<V > . (132) 

Thus, we obtain 

<.->=^ = ^^^^ (133) 

<x2>=— (n + -). (134) 
Similarly, the mean of can be readily calculated 

<p'- >=2m< ^ >=2m <T >=m <B. >= mhujin + -) . (135) 
2m 2 

Employing (133) and (135), we have 

< (Ap)2(Aa;)2 >= (n + )^^n^ . (136) 

Based on this result, we come to the conclusion that in the HO stationary 
states that actually have not been directly used, Heisenberg's uncertainty 
relation is satisfied and it is at the minimum for the ground state, n = 0. 

Problem 5.4 

Obtain the matrix elements of the operators a, a^, x, and p. 

Let us first find the matrix elements for the creation and annihilation 

operators, which arc very helpful for all the other operators. 
We shall use the relatinships (65) and (66), leading to 

< m I a I n >= \fn <m\n— \ >= \/n5m,n-i ■ (1^'^) 

Similarly for the creation operator we have the result 

< m \ \n >= \/n + l <m\n + l >= Vn + 16m,n+i ■ (1>^8) 

Let us proceed now with the calculation of the matrix elements of the 
position operator. For this, let us express this operator in terms of creation 
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and annihilation operators. Using the definitions (39) and (40), one can 
immediately prove that the position operator is given by 



^ (a + a^) . (139) 



2muj 



Employing this result, the matrix elements of the operator x can be 
readily calculated 



< 'm \ X \ n > = < m \ \ (a + a^) I n > 

V 2mu 



2mLjj 



Wn5m,n-i + \/n + l5m,n+i] ■ (140) 



Following the same procedure we can calculate the matrix elements of the 
momentum operator, just by taking into account that p is given in terms of 
the creation and annihilation operators as follows 



P = nr^(a*-«). (141) 



This leads us to 



<m\p\n>=i^ ^^^^^[Vn + lSm,n+i - VnSm,n-i] ■ (142) 

One can realize the ease of the calculations when the matrix elements of 
the creation and annihilation operators are used. Finally, we remark on the 
nondiagonality of the obtained matrix elements. This is not so much of a 
surprise because the employed representation is that of the number operator 
and none of the four operators do not commute with it. 



Problem 5.5 

Find the mean values of and for thelD HO and use them to 
calculate the mean (expectation) values of the kinetic and poten- 
tial energies. Compare the result with the virial theorem. 

First of all, let us obtain the mean value of x^. For this, we use eq. (139) 
that leads us to 

x'^ = -^(a'^ + {a'f + a'a + aa'). (143) 
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Recall that the creation and annihilation operators do not commute. Based 
on (143), we can calculate the mean value of 



< > = < n \ x'^ \ n > 



,n+2 



2mu) 

+ n Sn,n + (n + l) 6n,n] , (144) 
which shows that 

< >=< n\ x'^ \n>= --^—(2n + 1) . (145) 

2muj 

In order to calculate the mean value of we use (141) that helps us to 
express this operator in terms of the creation and annihilation operators 

f = -^(a^ + (at)2 - aat - a^a) . (146) 

This leads us to 

<p^ >=<n\p^ \n>= —^{2n + l) . (147) 
The latter result practically gives us the mean kinetic energy 

On the other hand, the mean value of the potential energy 

< V >=< ]^rrvj?x^ >= ^rnu? < x^ >= ^(2n + 1) , (149) 

where (145) has been used. 

We can see that these mean values are equal for any n, which confirms 
the quantum virial theorem, telling us that for a quadratic (HO) potential, 
the mean values of the kinetic and potential energies should be equal and 
therefore be half of the mean value of the total energy. 
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6. THE HYDROGEN ATOM 



Introduction 

In this chapter we shall study the hydrogen atom by solving the time- 
independent Schrodinger equation for the potential due to two charged par- 
ticles, the electron and the proton, and the Laplacian operator in spherical 
coordinates. Prom the mathematical viewpoint, the method of separation 
of variables will be employed, and a physical interpretation of the wavefunc- 
tion as solution of the Schrodinger equation in this important case will be 
provided, together with the interpretation of the quantum numbers and of 
the probability densities. 

The very small spatial scale of the hydrogen atom is a clue that the related 
physical phenomena enter the domain of applicability of the quantum me- 
chanics, for which the atomic processes have been a successful area since 
the early days of the quantum approaches. Quantum mechanics, as any 
other theoretical framework, gives relationships between observable quan- 
tities. Since the uncertainty principle leads to a substantial change in the 
understanding of observables at the conceptual level, it is important to have 
a clear idea on the notion of atomic observable. As a matter of fact, the 
real quantities on which quantum mechanics offers explicit answers and con- 
nections arc always probabilites. Instead of saying, for example, that the 
radius of the electron orbit in the fundamental state of the hydrogen atom is 
always 5.3 x 10~^^ m, quantum mechanics asserts that this is a truly mean 
radius (not in the measurable sense). Thus, if one performs an appropriate 
experiment, one gets, precisely as in the case of the common arrangement 
of macroscopic detectors probing macroscopic properties of the matter, ran- 
dom values around the mean value 5.3 x 10~^^ m. In other words, from the 
viewpoint of the experimental errors there is no essential difference with re- 
gard to the classical physics. The fundamental difference is in the procedure 
of calculating the mean values within the theoretical framework. 

As is known, for performing quantum-mechanical calculations, one needs 
a corresponding wave function Although * has no direct physical inter- 
pretation, the square modulus | ^' ^ calculated at an arbitrary position and 
given moment is proportional to the probability to find the particle in the 
infinitesimal neighbourhood of that point at the given time. The purpose of 
quantum mechanics is to determine for a specified particle in the prepared 
experimental conditions. 
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Before proceeding with the rigorous approaches of getting ^ for the 
hydrogen electron, we will argue on several general requirements regarding 
the wave function. First, the integral of | ^' p over all space should be finite 
if we really want to deal with a localizable electron. In addition, if 

/oo 
|*pcZy = 0, (1) 
-oo 

then the particle does not exist. | * P cannot be negative or complex 

because of simple mathematical reasons. In general, it is convenient to 
identify | ^' P with the probability P not just the proportionality. In order 
that I ^ p be equal to P one imposes 

/oo 
I * p czy = 1 , (2) 
-oo 

because 

/oo 
FdV = 1 (3) 
-oo 

is the mathematical way of saying that the particle exists at a point in 
space at any given moment. A wave function respecting eq. 2 is said to be 
normalized. Besides this, ^' should be single valued, because P has a unique 
value at a given point and given time. Another condition is that ^ and its 
partial first derivatives ^) ^ should be continuous at any arbitrary 
point. 

The Schrodinger equation is considered as the fundamental equation of 
nonrelativistic quantum mechanics in the same sense in which Newton's 
force law is the fundamental equation of motion of newtonian mechanics. 
Notice however that we have now a wave equation for a function ^' which is 
not directly measurable. 

Once the potential energy is given, one can solve the Schrodinger equa- 
tion for implying the knowledge of the probability density | ^ p as a 
function of x, y, z, t. In many cases of interest, the potential energy does not 
depend on time. Then, the Schrodinger equation simplifies considerably. 
Notice, for example, that for a ID free particle the wave function can be 
written 

= ^{x)e-^'^/^^* , (4) 
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i.e., ^{x,t) is the product of a time-dependent phase e^^*^/'*)* and a sta- 
tionary wave function '4>{x). 

In the general case, the stationary Schrodinger equation can be solved, 
under the aforementioned requirements, only for certain values of the energy 
E. This is not a mathematical difficulty, but merely a fundamental physi- 
cal feature. To solve the Schrodinger equation for a given system means 
to get the wave function ip, as a solution for which certain physical bound- 
ary condition hold and, in addition, as already mentioned, it is continuous 
together with its first derivative everywhere in space, is finite, and single 
valued. Thus, the quantization of energy occurs as a natural theoretical 
element in wave mechanics, whereas in practice as a universal phenomenon, 
characteristic for all stable microscopic systems. 



Schrodinger equation for the hydrogen atom 



In this section, we shall apply the Schrodinger equation to the hydrogen 

atom, about which one knows that it is formed of a positive nucleus/proton 
of charge +e and an electron of charge -e. The latter, being 1836 times 
smaller in mass than the proton, is by far more dynamic. 

If the interaction between two particles is of the type u{r) = u{\ r\ — f2 |), 
the problem of the motion is reduced both classically and quantum to the 
motion of a single particle in a field of spherical symmetry. Indeed, the 
Lagrangian 



L = -niir^ H — m2r2 



is transformed, using 



and 



in the Lagrangian 



where 



and 



R 



r = n — r2 



miri -I- m2r2 
mi + m2 



L = ]^MB? + ^fj,r^ - u{r) , 



M = mi + m2 



mim2 
mi + m2 



(5) 
(6) 

(7) 

(8) 

(9) 
(10) 
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On the other hand, the momentum is introduced through the Lagrange 
formula 

P=^ = Mk (11) 
dR 

and 

- dL 

p = —^ = mr , (12) 
or 

that allows to write the classical Hamilton function in the form 

Thus, one can obtain the hamiltonian operator for the corresponding 
quantum problem with commutators of the type 

[Pi,Pk] = -in5ik (14) 

and 

\PhPk] = -ifiSik . (15) 
These commutators implies a Hamiltonian operator of the form 

which is fundamental for the study of the hydrogen atom by means of the 
stationary Schrodinger equation 

HiIj = EiI; . (17) 

This form does not include relativistic effects, i.e., electron velocities close 

to the velocity of light in vacuum. 

The potential energy u{r) is the electrostatic one 

u = -- (18) 

47reor 

There arc two possibilities. The first is to express u as a function of the 
cartesian coordinates x,y,z, substituing r by a/x^ + j/^ + z"^. The second 
is to write the Schrodinger equation in spherical polar coordinates r,9,(f). 
Because of the obvious spherical symmetry of this case, we shall deal with the 
latter approach, which leads to considerable mathematical simplifications. 
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In spherical coordinates, the Schrodinger equation reads 



1 d f 29tp\ , 1 5 /^_- , ^ 5^V' , 2m 



r2 9r V"" dr)^ sin 9 86 89 ) ^ sin^ 9 9</)2 + ""^^ " ° 

(19) 

Substituing (18), and multiplying the whole equation by r^sin^^, one gets 



sin^^^ 
8r 



')8ip\ . ^8 f . ^8Tp\ 8'^'ib 2mr'^s\T?9 ( J\ , 



(20) 

This equation is a partial differential equation for the electron wavefunction 
^{r, 9, (p) 'within' the atomic hydrogen. Together with the various conditions 
that the wavefunction ■ip{r, 9, cf)) should fulfill [for example, '4){r, 9, (j)) should 
have a unique value at any spatial point {r,9,(j))], this equation specifies 
in a complete manner the behavior of the hydrogen electron. To sec the 
explicit behavior, we shall solve eq. 20 for 'ip[r,9,(f)) and we shall interpret 
appropriately the obtained results. 



Separation of variables in spherical coordinates 

The real usefulness of writing the hydrogen Schrodinger equation in spherical 
coordinates consists in the easy way of achieving the separation procedure 
in three independent equations, each of them being one-dimensional. The 
separation procedure is to seek the solutions for which the wavefunction 
tp{r, 9, (j)) has the form of a product of three functions, each of one of the 
three spherical variables, namely R{r), depending only on r; Q{9) depending 
only on 6, and that depends only on ^. This is quite similar to the 

separation of the Laplace equation. Thus 

i,{r,9,<l>)=R{r)@{9mct>) . (21) 

The R{r) function describes the differential variation of the electron wave- 
function ip along the vector radius coming out from the nucleus, with 9 and 
(f) assumed to be constant. The differential variation of -0 with the polar 
angle 9 along a meridian of an arbitrary sphere centered in the nucleus is 
described only by the function 6(6*) for constant r and (j). Finally, the func- 
tion $(^) describes how ^ varies with the azimuthal angle along a parallel 
of an arbitrary sphere centered at the nucleus, under the conditions that r 
and 9 are kept constant. 
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Using = i?G$, one can see that 



Obviously, the same type of formulas are maintained for the unmixed higher- 
order derivatives. Subtituting them in eq. 20, and after deviding by i?G$, 

we get 



sin^g d 
R dr 



{ 2dR\ sine d f . ,de\ 1 2mr^ sin'^ 6 f \ ^ 

['^r^deH^r^w^^^ iw+^j=« 

(25) 

The third term of this equation is a function of the angle only, while the 
other two terms are functions of r and 9. We rewrite now the previous 
equation in the form 

sm^e d { odR\ sine d / de\ 2mr^sin'^ef \ 18'^^ 

+— — — sin^— + + 



R dr\ drj e de\ 89 J \ATT€Qr J ^ dq 

(26) 

This equation can be correct only if the two sides are equal to the same 
constant, because they are functions of different variables. It is convenient 
to denote this (separation) constant by mf. The differential equation for 
the $ function is 

1 2 

If one substitutes rnf in the right hand side of eq. 26 and devides the 
resulting equation by sin^ 9, after regrouping the terms, the fllowing result 
is obtained 

I d f 2dR\ , 2mr2 f e'^ mf 1 d f . ^de\ 



Rdr\! dr)"" l^47reor+^j sin^ @ sin 9 d9 V'"" d9 ) - 

(28) 

Once again, we end up with an equation in which different variables occur 
in the two sides, thus forcing at equating of both sides to the same constant. 
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For reasons that will become clear later on, we shall denote this constant by 
l{l + 1). The equations for the functions Q{0) and R{r) reads 



mt 



d 



and 



,dR 



2mr^ 



+ E 



dr J fi y 47reor 
The equations 27, 29 and 30 are usually written in the form 



+ $ = 



1 d 

^^de 



sm 



dr V dr 



de ) 

Ti^ l47reor 



+ E 








(29) 
(30) 

(31) 
(32) 
(33) 



Each of these equations is an ordinary differential equation for a function 
of a single variable. In this way, the Schrodinger equation for the hydrogen 

electron, which initially was a partial differential equation for a function i]j of 
three variables, got a simple form of three ID ordinary differential equations 
for unknown functions of one variable. 



Interpreting the separation constants: the quantum 
numbers 



The solution for the azimuthal part 

Eq. 31 is readily solved leading to the following solution 



$(^) = A^e 



,imi() 



(34) 



where is the integration constant. One of the conditions that any wave- 
functions should fulfill is to have a unique value for any point in space. This 

applies to $ as a component of the full wavefunction iIj. One should notice 
that (f> and + 27r arc identical in the same meridional plane. Therefore, 



one should have '1>(0) = $(0 + 27r), i.e., A^e 



imiq 



^ im,(</'+27r)_ This can 



be fulfilled only if mi is zero or a positiv or negative integer (±1, ±2, ±3, ...). 
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mi is known as the magnetic quantum number of the atomic electron and 
is related to the direction of the projection of the orbital momentum L^. It 
comes into play whenever the effects of axial magnetic fields on the electron 
may show up. There is also a deep connection between mi and the orbital 
quantum number I, which in turn determines the modulus of the orbital 
momentum of the electron. 

The interpretation of the orbital number I does not miss some problems. 
Let us examine eq. 33 that corresponds to the radial wavefunction R{r). 
This equation rules only the radial motion of the electron, i.e., with the 
relative distance with respect to the nucleus along some guiding ellipses. 
However, the total energy of the electron E is also present. This energy 
includes the kineticelctron energy in its orbital motion that is not related 
to the radial motion. This contradiction can be eliminated by the following 
argument. The kinetic energy T has two parts: Tradiat due to the radial 
oscillatory motion and Torutah which is due to the closed orbital motion. 
The potential energy V of the electron is the electrostatic energy. Therefore, 
its total energy is 

E = Tradial + Tgrbital ~ ~, • (35) 

Substituting this expression of E in eq. 33 we get with some regrouping of 
the terms 



Id/ ^dR\ 2m 
dr \ dr J 



n^iji + 1) 

'-radial "r ^orbital r, o 



R = Q . (36) 



If the last two terms in parentheses compansates between themselves, we 
get a differential equation for the pure radial motion. Thus, we impose the 
condition 

_ nH{i + 1) 

J- orbital - 2j„^2 • 

However, the orbital kinetic energy of the electron is 

Torbital = ^niv^j.^n^i (38) 
and since the orbital momentum of the electron L is 

L = mVorHtair , (39) 
we can express the orbital kinetic energy in the form 

'^orbital = • (40) 
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Therefore, we have 



L2 



n'^i{i + 1) 



(41) 



2mr2 



and consequently 



L = ^Jl{l + l)h . 



(42) 



The interpretation of this result is that since the orbital quantum number I 
is constrained to take the values I = 0,1, 2,..., (n — 1), the electron can only 
have orbital momenta L specified by means of eq. 42. As in the case of the 
total energy E, the angular momentum is conserved and gets quantized. Its 
natural unit in quantum mechanics is ^ = /i/27r = 1.054 x 10~'^^ J.s. 

In the macroscopic planetary motion (putting aside the many-body fea- 
tures), the orbital quantum number is so large that any direct experimental 
detection is impossible. For example, an electron with I = 2 has an angular 
momentum L = 2.6 X 10~^^ J.s., whereas the terrestrial angular momentum 
is 2.7 X 10^° J.S.! 

A common notation for the angular momentum states is by means of 
the letter s for Z = 0, p for Z = 1, d for Z = 2, and so on. This alphabetic 
code comes from the empirical spectroscopic classification in terms of the 
so-called series, which was in use before the advent of quantum mechanics. 

The combination of the principal quantum number with the latter cor- 
responding to the angular momentum is another frequently used notation 
in atomic and molecular physics.. For example, a state for which n = 2 and 
Z = is a state 2s, while a state n = 4 and Z = 2 is a state id. 

On the other hand, for the interpretation of the magnetic quantum num- 
ber, we shall take into account, as we did for the linear momentum, that 
the orbital momentum is a vector operator and therefore one has to specify 
its direction, sense, and modulus. L, being a vector product, is perpendic- 
ular on the plane of rotation. The geometric rules of the vectorial products 
still hold, in particular the rule of the right hand: its direction and sense 
arc given by the right thumb whenever the other four fingers point at the 
direction of rotation. 

But what significance can be associated to a direction and sense in the 
limited space of the atomic hydrogen ? The answer may be quick if we 
think that the rotating electron is nothing but a one-electron loop current 
that considered as a magnetic dipole has a corresponding magnetic field. 
Consequently, an atomic electron will always interact with an applied mag- 
netic B. The magnetic quantum number mi specifies the spatial direction 
of L, which is determined by the component of L along the direction of the 
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external magnetic field. This effect is commonly known as tlie quantization 
of the space in a magnetic field. 

If we choose the direction of the magnetic field as the z axis, the com- 
ponent of L along this direction is 

= mih . (43) 

The possible values of nii for a given value of I, go from +/ to —I, passing 
through zero, so that there arc 2/ + 1 possible orientations of the angular 
momentum L in a magnetic field. When i = 0, can be only zero; when! = 
1, Lz can be h, 0, or —h; when I = 2, takes only one of the values 2h, U, 
0, — ^, or —2^, and so forth. It is worth mentioning that L cannot be put 
exactly parallel or anti-parallel to i?, because is always smaller than the 
modulus ^/l(^+T)h of the total orbital momentum. 

The spatial quantization of the orbital momentum for the hydrogen atom 
is shown in fig. 6.1 in a particular case. 




Fig. 6.1: The spatial quantization of the electron angular momentum for states 

1 = 2, L = Ven. 

One should consider the atom/electron characterized by a given mi as 
having the orientation of its angular momentum L determined relative to 
the external applied magnetic field. 

In the absence of the external magnetic field, the direction of the z axis 
is fully arbitrary. Therefore, the component of L in any arbitrary chosen 
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direction is mih; the external magnetic field offers a preferred reference 
direction from the experimental viewpoint. 

Why is quantized only the component ? The answer is related to 
the fact that L cannot be put along a direction in an arbitrary way. Its 
'vectorial arrow' moves always along a cone centered on the quantization 
axis such that its projection is mih. The reason why such a phenomenon 
occurs is due to the uncertainty principle. If L would be fixed in space, in 
such a way that L^, Ly and would have well-defined values, the electron 
would have to be confined to a well-defined plane. For example, if L would 
be fixed along the z direction, the electron tends to maintain itself in the 
plane xy (fig. 6.2a). 



Fig. 6.2: The uncertainty principle forbids a fixed direction in space of the 

angular momentum. 

This can occur only in the case in which the component Pz of the elec- 
tron momentum is 'infinitely' uncertain. This is however impossible if the 
electron is part of the hydrogen atom. But since in reality just the compo- 
nent Lz of L together with have well-defined values and | L | > | | , the 
electron is not constrained to a single plane (fig. 6.2b). If this would be the 
case, an uncertainty would exist in the coordinate z of the electron. The 
direction of L changes continuously (see fig. 6.3), so that the mean values 
of Lx and Ly are zero, although keeps all the time its value mih. 




fa 
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Fig. 6.3: The angular momentum displays a constant precession around the z axis. 

The solution for $ should also fulfill the normalization condition given 
by eq. 2. Thus, we have 



and substituting one gets 



/■2-7r 

/ I $ p # = 1 

h 

Ald4> = i. 



(44) 



(45) 



It follows that = l/\/2Tr, and thefore the normalized $ is 

1 



imicj) 



27r 



(46) 



Solution for the polar part 

The solution of the 6(^) equation is more complicated. It is expressed in 
terms of the associated Legendre polynomials 



Prix) = (-i)-<(i-x2)-</2^p,(x) = (-1)-' 



2Hl dx'^i+^ 



(47) 
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Their orthogonality relationship is 



' [P^ (cos9)]''dcos9 = 7^^^ . (48) 

1 ' ^ 2l + l{l-mi)\ ^ ' 



For the case of quantum mechanics, Q{9) is given by the normalized associ- 
ated Legendre polynomials. Thus, if 

Q{9) = AePj;^\cose) , (49) 

then the normalization condition is 

nl 



/_ 



Al[Pp{cos9)fdcos9 = 1 . (50) 



Therefore, the normalization constant for the polar part is given by 



21 + l il-miy. 

and consequently, the function @{9) already normalized reads 



For our purposes here, the most important property of these functions 
is that they exist only when the constant I is an integer number greater or 
at least equal to | m; |, which is the absolute value of mi. This condition 
can be written in the form of the set of values available for mi 

mi = 0,±l,±2,...,±l . (53) 

Unification of the azimuthal and polar parts: spherical har- 
monics 

The solutions of the azimuthal and polar parts can be unified within spher- 
ical harmonics functions that depend on both (p and 0. This simplifies the 
algebraic manipulations of the full wave functions il){r,9,4'). Spherical har- 
monics are introduced as follows 



Yr{9,<t>) = (-^rf-^J^,Pri^os9)e^'^^f . (54) 
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The supplementary factor (—1)"^' does not produce any problem because 
the Schrodinger equation is linear and homogeneous. This factor is added 
for the sake of convenience in angular momentum studies. It is known as the 
Condon-Shortley phase factor and its effect is to introduce an alternance of 
the signs ± for the spherical harmonics. 



Solution for the radial part 

The solution for the radial part R{r) of the wave function ip of the hydrogen 
atom is somewhat more complicated. It is here where significant differences 
with respect to the electrostatic Laplace equation do occur. The final result 
is expressed analytically in terms of the associated Laguerre polynomials 
(Schrodinger 1926). The radial equation can be solved in exact way only 
when E is positive or for one of the following negative values £'„ (in which 
cases, the electron is in a bound stationary state within atomic hydrogen) 



327r2e2n2 W 

where n is an integer number called the principal quantum number. It 
gives the quantization of the electron energy in the hydrogen atom. This 
discrete atomic spectrum has been first obtained in 1913 by Bohr using semi- 
empirical quantization methods and next by Pauli and Schrodinger almost 
simultaneously in 1926. 

Another condition that should be satisfied to solve the radial equation is 
that n have to be strictly bigger than /. Its lowest value is Z + 1 for a givem 
I. Vice versa, the condition on I is 

Z = 0,l,2,...,(n-1) (56) 

for given n. 

The radial equation can be written in the form 



nd'^R „ dR 
r -r^ + ^r— + 
cir^ ar 



2mE o 2me 



,2 



-r 



-m+i) 



R = Q , (57) 



Dividing by and using the substitution x(r) = rR to eliminate the first 
derivative one gets the standard form of the radial Schrodinger equation 
displaying the effective potential U{r) = — const /r + 1(1 + l)/r'^ (actually, 
electrostatic potential plus quantized centrifugal barrier). These are nec- 
essary mathematical steps in order to discuss a new boundary condition. 
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since the spectrum is obtained by means of the R equation. The difference 
between a radial Schrodinger equation and a fuU-Une one is that a suppU- 
mentary boundary condition should be imposed at the origin (r = 0). The 
coulombian potential belongs to a class of potentials that arc called weak 
singular for which lim^^o = U{r)r'^ = 0. In these cases, one tries solutions 
of the type x oc r'^, implying u{u — 1) = + 1), so that the solutions are 
1/1=1 + 1 and 1^2 = —I, just as in electrostatics. The negative solution is 
eliminated for / ^ because it leads to a divergent normalization constant, 
nor did it respect the normalization at the delta function for the continuous 
part of the spectrum. On the other hand, the particular case f 2 = is elmin- 
inated because the mean kinetic energy is not finite. The final conclusion is 
that x(0) = for any I. 

Going back to the analysis of the radial equation for R, first thing to do 
is to write it in nondimensional variables. This is performed by noticing that 
the only space and time scales that one can form on combining the three 
fundamental constants entering this problem, namely e^, rrif, and h are the 
Bohr radius qq = H^/me'^ = 0.529-10-^ cm. and to = h^jme^ = 0.24210-^^ 
sec, usually known as atomic units. Employing these units, one gets 



d^R 2 dR 



H T- + 



2E + - 



2 ^(^ + 1) 



i? = , (58) 



where we are especially interested in the discrete part of the spectrum (E < 
0). The notations n = l/\/—E and p = 2r/n leads us to 



d^R 2 dR 



+ — r- + 



dp^ p dp 



n 1 1(1 + 1) 



P 4: p2 



R = . (59) 



For p 00, this equation reduces to = -f , having solutions R oc e"^^/^. 
Because of the normalization condition only the decaying exponential is 
acceptable. On the other hand, the asymptotics at zero, as we already com- 
mented on, should be i? oc pK Therefore, we can write i? as a product 
of three radial functions R = p^e^''^'^F{p), of which the first two give the 
asymptotic behaviors, whereas the third is the radial function in the inter- 
mediate region. The latter function is of most interest because its features 
determine the energy spectrum. The equation for F is 

d"^ F dF 
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This is a particular case of confluent hyper geometric equation for which 
the two 'hyper'geometric parameters depend on the pair of quantum num- 
bers n, I. It can be identified as the equation for the associated Laguerre 
polynomials L'^^^{p). Thus, the normalized form of R is 

where the following Laguerre normalization condition has been used 

We have now the solutions of all the equations depending on a single 
variable and therefore we can build the wave function for any electronic 
state of the hydrogen atom. The full wave function reads 

ilj{r, e, 4>) = MH{arye-''''/^L^^+l{ar)Pl^' {cose)e'"'"^ , (63) 

where Mh = -^\[^^^^^^^ ^nd a = 2/nao. 

Using the spherical harmonics, the solution is written as follows 

^irA<l^) = -;|y^f^^("-)'^""'^^'^n|/K)l^r(e,'/') • (64) 

The latter formula may be considered as the final result for the Schrodinger 
solution of the hydrogen atom for any stationary electron state. Indeed, one 
can see explicitly both the asmptotic dependence and the two orthogonal 
and complete sets of functions, i.e., the associated Laguerre polynomials 
and the spherical harmonics that correspond to this particular case of lin- 
ear partial second-order differential equation. The parabolic coordinates 
= r(l — cos 6*), rj = r(l -|- cos^), (p = ^]) are another coordinate system in 
which the Schrodinger hydrogen equation is separable (E. Schrodinger, Ann. 
Physik 80, 437, 1926; P.S. Epstein, Phys. Rev. 28, 695, 1926; I. Waller, 
Zf. Physik 38, 635, 1926). The final solution in this case is expressed as the 
product of factors of asymptotic nature, azimuthal harmonics, and two sets 
of associate Laguerre polynomials in the variables ^ and rj, respectively. The 
energy spectrum (— 1/n^) and the degeneracy (n^) of course do not depend 
on the coordinate system. 
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Electronic probability density 



In the Bohr model of the hydrogen atom, the electron rotates around the nu- 
cleus on circular or elliptic trajectories. It is possible to think of appropriate 
experiments allowing to "see" that the electron moves within experimental 
errors at the predicted radii r = n^oo (where n is the principal quantum 
number labeling the orbit and oq = 0.53 A is the Bohr radius) in the equa- 
torial plane 6 = 90°, whereas the azimuthal angle may vary according to the 
specific experimental conditions. 

The more rigorous quantum theory changes the conclusions of the Bohr 
model in at least two important aspects. First, one cannot speak about exact 
values of r, 6, cp, but only of relative probabilities to find the electron within 
an infinitesimal given region of space. This feature is a consequence of the 
wave nature of the electron. Secondly, the electron does not move around 
the nucleus in the classical conventional way because the probability density 
I ■0 P does not depend on time but can vary substantially as a function of 
the relative position of the infinitesimal region. 

The hydrogenic electron wave function il^'isip = R@^, where R = Rni{r) 
describes the way ^ changes with r when the principal and orbital quantum 
numbers have the values n and I, respectively. @ = Qimi (^) describes in 
turn how ijj varies with 9 when the orbital and magnetic quantum numbers 
have the values I and m/ , respectively. Finally, $ = ^rni {^) gives the change 
of with (j) when the magnetic quantum number has the value mi. The 
probability density | V' P can be written 

I V P=| R n Pi * ? ■ (65) 

Notice that the probability density | $ which measures the possibility 
to find the electron at a given azimuthal angle (f), is a constant (does not 
depend on (f>). Therefore, the electronic probability density is symmetric 
with respect to the z axis and independent on the magnetic substates (at 
least until an external magnetic field is applied) . Consequently, the electron 
has an equal probability to be found in any azimuthal direction. The radial 
part R of the wave function, contrary to not only varies with r, but it 
does it differently for any different combination of quantum numbers n and 
I. Fig. 6.4 shows plots of as a function of r for the states Is, 2s, and 
2p. R is maximum at the center of the nucleus (r = 0) for all the s states, 
whereas it is zero at r = for all the states of nonzero angular momentum. 
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Fig. 6.4: Approximate plots of the radial functions Ris, i?2s, R2p] (ao — 0.53 A). 

A 2 2 
dP/dr=47ir R 




Fig. 6.5: Probability density of finding the hydrogenic electron between r and 
r + dr with respect to the nucleus for the states Is, 2s, 2p. 

The electronic probability density at the point r, 9, (p is proportional to 
ip P, but the real probability in the infinitesimal volume element dV is 
■0 p dV. In spherical coordinates 

dV = r'^ sin edrd9d(l) , (66) 
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and since B and $ are normalized functions, the real numerical probability 
P{r)dr to find the electron at a relative distance with respect to the nucleus 
between r and r + is 

P{r)dr = r"^ \ R\'^ dr r \ smOdO /^'' | $ p 

Jo Jo 

= I p dr (67) 

P(r) is displayed in fig. 6.5 for the same states for which the radial func- 
tions R appear in fig. 6.4. In principle, the curves are quite different. We 
immediately see that P{r) is not maximal in the nucleus for the states s, as 
happens for R. Instead, their maxima are encountered at a finite distance 
from the nucleus. The most probable value of r for a Is electron is exactly 

00, the Bohr radius. However, the mean value of r for a Is electron is l.Sao- 
At first sight this might look strange, because the energy levels are the same 
both in quantum mechanics and in Bohr's model. This apparent unmatch- 
ing is eliminitated if one takes into account that the electron energy depends 
on 1/r and not on r, and the mean value of 1/r for a Is electron is exactly 
I/oq. 

The function 9 varies with the polar angle 9 for all the quantum num- 
bers I and mi, unless I = mi = 0, which are the s states. The probability 
density | 6 ^ for a s state is a constant (1/2). This means that since | ^ P 
is also a constant, the electronic probability density | P li^s the same 
value for a given r value, not depending on the direction. In other states, 
the electrons present an angular behavior that in many cases may be quite 
complicated. This can be seen in fig. 6. 5, where the electronic probability 
densities for different atomic states are displayed as a function of r and 6. 
(The plotted term is | V' P and not | V" P dV). Because | -0 p is independent 
of (j), a three-dimensional representation of | -0 p can be obtained by rotat- 
ing a particular representation around a vertical axis. This can prove that 
the probability densities for the s states have spherical symmetry, while all 
the other states do not possess it. In this way, one can get more or less 
pronounced lobes of characteristic forms depending on state. These lobes 
are quite important in chemistry for specifying the atomic interaction in the 
molecular bulk. 

6N. Note: 

1. In 1933, E. Schrodinger has been awarded the Nobel Prize in Physics (to- 
gether with Dirac) for the "discovery of new productive forms of atomic the- 
ory" . Schrodinger wrote a remarkable series of four papers "Quantisierung 
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als Eigenwertproblem" ["Quantization as an eigenvalue problem"] (I-IV, re- 
ceived by Annalen der Physik on 27 January, 23 February, 10 May and 21 
June 1926, respectively). 



6P. Problems 



Problem 6.1 - Obtain the formulas for the stable orbits and the energy 
levels of the electron in the atomic hydrogen using only arguments based on 
the de Broglie wavelength associated to the electron and the empirical value 
5.3 • 10~^^ m for the Bohr radius. 

Solution: The electron wavelength is given by A = whereas if we 
equate the electric force and the centripetal force = we obtain 

the electron 'velocity' v = -?= .Thus, the wavelength of the electron 



is A = ^y ^m^ - If we now use the value 5.3 x 10~^^m for the radius 

r of the electron orbit, we can see that the wavelength of the electron is 
A = 33 X 10~^^ m. But this is exactly the same value as of the circumference 
of the orbit, 27rr = 33 x 10~^^ m. One may say that the electron orbit 
in the atomic hydrogen corresponds to a wave "closing into itself (i.e., 
stationary). This fact can be compared to the vibrations of a metallic ring. If 
the wavelengths are multiples of the circumference, the ring goes on with its 
vibrations for a long time with very small dissipation If, on the other hand, 
the number of wavelengths making a circumference is not an integer, the 
interference of the waves is negative and they dissapear in a short period of 
time. One may say that the electron will rotate around the nucleus without 
radiating its energy for an infinite time as far as its orbit contains an integer 
number of de Broglie wavelengths. Thus, the stability/stationary condition 
is 



where r„ is the radius of the electron orbit containing n wavelengths. Sub- 
stituting A, we have 





and therefore the stationary electron orbits are 




irme- 
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To get the energy levels, we use E = T + V and substituting the kinetic 
and potential energies leads to 

^ 1 2 e2 
E = -mv 



or equivalently 



e2 



Plugging the value of r„ into the latter equation, we get 

me^ / 1 

En = — 



Problem 6.2 - Unsold's theorem tells that for any value of the orbital 
number I, the probability densities, summed over all possible substates, from 
mi = —I to mi = +1 give a constant that is independent of the angles 6 and 
0, i.e. 

+/ 

E I ^im Pi ?=ct. 

mi=-l 

This theorem shows that any atom or ion with closed (occupied) sublevels 
has a spherically-symmetric charge distribution. Check Unsold's theorem for 
^ = 0, i = 1, and I = 2. 

Solution: For I = 0, 9oo = l/V2 and $o = l/V^, so that 

I eo,o I'l $0 1'= ■ 

' 47r 

For Z = 1, we have 

+1 

E I ®imi I'l M I'l $-1 1' + 1 ei,o Pi $0 1' + 1 ei,i |2| $1 



|2 

mi=-l 



On the other hand, the wave functions are given by 01 _i = iV3/2)sin9, 
= (l/^/2^)e-^<^, 01,0 = {V6/2)cose, $o = 1/^2^, 61,1 = {V3/2)sine, 
$1 = (l/\/27r)e*'''' , which plugged into the previous equation give 



m;=-l 
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and again we've got a constant. 
For Z = 2, we have 

+2 

E I ®lmi I'l M ©2,-2 I'l $-2 I'l 62,-1 I'l $-1 I' 

mi=-2 

+ I ©2,0 I'l ^0 I' + I ©2,1 I'l ^1 I' + I 62,2 I'l ^2 I' , 

and the wave functions are 62,-2 = (\/T5/4) sin^^, <l>_2 = (l/\/27r)e~^*'^, 
62,-1 = {Vl5/2)sinecose, $_i' = (l/v^je-*"^, 62,0 = {VW/A){3cos^e - 
1), $0 = l/\/2^, 62,1 = (^/T5/2)sin^cos^, $1 = (l/^/2^)e^<^, 62,2 = 
{■\/l5 / 4) sin^ 9, $2 = {1/V^)e'^'''^, Plugging them into the previous equa- 
tion give 




which again fulfills Unsold's theorem. 



Problem 6.3 - The probability to find an atomic electron whose radial wave 
functions is that of the ground state Rio{r) outside a sphere of Bohr radius 
ao centered on the nucleus is 



roc 



Rio{r) P r'^dr . 



Obtain the probability to find the electron in the ground state at a distance 
from the nucleus bigger than cq- 

Solution: The radial wave function corresponding to the ground state is 

2 



Rio{r) 



3/2 



-r/ao 



Substituting it in the integral, we get \ R{r) p r'^dr = r^e ^r/ao^^ 



roc 
Jao 



R{r) 1^ r^dr 



2 2 4 



ao 



This leads us to 



roc 
Jao 



R{r) p r^dr 



68% !! , 



which is the result asked for in this problem. 
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7. QUANTUM SCATTERING 

Introduction 

One usually begins the quantum theory of scattering by referring to results 
already known from the classical scattering in central fields with some sim- 
plifying assumptions helping to avoid unnecessary calculations in getting 
basic results. It is generally known that studying scatterings in the labo- 
ratory provides information on the distribution of matter in the target and 
other details of the interaction between the incident beam and the target. 
The hypotheses that we shall assume correct in the following are 

i) The particles are spinless. This, of course, does not mean that spin 
effects are not important in quantum scatterings. 

ii) We shall study only elastic scattering for which the internal structure 
of the particles is not taken into account. 

iii) The target is sufficiently thin to neglect multiple scatterings. 

iv) The interactions are described by a potential that depends only on 
the relative distance between the particles (central potential). 

These hypotheses eliminate some quantum effects that are merely details. 
They also represent conditions for getting the quantum analogs of basic 
classical results. We now define 

where is the solid angle infinitesimal element, Iq is the number of inci- 
dent particles per unit transverse area, and IdQ is the number of scattered 
particles in the solid angle element. 

Employing these well-known concepts, together with the asymptotic no- 
tion of impact parameter b associated to each classical incident particle, one 
gets in classical mechanics the following important formula 

da b db 

dn ~ ih^'d^' ■ ^ ' 

If one wants to study the scattering phenomenology in quantum ter- 
minology, one should investigate the time evolution of a 'scattering' wave 
packet. Let Fi be the flux of incident particles, i.e., the number of parti- 
cles per unit of time passing through the unit of transverse surface onto the 
propagation axis. An appropriate detector configuration is usually placed 
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far away from the effective interaction region, 'seeing' a solid angle dO. of 
that region. In general, the number of particles dn/dt scattered per unit of 
time in dQ in the direction [9, f) is detected. 




Inc . flux F i 

Fig. 7.1 



dn/dt is proportional to dO, and Fj. Let us call a{9,ip) the coefficient of 
proportionality between dn and FidQ,: 

dn = a{e, (p)Fidn , (3) 

which is by definition the differential cross section. 

The number of particles per unit of time reaching the detector is equal 
to the number of particles crossing the surface a{6, (p)dQ., which is perpen- 
dicular to the beam axis. The total section is by definition 

a = J a{e, ip)dn . (4) 

To further simplify the calculation, we choose the z axis along the inci- 
dent beam direction. 

On the negative side of the axis, for large negative t, the particle is prac- 
tically free: it is not affected by V{r) and its state can be represented by 
plane waves. Therefore, the wave function contains terms of the form e^^^, 
where k is the constant ocurring in the Helmholtz equation. By analogy 
with optics, the form of the scattered wave is 

pikr 

m = — . (5) 
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Indeed 

(V^ + k'^)e'''^ 7^ (6) 

and 

(V^ + A;^)— =0 (7) 



r 

for r > ro, where vq is any positive number. 

We assume that the motion of the particle is described by the Hamilto- 
nian 

2 

H = ^ + V = Ho + V . (8) 

V is different of zero only in a small neighbourhood close to the origin. 
A wave packet at t = can be written 

V^(r, 0) = [ ^(k) exp[zk • (r - ro)]d=*k , (9) 

(27r)2 J 

where V' is a function that is nonzero in a 'width' Ak centered on ko- We 

also assume that ko is antiparallcl to Tq. In order to see quantitatively 
what happens to the wave packet when scatters the target, one can use 
the expansion of V'(r, 0) in the eigenfunctions V'n(r) of i?, i.e., ip{r,^) = 
SnCnV'n(r)- Thus, the wave packet at time t is 

ip{Y,t) = ^c„(^n(r)exp(-^£;ni) • (10) 

n 

This is an eigenfunction of the operator Hq, not of H, but we can sub- 
stitute these eigenfunctions by eigenfunctions of H, which we denote by 
V'fe^^ (r) . The asymptotic form of the latter is 

V'i+)(r)~e'''-'- + f(r)^ , (11) 
where, as usually p = ^k and E = 

This corresponds to a plane wave of the incident beam type and a di- 
vergent spherical wave resulting from the interaction between the incident 
beam and the target. One can expand ip{r, 0) in plane waves and V'fc(r) 

^(r, 0) = y (^(k) exp(-ik • ro)M^)d^k , (12) 

where hu = k^^. The divergent spherical wave does not contribute to the 
initial wave packet because it is an additive part. 
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Scattering of a wave packet 



Any wave is dispersed during its propagation. This is why one cannot ignore 
the effect of the divergent wave from this viewpoint. One can make use of 
the following trick 



= 2^^' = 2^t^° + " = ^^^^ • k - + (k - ko)^] , (13) 

pentru a neglija ultimul termen in parantczc. Substituting uj in tp, we ask 
that ^(k - ko)^T < 1, where T ~ Therefore 

« 1 . (14) 

This condition tells us that the wave packet docs not disperse significantly 
even when it moves over amacroscopic distance tq. 

Choosing the direction of the vector k of the incident wave along one of 
the three cartesian directions (we use the z one), we can write in spherical 
coordinates the following important formula 

V^fc(r,g,y) -e^^^ + ^^^-y"^ . 

Since the Hamiltonian H, up to now not considered as an operator (the 
class of the results presented are the same both at the classical and quantum 
level), is invariant under z rotations, we can choose boundary conditions of 
spherical symmetry too. Thus 

This type of functions are known as scattering wave functions. The coeffi- 
cient f{6) of the spherical wave is known as the scattering amplitude. It is 
a basic concept in the formal theory of quantum scatterings. 



Probability amplitude in scattering 

We write the Schrodinger equation as follows 

ih^ = -^V''^ + Vir,t)^ . (15) 

Recall that the expression 

P(r,t) = ^*(r,t)^(r,t) = |V'(r,t)p (16) 
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can be interpreted, cf. Max Born, as a probability density under normaliza- 
tion conditions of the type 



/I 



|V(r,t)|^d^r = l . (17) 

This normalization integral should be time independent. This can be 
noted by writing 

and from Schrodinger's equation 
one gets 



2m Jn 2m Jci 



= ^ [ [r^^-{^r)^]ndA, (20) 

2m J A 

where the Green theorem has been used to evaluate the volume integral. 
dA is the infinitesimal surface clement on the boundary of the integration 
region and [ denotes the component along the normal direction to the 
surface element dA. 
Defining 

S(r,t) = ^[rVV-(VV'*)V] , (21) 

we get 

7 = 4/ P{r,t)d^r = - / V • SdV = - / SndA , (22) 
at Jn Jn J A 

for well-bahaved wave packets (not funny asymptotically) so that the nor- 
malization integral converges. The surface integral is zero when Q, covers the 
whole space. One can prove (see P. Dennery & A. Krzywicki, Mathematical 
methods for physicists) that the surface integral is zero. Therefore, the nor- 
malization integral is constant in time and the initial condition holds. From 
the same equation for S, we get 

^^ + V-S(r,t) = 0, (23) 
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which is the continuity equation for the density flux P and the current 
density S in the absence of any type of sources or sinks. If we interpret ^ V 
as a sort of velocity 'operator' (as for time, it is difficult to speak rigorously 
about a velocity operator in quantum mechanics!), then 

S(r,t) = Re(Tp*—Vilj) . (24) 

To calculate the quantum current density for a scattering wave function 
is a tricky and inspiring (not illustrative) exercise! The final result is jV = 
-^|/(0)P, where the direction 9 = should not be included. 



Green's function in scattering theory 

Another way of writing the Schrodinger equation at hand is (~^V^+y)V' = 
Ei;, or (V2 + A;2)-0 = [7-0, where k'^ = §i J7 = 

It follows that it is more convenient to put this equation in an integral 
form. This can be done if we consider [/■!/; in the right hand side of the 
equation as a inhomogeneity. This allows to build the solution by means of 
Green's function (integral kernel), which, by definition, is the solution of 

(V2 + A;2)G(r,r') = (5(r-r') . (25) 

One can write now the Schrodinger solution as the sum of the homogeneous 
equation and the inhomogeneous one of Green's type 

'0(r) = A(r) - J G(r,r )U(r )'0(r )d^r' . (26) 

We seek now a G function in the form of a product of linear independent 
functions, for example, plane waves 

G(r,r' = J ^(q)e^^-('-'"')dg . (27) 



Using eq. 25, we have 

J A(q)(fc2 _ g2)eiq-(r-r')^^ ^ g^j. _ ^'^ ^ (28) 

which turns in an identity if 

A(q) = (27r)-3(fe2_52^-i _ (29) 
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Thus 

where R = |r — r'|. After performing a calculation of complex variable 
we get 

1 ikr 

Gir) = --— . (31) 
47r r 

This function is not determined univoquely since the Green function can 
be any solution of the eq. The right particular solution is chosen by 
imposing boundary conditions on the eigenfunctions "0/0 (r)- 

The Green function obtained in this way is 

G("') = -(j^) . (32) 

Thus, we finally get the integral equation for the scattering wave function 



tl^ik, r) = ip{k, r) - — 2 / rUir'Mk, r)dr , (33) 

27rn, J r — r 

where (/? is a solution of the Helmholtz equation. Noticing that |r — r'| = R, 
then 

(V^ + k'^)i; = (V^ + k^)[<f + J G{r,r')U{r')'il^{r')d^r'] (34) 

and assuming that we can change the order of operations and put the V 
operator inside the integral, we get 

(V^ + k^)ij = y"(V^ + k^)G{r, r')Uir')ij{r')d^r' = C/(r)V'(r) , (35) 

1 ikR 

which shows us that G{R) = -^^-^ is indeed a solution. 

Optical theorem 

The total cross section is given by 

-£dn . (36) 



See problem 7.1. 
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Let us express now f{6) as a function of the phase shift Si{k) = e^'^^'-^^^ 
in the form 

1 oo 

/(^) = r E(2^ + sin(5Kfc)PKcos0) . (37) 



Then 



/I OC) 

1=0 

[ y + ^^'^'"''^''^ sin<5p(A;)PHcos0)] . (38) 



/'=0 

47r 



Using now / P/(cos 0)P;/(cos 0) = 2^r[Sii', we get 



oo 



(Ttot = ^Y.i2l + l)sin5i{kf . (39) 



Of interest is the relationship 

1 oo 1 oo 

Im/(0) = -^(2Z + l)Im[e^^'«sin<5KA;)]P,(l) = - ^(2^ + 1) sin .JKfc)^ = 

1=0 1=0 

i-^ot , (40) 

which is known as the optical theorem. Its physical significance is related to 
the fact that the interference of the incident wave with the dispersed wave 
at zero/forward angle produces the "getting out" of the particle from the 
incident wave, allowing in this way the conservation of the probability. 



Born approximation 

Let us consider the situation of Fig. 7.2: 
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Fig. 7.2 



The observation point M is far away from P, which is in the range of the 
potential U. The geometrical conditions are r ^ L, r' ^ /. The segment MP 
that corresponds to |r — r'| is in the aforementioned geometrical conditions 
approxiamtely equal to the projection of MP onto MO 

|r — r'l ~ r — u • r' , (41) 

where u is a unit vector (versor) in the r direction. Then, for large r 

1 „ife|r-r'| 1 ikr 

G = -^^ J. c^r-^oo e-^"- . (42) 

47r |r — r'l 47r r 

We now substitute G in the integral expression for the scattering wave func- 
tion ^ 

V,(r) = e'''' - / e-^'="-'C/(r')V'(r')rfV' . (43) 

The latter is already not a function of the distance r = OM, but only of 9 
and ip. Thus 

/(0,^) = -i- J e-^*'" ■''U{r')^{r')(fr' . (44) 

We define now the incident wave vector ki as a vector of modulus k directed 
along the polar axis of the beam. Then e*'^^ = e*'*' '". Similarly, kj, of 
modulus k and of direction fixed by 6 and (p, is called the shifted wave 
vector in the direction {0,ip): kd = ku. 

The momentum transfer in the direction {6, ip) is introduced as the vec- 
torial difference K = kd — ki. 
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Fig. 7.3 

Hence we can write the integral equation in the form 

V;(r) = e^'^--'" + J G{r,r')U{r')i;{r')d\' (45) 

One can try to solve this equation iteratively. Putting r r';r' r" , 
we can write 

^(r') = e'^' "' + J G{r' y)U{v")^{T")(fr" . (46) 



Substituting in 45, we get 



^(r) = e^'^'-^ + j G(r,r')C/(r')e*'^'-'''dV+ 

j J G{r,r')U{r')G{r',r")U{r")i;{r")d^r"d^r' . (47) 

The first two terms in the right hand side are known and it is only 
the third one that includes the unknown function ipir). We can repeat the 
procedure: substituting r by r", and r' by r'", we get ip{r") , that we can 



reintroduce in the eq. 47 



^^r) = e'^'-" + J G(r,r')?7(r')e^'^'-"' + 

J J G{r,r')U{r')G{r' ,r")U{r'Y^'-''" d:^r'd^r"+ 

G(r, r')f/(r )G(r , r")U{r")e'^'-''" G{r" ,r"')U{r"')i^{r"') . (48) 

The first three terms are now known and the unknown function VCr) has 
been sent to the fourth term. In this way, by succesive iterations we can 
build the stationary dispersed wave function. Notice that each term of the 
series expansion has one more power in the potential with respect to the 
previous one. We can go on until we get a negligible expression in the right 
hand side, obtaining 'ip{r) as a function of only known quantities. 
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Substituting the expression of ^p{r) in f(9,ip), we get the expansion in 
Born series of the scattering amphtude. In first order in U, one should 
replace tpi'"^') by e*^' "" in the right hand side to get 

Air J Att J 

e-^^""'[/(r')dV (49) 

K is the momentum transfer vector. Thus, the differential cross section is 
simply related to the potential, V{r) = ^U{r). Since a{9,(p) = \ f{6,(p)\'^, 
the result is 

^T^'^HO,^) = ^1 / e-K-y(r)dVp (50) 

The direction and modulus of K depends on the modulus k of kj and kd as 
well as on the scattering direction (9, </?). For given 9 and it is a function 
of k, the energy of the incident beam. Analogously, for a given energy, a^^^ 
is a function of 9 and if. Born's approximation allows one to get information 
on the potential V{r) from the dependence of the differential cross section 
on the scattering direction and the incident energy. 

7N. Note - The following paper of Born was practically the first dealing 

with quantum scattering: 

M. Born, "Quantenmechanik der Stossvorgange" ["Quantum mechanics of 
scattering processes "], Zf. f. Physik 37, 863-867 (1926) 




7P. Problems 

Problem 7.1 

Calculus of complex variable for the scattering Green function 

We recall that we already obtained the result 

G{r, r') = -^2^ / fcl'I^a d'^Q , cu i? = |r — r'|. Since cfq = sin 9 dqdOdcp, 

we get after integrating in angular variables 
fjr^ „/^ _ i roo (e-'''«-e''^fi) , 

^y^' ^ ) - J-oo — — ^otg • 
Putting C = 4^p^, we separate the integral in two parts 

Let us make now q —qm. the first integral 
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i-oo fc^-(-q)^ i-Qjdi-Q) = Joo fe^^t^^ = - J-oo k^ldQ , 

so that 

G(r,r') = -2C(/_- 
Substituting C, leads to 

In this form, the integral can be calculated by means of the theorem of 
residues of its poles. Notice the presence of simple poles at q =i k. 





Fig. 7.4: Contour rules around the poles for G+ and G_ 

We use the contour of fig. 7.4 encircling the poles as shown, because in 
this way we get the physically correct effect from the theorem of residues 
G(r) = -^^ (ImA;>0), 
GW = -^^ (ImA:<0). 

The solution of interest is the first one, because it provides divergent 
waves, whereas the latter solution holds for convergent waves (propagating 
towards the target). Moreover, the linear combination 

i limg^o[Gfc+ie + Gk-ie] = —^^^7^ 
corresponds to stationary waves. 

The formal calculation of the integral can be performed by taking — q^ ^ 
k"^ + ie- q^, so that: w^^q j^r^^^^dq . This is possible 

for > 0. This is why the contour for the calculation will be placed 
in the upper half plane. Thus, the poles of the integrand are located at 
q = ±Vk'^ + ie ~ ±{k + ^). The procedure of taking the limit e ^ should 
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be applied after calculating the integral. 



Problem 7.2 
Asymptotic form of the radial function 

As we have already seen in the chapter Hydrogen atom the radial part 
of the Schrodinger equation can be written 

ii. + 'Ti)Rnlmir) - |?[F(r) + iigi^]i?„,^(r) + ^-^R^Ur) = . 
n,l,m are the spherical quantum numbers. For the sake of convenience of 
writing we shall discard them hereafter. R is the radial wave function (i.e., 
depends only on r). We assume that the potential goes to zero stronger 
than 1/r, and that lim,r_+o ^^^(^) = 0. 

Using u(r) = rR, since + f = ^J^n, we have 

i.u+p[E-Vir)-^-^]u = 0. 
Notice that the potential displays a supplementary term 

V(r) V{r) + , 
which corresponds to a repulsive centrifugal barrier. For a free particle 
V{r) = 0, and the equation becomes 

[i.+^^i)-'-^]R+eR=o. 

Introducing the variable p = kr, we get 

dp^ p dp 

The solutions are the so-called spherical Bessel functions. The regular solu- 
tion is 

while the irregular one 

niip) = -{-P)\\i)\^) . 

For large /?, the functions of interest are the spherical Hankel functions 
h^'\p)=ji{p)+ini{p) §i hf\p) = [h^^\p)r . 
The behaviour for p^ I is oi special interest 

-sin(p- ^) (51) 
P 2 

1 l-K 

Mp)- — cos(p-— ). (52) 
P 2 

Then 

~ _ie»(p-W2) 
I p 

The solution regular at the origin is Ri{r) = ji{kr) . 
The asymptotic form is (using eq. |5^) 
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Problem 7.3 

Born approximation for Yukawa potentials 

Let us consider the potential of the form 

— av 

V{v) = , (53) 

where Vq and a are real constants and a is positive. The potential is ei- 
ther attractive or repulsive depending on the sign of Vb; the larger |Vb|, the 
stronger the potential. Wc assume that |Vb| is sufficiently small that Born's 
approximation holds. According to a previous formula, the scattering am- 
plitude is given by 

Since this potential depends only on r, the angular integrals are trivial 
leading to the form 

= ^2^^/-sin|K|r^rdr . 
Thus, wc obtain 

f''KeM = -'-^^- 

From the figure we can notice that |K| = 2fesin |. Therefore 
^(B)(0) = l!!g^ . 

The total cross section is obtained by integrating 
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8. PARTIAL WAVES 



Introduction 

The partial waves method is quite general and applies to particles interact- 
ing in very small spatial regions with another one, which is usually known 
as scattering center because of its physical characteristics, (for example, 
because it can be considered as fixed). Beyond the interaction region, the 
interaction between the two particles is usually negligible. Under this cir- 
cumstances, it is possible to describe the scattered particle by means of the 
Hamiltonian 

H = Ho + V , (1) 

where Hq corresponds to the free particle Hamiltonian. Our problem is to 
solve the equation 

{Ho + V)\^p)=E\^) . (2) 

Obviously, the spectrum will be continuous since we study the case of 
elastic scattering. The solution will be 

I ^) = I ^)+ I ^) ■ (3) 

It is easy to see that for F = one can obtain the solution | (p), i.e., 
the solution corresponding to the free particle. It is worth noting that in a 
certain sense the operator ^-Hq '^^ anomalous, because it has a continuum 
of poles on the real axis at positions coinciding with the eigenvalues of Hq. 
To get out of this trouble, it is common to produce a small shift in the 
imaginary direction (±ie) of the cut on the real axis 

This equation is known as the Lippmann-Schwinger equation. Finally, 
the shift of the poles is performed in the positive sense of the imaginary 
axis because in this case the causality principle holds (cf. Feynman). Let us 
consider the x representation 

(x I V±) = (x I 0) + I d^x (x| ^_^^^.J x^ (x I y I . (5) 

The first term on the right hand side corresponds to a free particle, 
while the second one is interpreted as a spherical wave getting out from the 
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scattering center. The kernel of the previous integral can be considered as 
a Green function (also called propagator in quantum mechanics). It is a 
simple matter to calculate it 

G.(x,x') = ^l/x| |x\ = -i-f""'T','. (6) 

where E = fi^k"^ /2m. Writing the wave function as a plane wave plus a 
divergent spherical one (up to a constant factor), 



(x| V+) = e^-^ + ^/(k,k') . (7) 

the quantity /(k, k ) is known as the scattering amplitude and is explicitly 

1 Im 

/(k,k') = --(27r)^^(k'|y|V'+). (8) 

Let us now define an operator T such that 

T\ct>) = V\^l;+) (9) 

If we multiply the Lippmann-Schwinger equation by V and make use of 
the previous definition, we get 

r|0) = y|0) + F— — ^— -T|0) . (10) 

Iterating this equation (as in perturbation theory) we can get the Born 
approximation and its higher-order corrections. 



Partial waves method 

Let us now consider the case of a central potential. In this case, using the 
definition (9) , it is found that the operator T commutes with and L; it is 
said that T is a scalar operator. To simplify the calculations it is convenient 

to use spherical coordinates, because of the symmetry of the problem that 
turns the T operator diagonal. Let us see now a more explicit form of the 
scattering amplitude 

/(k,k') = const. J2 j j dE {\l \ E l' m ) {E l' m \ T \ Elm){Elm \ k) , 

Iml'm' 

(11) 
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where const. 



47r h- 

/(k,k' 



1 2m^27r)3. After some calculation, one gets 



(12) 



Choosing the coordinate system such that the vector k have the same 
direction with the z axis, one infers that only the spherical harmonics of 
m = will contribute to the scattering amplitude. If we define by 9 the 
angle between k and k , we will get 



>l°(k') 



121 + I 
47r 



PiicosO) . 



Employing the following definition 

fi{k) ^ 

eq. (12) can be written as follows 



irTiiE) 



/(k,k') = f{e) = J^{2i + i)fi{k)Pi{cose) . 



(13) 



(14) 



(15) 



For fi{k) a simple interpretation can be provided, which is based on 
the expansion of a plane wave in spherical waves. Thus, we can write the 
function (x | ^"'") for large values of r in the following form 



(x I V'^) 



(27r)3/2 



Y,{2l + l)Pi{cose) 



1 



(27r)3/2 



' „ikr pi{kr—liT)\ „ikr 

-]+J2i'^i + i)fiik)Pi{cose) — 
J I ^ . 



2ikr 



(2^)3/2 Z.i^' + ^) 2ik 



pikr pi{kr—liT) 

[l + 2ikfi{k)] 

r r 



(16) 



This expression can be interpreted as follows. The two exponential terms 

correspond to spherical waves: the first to a divergent wave, and the latter to 
a convergent one. Moreover, the scattering effect is conveniently displayed 
in the coefficient of the divergent wave, which is unity when there are no 
scattering centers. 
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Phase shifts 



We consider now a surface enclosing the scattering center. Assuming that 
there is no creation and annihilation of particles, one has 



/ 



j • = , (17) 

where the integration region is the aforementioned surface, and j is the 
probability current density. Moreover, because of the conservation of the 
orbital momentum, the latter equation should hold for each partial wave. 
The theoretical formulation of the problem does not change if one assumes 
the wave packet as a flux of noninteracting particles propagating through a 
region of central potential for which the angular momentum of each particle 
is conserved, so that the 'particle' content of the wave packet really does not 
change. Thus, one may think even intuitevely that only phase factor effects 
can be introduced under these circumstances. Thus, if one defines 

Si{k) = l + 2ikfi{k) (18) 

we should have 

I Si{k) 1= 1 . (19) 

These results can be interpreted using the conservation of probabilities. 
They are natural and expected because wc assumed that there is no creation 
and annihilation of particles. Therefore, the effects of the scattering center 
is reduced to adding a phase factor in the components of the divergent wave. 
Taking into account the unitarity of the phase factor, we can write it in the 
form 

Si = e^'^' , (20) 

where Si is a real function of k. Taking into account the definition (18), we 
can write 

■'^ 2ik k kcot{5i)-ik ' ^ ^ 

The total cross section has the following form 

(Ttotal = J I f{0)\^dn = 

^J^"d<P fd{cos{9))J2Y.('^^ + + ^y^' sm{6i)e'^i' sin(5^0^i^;' 



I I 
An 



^(2/ + 1) sin 2(5,0 • (22) 
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Getting the phase shifts 

Let us consider now a potential V that is zero for r > R, where the parameter 
R is known as the range of the potential. Thus, the region r > R corresponds 
to a spherical unperturbed/free wave. On the other hand, the general form 
of the expansion of a plane wave in spherical ones is 

{^\^^) = -^:^Y.''{'^l + m{r)Pi{cos0) {r>R), (23) 

where the coefficient Ai is by definition 

Ai = hf^ (kr) + hf^ (kr) . (24) 
h\^^ and hf'^ are the spherical Hankel functions whose asymptotic forms 



are the following 



.(1) 



gj(fer--/7r/2) 

ikr 

g-i(A;r-«7r/2) 

ikr 



Inspecting the following asymptotic form of the expression (23) 



1 



(2vr)3/2 ^ 
one can see that 



J,kr 



2ikr 



.(2) 



,—i{kr—l-K) 

2ikr 

1 

2 ■ 



This allows to write the radial wave function for r > i? in the form 
Ai = e^*"^' [cos 6iji{kr) — sin 5ini{kr)] . 



(25) 

(26) 
(27) 



Using the latter equation, we can get the logarithmic derivative inr = R, 
i.e., at the boundary of the potential range 



r dA 



A, dr 



kR 



ji cos 6i — Ui {kR) sin 5i 
ii cos (5/ - ni{kR) sin (5; 



(28) 



is the derivative of ji with respect to r evaluated at r = i?. Another 
important result that can be obtained from the knowledge of the previous 
one is the phase shift 



tan (J; = 



kRj'iikR) - pijijkR) 
kRn'i{kR) - /3ini{kR) 



(29) 
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To get the complete solution of the problem in this case, it is necessary 
to make the calculations for r < R, i.e., within the range of the potential. 
For a central potential, the 3D Schrodinger equation reads 

(Pui /,2 2m,, + . . 

-3^+{'^-^V-±^)u, = 0. (30) 

where ui = rAi{r) is constrained by the boundary condition ui \r=o = 
0. Thus, one can calculate the logarithmic derivative, which, taking into 
account the continuity of the log-derivative (equivalent to the continuity 
condition of the derivative at a discontinuity point) leads to 

A \in= A lout ■ (31) 



An example: scattering on a hard sphere 

Let us now consider an important illustrative case, that of the hard sphere 
potential 

y=in (32) 
] r > R . ^ ^ 

It is known that a particle cannot penetrate into a region where the 
potential is infinite. Therefore, the wave function should be zero at r = i?. 
Since we deal with an impenetrable sphere we also have 

Ai{r) \r=R= . (33) 

Thus, from eq. (27), we get 

ni{kR) 

One can see that the phase shift calculation is an easy one for any I. In 
the I = case (s wave scattering), we have 

6i = -kR 

and from eq. (27) 

. , , sinkr ^ cos kr . ^ 1 • /, c- n /o^x 
Ai=o[r) ~ cosOo H smOo = ^ sm(fcr + Oq) . (35) 
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We immediately see that tliere is an additional phase contribution with 
regard to the motion of the free particle. It is also clear that in more general 
cases the various waves will have different phase shifts leading to a transient 
distortion of the scattered wave packet. At small energies, i.e., kR « 1, 
the spherical Bessel functions (entering the formulas for the spherical Hankel 
functions) are the following 

jlikr) ~ , [^''^\„ (36) 
^ ' (2Z + 1)!! ^ ' 



(2^-1)!! 
{kr) 
leading to 



Mkr) ~ -hlZMTr ' (37) 



= (2/ + l)[(2Z-l)!!p • 

From this formula, one can see that a substantial contribution to the 
phase shift is given by the / = waves. Moreover, since Sq = —kR the cross 
section is obtained as follows 



total 



[ = AttR^ . (39) 



One can see that the total scattering cross section is four times bigger 
than the classical one and coincides with the total area of the impenetrable 

sphere. For large values of the incident energy, one can work in the hypoth- 
esis that all values of Z up to a maximum value Imax ~ kR contribute to the 
total cross section 

<Ttotal = T2 E + • (40) 

1=0 

In this way, from eq. (34), we have 

9 c- tan2(5, IjiikR)? . 9 „ ^7r\ , 

sin^Si = ^ = , — „ ~ sin^ I kR , (41 

l + tan2,5, \ji{kR)]^ + ImikR)]^ \ 2) ' ^ ' 

where the expressions 



jlikr) ~ — sin (kr — — 
■^'^ ' kr \ 2 

nAkr) ~ — — cos (kr — . 

kr V 2 / 
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have been used. 

Inspection of 5i shows a negative jump of | whenever I is augmented by 
a unity. Thus, it is clear that sin^(5; + sin^(5i+i = 1 holds. Approximating 
siv?5i by its mean value ^ over a period and using the sum of odd numbers, 
one gets 

<y total = %ikRf\ = 2^^' • (42) 

Once again the quantum- mechanical result, although quite similar to the 
corresponding classical result is nevertheless different. What might be the 
origin of the factor of two that makes the difference ? To get an explanation, 
we first separate eq. (15) in two parts 

, l=kR ■ l=kR 

m = ^ E m+^)e'''^Pi cos(^)+^ E (2^+1)^^ cos(0) = /,efi+/shadow ■ 

(43) 

Calculation of / \ f j-efiP^^ gives 

/<2 Imax j. 1 l2 
l/refll'rf^ = i^E {2l + lf[Picos{e)fd{cose) = ^ = nR^ . 

(44) 

Analysing now /g^adow small angles, we get 

/shadow - i E(2^ + ^)Jom - ik j\jo{kb9)db = 'IIl^ . (45) 

This formula is rather well known in optics. It corresponds to the Praun- 
hofer diffraction. Employing the change of variable z = kRd one can calcu- 
late the integral / |/ shadow l^'^^ 

/ l/shadowl''^^ ~ '^''^^ r^^^"^' ~ ■ (^^^ 

Finally, neglecting the interference between /j.g£ and / g^adow (^i^^^e 
the phase oscillates between 25;_|_i = 25i — tt), one gets the result (42). The 
label 'shadow' for one of the terms is easily explained if one thinks of the 
wavy behaviour of the scattered particle (from the physical viewpoint there 
is no difference between a wave packet and a particle in this case). Its origin 
can be traced back to the backward-scattered components of the wave packet 
leading to a phase shift with respect to the incident waves and destructive 
interference. 
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Coulomb scattering 



In this section we briefly consider the Coulomb scattering in the quantum- 
mechanical approach. For this case, the Schrodinger equation is 

(-|^V^-^^)v'(r) = i?^(r), ^>0, (47) 

where m is the reduced mass of the system, E > since we deal with the 
simple scattering case where no kind of bound states are allowed to form. 
The previous equation is equivalent to the following expression (for adequate 
values of the constants k and 7) 

V2 + A;2 + ^)V'(r)=0. (48) 

If we do not consider the centrifugal barrier, i.e., we look only to the s 
waves, we really deal with a pure coulombian interaction, for which one can 
propose a solution of the following form 

^(r) = e^'^-'-x(tx) , (49) 

where 

u = ikr{l — cos 6) = ik{r — z) = ikw , 
\s. - r = kz . 

V'(r) is the complete solution of the Schrodinger equation with an asymptotic 

'physical' behaviour to which a plane wave e^^'^ and a spherical wave are 
expected to contribute r~^e**^^ are expected to contribute. Defining new 
variables 

z = z w = r — z \ = 6 , 



and by employing of previous relationships, eq. (48) takes the form 

X{u) = . (50) 



To solve this equation, one should first study its asymptotic behaviour. 
Since we have already tackled this issue, we merely present the asymptotic 
normalized wave function that is the final result of all previous calculations 

V'k(r) = U^-r-,lnikr-^..r)] + fc{k , 9)e^l^^+'r^-^''^^ \ ^^^^ 
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As one can see, this wave function displays terms that turns it quite dif- 
ferent from the form in eq. (7). This is due to the fact that the Coulomb po- 
tential is of infinite range. Performing the exact calculation for the Coulomb 
scattering amplitude is not an easy matter. Here we give only the final result 
for the 'normalized' wave function 

1 / Ti 111.^ M Otf^) T gj[fcr+7/n2fcr] \ 



(27r)3/2 \^ 5i(7)2A:sin(e/2)2 r j ' 

(52) 

where 51(7) = r(iij^) ■ 

In addition, we reduce the partial wave analysis to a clear cut presenta- 
tion of the results, of which some have already been mentioned. First of all, 
we write the wave function V'(r) in (49) as follows 

^(r) = e^^-^x(^^) = [ e"*fT-Hl - ty'^^dt , (53) 

Jc 

where A is a 'normalization' constant, while all the integral part is the 
inverse Laplace transform of the direct transform of eq. (50). A convenient 
form of the latter equation is 

V'(r) = A [ 6**^-^(1 - t)e'''''\l - t)d{t, -f)dt (54) 
JC 

where 

d(t,7)=t*''-'(l-t)-'^-' . (55) 
Within the partial wave analysis we proceed by writing 

00 

V^(r) = ^(2/ + l)i^Pi{cose)Ai{kr) , (56) 

1=0 

where 

Ai{kr) = a[ e^'^'jilkril - i)](l - t)d{t, 7) . (57) 

Applying the relationships between the spherical Bessel functions and 
the Hankel functions, we get 

Ai{kr) = (kr) + ^ (kr) . (58) 

We shall not sketch here how these coefficients are obtained (this is quite 
messy). They are 

A\^\kr) = (59) 
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^''"^ (60) 

where 

Calculation of the Coulomb scattering amplitude 

If we perform the Laplace transform of eq. (50), we get 

XM =A I e"*f^-^(l - t)-'^dt . (62) 
Jc 

The contour C goes from — oo to oo on the real axis and closes through 
the upper half-plane. There are two poles in this case at f = and t = 1. 
By the change of variable s = ut, we get 



X(n) = A f e's''^-\u - s)'"'^ . (63) 



xiu) should be regular in zero. Indeed, we get 

;^(0) = (-1)-'^A I —ds . = (-\)-'^A2m (64) 

Performing now the limit u — oo, let's do an infinitesimal shift to avoid 
the location of the poles on the contour. Moreover, by the change of variable 
u ~ "■^^^^J i^isX this expression goes to zero when u oo. Thus, 

we can expand {u — s) in power series of ^ for the pole with s = 0. This 
expansion is not the right one in s = 1, because in this case s = — so+i(«;±£). 
It comes out that ^ = 1— ^^"^^^•^ tends to 1 when k ^ oo. If instead we do 
the change of variable s' = s — u, we get rid of this difficulty 

X{u) = (^[e"s*T-i(« - s)-^^]ds + [e"'+"(-s')^^(« + s')'^-^]ds'^ . 

(65) 

Expanding the power series, it is easy to calculate the previous integrals, 
but one should take the limit - ^ in the result in order to get the correct 
asymptotic forms for the Coulomb scattering 

X{u) ~ 21,1 A [n-^V(7) - (-n)^^-^e"52(7) 
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27r32(7) = i I e's-'^ds . (66) 

After all this chain of variable changes, we get back to the original s one 
to obtain 

{u*y^ = {-iy^[k{r - z)p = eW2ei7lnfc{r-z) 

The calculation of x, once effected, is equivalent with having V'k(r) start- 
ing from (49). 

Eikonal approximation 

We shall briefly expound on the eikonal approximation whose philosophy 
is the same to that used when one wants to pass from the wave optics to 
the geometrical optics. Therefore, it is the right approximation when the 
potential varies slowly over distances comparable to to the wavelength of 
the scattered wave packet, i.e., for the case E » \V\. Thus, this approx- 
imation may be considered as a quasiclassical one. First, we propose that 
the quasiclassical wave function has the known form 

^ ^ e^s^^)!^ , (68) 
where S satisfies the Hamilton-Jacobi equation, having the solution 



n J-00 I ^ ' 



1/2 

dz' -I- constant . (69) 



The additive constant is chosen in such a way to fulfill 

%^kz for y ^ . (70) 

n 

The term multiplying the potential can be interpreted as a change of 
phase of of the wave packet, having the following explicit form 

Within the method of partial waves, the eikonal approximation has the 
following application. We know it is correct at high energies, where many 
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partial waves do contribute to the scattering. Thus, we can consider I as 
a continuous variable and by analogy to classical mechanics we let I = bk. 
Moreover, as we already mentioned Imax = kR, which plugged into eq. (15) 
leads to 

f{e) = -ik j bJo{kbe)[e^''^^^^ -l]db. (72) 



8P. Problems 

Problem 8.1 

Obtain the phase shift and the differential cross section at small angles for 
a scattering centre of potential U{r) = It should be taken into account 
that for low-angle scattering the main contribution is given by the partial 
waves of large I. 

Solution: 

Solving the equation 



R, + 



k' 



l{l + 1) 2ma 







with the boundary conditions Ri{0) = 0, -R/(oo) = A'^, where A'^ is a finite 
number, we get 

Ri{r) = A^hikr) , 

and I is the first modified Bessel function. 



where A = 



(1 + 



1\2 _j_ 2ma 



To determine 5i, one should use the asymptotic expression of ly. 



Therefore 
Si = - 



X-l 



TT 

2 



(1 + 



+ 



2ma 



1/2 
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The condition of large I leads us to 

nma 



Si 



{21 + i)n 



2 ' 



whence one can see that \5i\ <^ 1 for large /. 

From the general expression of the scattering amplitude 



-. oo 

m = —J2{2l + l)mcos9){e''''-l), 



.=0 

at small angles one gets e^"^' 1 + 2iSi , so that 

oo 1 

^PKcos0) = — ^ 
J^o 2sm| 

Thus 

' kh'' 2 sin I 

The final result is 

da _ ■n^a^m 
d0 ~ ~We''^^2 ■ 
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